Linda Leone,Gianmattia Sgueglia,Salvatore La Gatta,Marco Chino,Flavia Nastri,Angela Lombardi
Enzymatic and Bioinspired Systems for Hydrogen Production
hydrogen,hydrogenases,hydrogen-evolving catalysts,bioinspired catalysis,artificial enzymes,small-molecule complexes,sustainable energy production,biofuel cells,photocatalytic hydrogen
The extraordinary potential of hydrogen as a clean and sustainable fuel has sparked the interest of the scientific community to find environmentally friendly methods for its production. Biological catalysts are the most attractive solution, as they usually operate under mild conditions and do not produce carbon-containing byproducts. Hydrogenases promote reversible proton reduction to hydrogen in a variety of anoxic bacteria and algae, displaying unparallel catalytic performances. Attempts to use these sophisticated enzymes in scalable hydrogen production have been hampered by limitations associated with their production and stability. Inspired by nature, significant efforts have been made in the development of artificial systems able to promote the hydrogen evolution reaction, via either electrochemical or light-driven catalysis. Starting from small-molecule coordination compounds, peptide- and protein-based architectures have been constructed around the catalytic center with the aim of reproducing hydrogenase function into robust, efficient, and cost-effective catalysts. In this review, we first provide an overview of the structural and functional properties of hydrogenases, along with their integration in devices for hydrogen and energy production. Then, we describe the most recent advances in the development of homogeneous hydrogen evolution catalysts envisioned to mimic hydrogenases.
Enzymatic and Bioinspired Systems for Hydrogen Production

The extraordinary potential of hydrogen as a clean and sustainable fuel has sparked the interest of the scientific community to find environmentally friendly methods for its production. Biological catalysts are the most attractive solution, as they usually operate under mild conditions and do not produce carbon-containing byproducts. Hydrogenases promote reversible proton reduction to hydrogen in a variety of anoxic bacteria and algae, displaying unparallel catalytic performances. Attempts to use these sophisticated enzymes in scalable hydrogen production have been hampered by limitations associated with their production and stability. Inspired by nature, significant efforts have been made in the development of artificial systems able to promote the hydrogen evolution reaction, via either electrochemical or light-driven catalysis. Starting from small-molecule coordination compounds, peptide- and protein-based architectures have been constructed around the catalytic center with the aim of reproducing hydrogenase function into robust, efficient, and cost-effective catalysts. In this review, we first provide an overview of the structural and functional properties of hydrogenases, along with their integration in devices for hydrogen and energy production. Then, we describe the most recent advances in the development of homogeneous hydrogen evolution catalysts envisioned to mimic hydrogenases.

Over the last century, modern society has reached a peak of energy consumption, required to withstand the technological progress and population growth. Fossil fuel is currently the primary source of energy, despite its impact on the environment and limited supply. Hence, developing sustainable alternatives is definitely the most urgent challenge for humanity [1]. Hydrogen (H2) represents a promising fuel to manage the existing energetic crisis, owing to its exceedingly high energy density (~120 MJ/kg), which surpasses carbon-based sources as methane and natural gas (~50 MJ/kg) [2]. Most importantly, hydrogen can be converted into electricity in a fuel cell or combusted into an engine, producing only water as the byproduct. Despite no carbon products being formed through its exploitation, considering hydrogen as a clean fuel is not straightforward, since several issues need to be addressed concerning the environmental impact of its production, transportation, and storage [3]. Being gaseous in ambient conditions, hydrogen needs to be stored under high pressure or liquified. These processes require additional energy, which may not come from renewable sources. Furthermore, despite being the most abundant element on Earth, hydrogen is readily available only in combination with other elements, but not in the elementary form. Currently, about ~96% of H2 is produced on an industrial scale from nonrenewable fossil fuels, mostly via the steam-reforming process, where water steam reacts with hydrocarbons to give H2 and CO2 [4]. The significant emissions of greenhouse gases, coupled to excessive energy requirement, are some of the inherent limitations of these methods. Therefore, there is an urgent need to design H2 production systems that are both effective and sustainable. Electrolysis of water can be considered as a sustainable alternative, especially when using solar or wind energy as source of electric power [5,6]. However, this method is still not economically viable, due to the high energy consumption and slow hydrogen evolution rates. To increase the efficiency of such a process, much research has been devoted to developing proficient and low-cost electrocatalysts for hydrogen evolution. Biological catalysts are the most environmentally friendly alternatives, as they usually operate under mild reaction conditions and typically do not produce carbon-containing byproducts. Biological pathways include those exploiting wastewater or organic solid waste as substrates [7,8]. Among them, anaerobic fermentation is the most studied, involving the conversion of organic molecules from waste by microorganisms such as anoxic bacteria [9,10,11]. Photobiological processes are, in principle, even more advantageous, as they involve the use of photosynthesis and the natural water-splitting process, in which the microbes absorb energy from the sun and evolve H2 directly from water, without an intermediate biomass stage [12,13]. Such a process requires the synergic work of a photosynthetic system and a hydrogenase enzyme, which are both found in cyanobacteria and green algae [14,15]. Albeit promising from a sustainability perspective, photobiological H2 production is limited by the low overall efficiency of natural photosynthesis. The challenge of converting solar energy more efficiently than natural organisms has inspired the development of hybrid systems comprising the use of synthetic photosensitizers coupled to hydrogenase enzymes [16,17,18,19]. However, scaling up these systems has proven to be very difficult due to the need for large-scale expression and purification of the enzymes. Furthermore, a significant drawback is the high O2 sensitivity of hydrogenases, particularly algal [FeFe]-hydrogenases [20]. A possible strategy to overcome the limitations associated with the low availability of hydrogenases is to build entirely artificial systems, which require the construction of robust, efficient, and easily viable mimetics of hydrogenase enzymes [21]. Following this path, a wide variety of molecular and biomolecular catalysts have been developed over the last decades, displaying either electrochemical or light-driven hydrogen-evolving activity [22]. In this review, we first describe the structural and functional properties of hydrogenases and summarize the most recent advances in the development of devices for hydrogen and energy production based on these natural enzymes. Then, we provide an overview of the artificial homogeneous catalysts envisioned to mimic hydrogenase function, including small-molecule, peptide-based, and protein-based analogues.

Hydrogenases are a class of widely distributed metalloenzymes able to catalyze both hydrogen production from protons and electrons and the reverse reaction, i.e., oxidation of molecular hydrogen to protons (H2 ⇌ 2H+ + 2e−) [23]. The relevance of this process to hydrogen synthesis has made them privileged biomolecular systems to study, with the initial goal of obtaining catalysts and devices that could be used for hydrogen evolution and, ultimately, for sustainable energy production. In biology, hydrogenases are required by some organisms, such as archaea, bacteria, and green algae, for H2-based respiration and are involved in maintaining the redox balance inside cells, by using reducing equivalents to reduce protons where necessary. Hydrogenases are usually located in the cytoplasm or periplasm, in either membrane-bound or soluble form. Several hydrogenases, known as regulatory hydrogenases, almost devoid of catalytic activity, are instead devoted to simply sensing the presence of hydrogen and precisely controlling the expression of genes coding for active hydrogenases and other associated proteins [24]. Most of these metalloenzymes are sensitive to the presence of molecular oxygen [25], which inactivate them via different mechanisms, depending on the specific hydrogenase type. This weakness is probably inherited, due to the origin of these enzymes in anaerobic organisms. Despite this, their overall structure and the nature of the active site cofactor have been finely tuned by evolution, resulting in impressively efficient catalysts with very high turnover rates that use exclusively Earth-abundant metals. It is interesting to note that hydrogenases display remarkable diversity in active site composition, which is the result of convergent evolution. Depending on the nature of metal ions housed within the active site, hydrogenases are classified into three distinct classes: [FeFe]-hydrogenases, [NiFe]-hydrogenases and [Fe]-hydrogenases (also called Hmd). Their peculiar features are briefly described below.

[FeFe]-hydrogenases (FFHs) represent the most active group of hydrogenases, found in both prokaryotes and eukaryotes [26]. They are involved in hydrogen metabolism of several microorganisms, ranging from green algae to sulfate-reducing bacteria and many others. Putative FFH sequences have been identified in the genome of many organisms [27], and various classifications have been proposed [28]. A structural feature that is common to all FFHs is the so-called H-domain, a conserved ~40 kDa domain hosting the active site, where hydrogen chemistry takes place. Catalytic activity requires a remarkably complex metal cofactor known as the H-cluster [29]. This metal cofactor is itself formed by two redox-coupled iron sites. One is an iron–sulfur cubane cluster, named [4Fe4S]H, similar to those found in some ferredoxins and other electron transfer proteins [30]. Other [4Fe4S] clusters are also present in FFHs and are probably involved in electron transfer. The other iron site, named [2Fe]H, is the cofactor where hydrogen evolution/splitting occurs. The [2Fe]H cluster (Figure 1a) is much more unusual in terms of both structure and iron coordination sphere composition, containing two iron ions connected by two bridging thiolate sulfur atoms coming from a 2-aza-propane-(1,3)-dithiolate (ADT) group and by a CO molecule. Additionally, another CO molecule and a CN− anion are coordinated to each iron atom. One of the two irons (called proximal iron Fep) is also bonded to a cysteinyl thiolate sulfur, which is shared with the [4Fe4S]H cofactor, making it 6-coordinated, while the other (distal iron or Fed), maintains one vacancy to coordinate hydrogen. Although some of them have been found to be more resilient, FFHs tend to be very O2-sensitive and are usually destroyed or inactivated even by trace amounts of oxygen [31,32] unless protecting conditions or reagents are employed [33,34]. The catalytic cycle of FFHs is quite complex and involves several intermediates, some of which have been extensively characterized by spectroscopic techniques; however, some aspects of the cycle remain to be clarified. The description of the catalytic cycle is beyond the scope of this review; thus, readers can refer to exhaustive papers [35,36] for a current overview of the proposed catalytic intermediates and their characterization.

[NiFe]-hydrogenases (NFHs) are the most common type of hydrogenases [37], found in archaea, bacteria, and even eukaryotes. As the name suggests, their metal cofactor contains an iron and a nickel ion, but there are also other differences in metal coordination spheres compared to FFHs [38]. NFHs generally exhibit lower activity and reversibility compared to FFHs, and they are more commonly employed for hydrogen splitting rather than hydrogen evolution. Most NFHs characterized so far share a similar heterodimeric structure formed by a larger (named alpha, ~60 kDa) and a smaller (named beta, ~30 kDa) subunit. The larger subunit contains the active NiFe site, while the smaller subunit houses a variable number of [4Fe4S] clusters necessary for electron transfer to and from the NiFe site. NFHs are commonly classified into four groups [39,40] depending on differences in active site composition, location in the cellular space, and function. There exists also a subgroup of NFHs with high O2 tolerance, known as [NiFeSe]-hydrogenases (NFSHs), in which one of the cysteine residues bonded to the nickel is replaced by a selenocysteine [41,42]. NFHs are usually more tolerant to oxygen than FFHs, being merely reversibly deactivated by oxygen, while some even retain their catalytic activity in air. Interestingly, O2-tolerant NFHs often contain a unique metal cofactor, a [4Fe-3S] 6-Cys cluster, in place of the proximal [4Fe-4S] 4-Cys cluster, which seems to be implicated in enzyme recovery [43,44]. As mentioned, the metal cofactor responsible for catalytic activity differs significantly from that of FFHs [45]. In NFHs, the nickel atom is coordinated by four cysteine thiolate ligands, two of which are also coordinated to the iron atom. The latter is also coordinated by a CO molecule and two CN− anions, similarly to the FFH [2Fe]H cluster (Figure 1b). The catalytic cycle of NFHs is still under investigation; however, thanks to much experimental effort, many of the main events and intermediates have been determined [46,47,48]. The first NFH reaction mechanism suggested heterolytic cleavage of the H–H bond [49]. Nevertheless, a possible homolytic cleavage mechanism cannot be excluded, as evidenced by subsequent computational investigations, which resulted in comparable energy barriers for the two mechanisms [50].

[Fe]-hydrogenases (FHs, also known as H2-forming methylenetetrahydromethanopterin dehydrogenases or Hmds), are significantly rarer than FFHs and NFHs, having been found only in methanogenic archaea, particularly when nickel availability is decreased [51]. They catalyze the stereospecific reversible hydride transfer between H2 and N5,N10-methenyltetrahydromethanopterin, the latter being a necessary cofactor for the process. This reaction is part of the hydrogenotrophic methanogenesis pathway in methanogens, which extracts reducing equivalents from hydrogen to reduce CO2 to methane gas. From a structural viewpoint, FH is a homodimeric protein formed by two 37 kDa subunits, joined by a shared alpha-helical bundle formed by the C-terminal regions, while the N-terminal region adopts a structure derived from the Rossman fold [52]. The metal cofactor present in FHs is unique and greatly different from those present in FFHs and FNHs, since it does not contain an iron sulfur cluster or other polynuclear clusters. Indeed, FHs were initially thought to be metal-free hydrogenases. However, it is now accepted that FHs contain a metal cofactor constituted by an iron ion bound through a pyridine nitrogen and an unusual acyl-iron bond to guanylylpyridinol (GP), a molecule constituting a guanine monophosphate unit bonded to pyridinol, in addition to two CO molecules and a Cys thiolate sulfur (Figure 1c). This cofactor is both heat- and light-sensitive, rapidly decomposing into guanylylpyridinol, CO, and free iron [53], which makes its manipulation and study somewhat cumbersome, even though the enzyme is O2-tolerant. A mechanism has been proposed for FHs on the basis of both experimental investigations [54,55] and computational studies [56,57].

Scientific enquiry into the inner workings of hydrogenases is greatly motivated by the potentially useful application of these enzymes and of systems derived from them [58]. Since hydrogenases can catalyze both hydrogen splitting and hydrogen formation with high turnover and low overpotential, they are ideally suited for the construction of devices designed to produce hydrogen or exploit it as an energy source. To accomplish these tasks, hydrogenases should be incorporated or deposited on the electrodes of devices specifically designed to perform and electrically couple the required reactions [59]. The choice of electrode materials and incorporation of hydrogenases into the electrode are not trivial. However, several advances in this field have been recently accomplished by exploiting the properties of advanced materials [60,61,62,63,64]. NFSHs are usually preferred to other types of hydrogenases for practical applications, owing to their high activity combined with O2 tolerance. The setup and detailed nature of electrochemical devices depends on the specific process desired. Two of the major types of devices, photoelectrochemical cells and biofuel cells, are discussed in further detail below. We are aware that it is impossible to exhaustively cover the huge amount of work in the field. Thus, several case studies are selected to focus on the most recent and promising results on the incorporation of hydrogenases into these kinds of electrochemical devices.

Hydrogen production can be performed by devices in which proton reduction to hydrogen, catalyzed by hydrogenase deposited onto the cathode, is electrochemically coupled to oxidation of another chemical species. When that species is water, the whole process produces H2 and O2 and is named water splitting. Water splitting can also be fueled by energy supplied as light [65], in which case it effectively represents a kind of artificial photosynthesis that can be performed through specialized devices named photoelectrochemical cells (PECs, Figure 2a). Notably, these devices require conjugation of catalytic activity to light absorption, which is commonly accomplished through the use of proper materials, such as semiconductors [66], or by functionalizing either the catalyst or the electrode surface with a suitable dye [67]. In such devices, an oxygen-evolving catalyst (OEC) is bound to the photoanode, while the photocathode is functionalized with a hydrogenase enzyme or, in general, a hydrogen-evolving catalyst (HEC). Natural enzymes have been successfully immobilized on the surface of electrodes in various electrochemical devices, which are known as bio-photoelectrochemical cells (BPECs). These include a number of electrocatalytic biomolecules such as Photosystem I (PSI), Photosystem II (PSII) [68], bacterial photosynthetic reaction centers (RCs) [69], glucose oxidase (GO) [70], bilirubin oxidase (BO) [71], and laccase [72]. Even thylakoid membranes themselves have been incorporated in tailor-made devices [73]. Some recent and notable examples regarding the incorporation of hydrogenases in such devices are briefly described below. Sokol et al. [74] proposed a remarkable design for a PEC based on Photosystem II (PSII) and a [NiFeSe]-hydrogenase that successfully accomplishes unassisted water splitting. In their device, the photoanode is formed by PSII wired to dye-sensitized titanium dioxide, while the cathode is wired to a [NiFeSe]-hydrogenase. The dye employed was a green light-absorbing diketopyrrolopyrrole, while electron transfer between the dye and PSII was enabled by a redox polymer, poly(1-vinylimidazole-co-allylamine)-Os(bipy)2Cl. In a previously designed similar PEC [75], in which light was absorbed only by PSII, a large external potential was needed to sustain the process, because of the low electrochemical potential of electrons produced by the PSII. This issue was resolved by employing the green light-absorbing dye. TiO2 is one of the most commonly employed semiconductor materials for the electrodes in electrochemical devices and is also suited to the deposition of biocatalysts thanks to its porous surface, which is believed to facilitate the adsorption of proteins [76]. The efficiency of TiO2 light absorption can also be enhanced by manufacturing TiO2-based inverse opals (IO), materials with a highly organized internal arrangement giving rise to periodic modulation of the dielectric constant [77,78]. Nam et al. [79] developed a novel Si-based photocathode featuring a hierarchically structured inverse opal (IO)-TiO2 layer, optimized for integration of a [NiFeSe]-hydrogenase. The photocathode allowed stable and high hydrogenase load, and it was tested when coupled to both inorganic and biological photoanodes, based on BiVO4 and PSII, respectively, to perform the overall water-splitting process. Coupling to the BiVO4-based photocathode yielded a system capable of unassisted water splitting, while coupling to the PSII required application of an external voltage, albeit lower than reported by Sokol et al. Perovskites have also been extensively explored as semiconductor materials for various applications including BPECs. These are described by the general formula ABX3, where A and B are cations of different size, and X is oxygen or a halide [80]. Perovskite-based materials have been rapidly improving in the last decade and have enabled the creation of photovoltaic cells with impressive power conversion efficiencies, now surpassing 29% [81]. In this context, Moore et al. [82] designed a sophisticated photocathode based on an organic–inorganic lead halide perovskite. Despite their high efficiency, these materials are sensitive to heat, air, and moisture; therefore, their employment is quite challenging. The authors were able to protect the perovskite layer by encapsulation, using a eutectic alloy, metal foil, and epoxy resin, reporting for the first time the integration of an NFSH with this kind of material. By using a BiVO4-based photoanode, they were also able to accomplish bias-free water splitting with higher solar-to-hydrogen efficiency than those illustrated above [82]. Lastly, Tian et al. [83] developed a PEC using ZnO-protected Cu2O, a material never used before for a hydrogenase photocathode. Cu2O presents favorable properties as a p-type semiconductor, albeit unstable. To increase its stability, the authors employed a protective ZnO layer, which also allows electron flow from the Cu2O layer. An FFH was then immobilized on the ZnO layer, and the electrode was integrated into a testing PEC device, confirming that, when illuminated, the photocathode evolved hydrogen. The protecting ZnO layer showed signs of degradation after prolonged operation, resulting in the formation of nano-sticks.

Hydrogen oxidation to protons can be coupled to the reduction of other species in so-called fuel cells (FCs), which exploit the overall chemical reaction promoted by catalysts deposited onto the electrodes to generate electrical energy (Figure 2b) [84]. Conceptually, the counterpart of a PEC performing water splitting is an FC exploiting hydrogen oxidation to protons and oxygen reduction to water. When the catalysts used by the apparatus are enzymes, the device is known as a biofuel cell (BFC), in which hydrogenases are usually wired to the anode to catalyze the oxidation of hydrogen. A number of interesting hydrogenase-based BFC designs have been recently devised [58]. Gentil et al. [85] successfully immobilized an NFSH on carbon nanotube-modified electrodes using a rational design strategy. To this end, they grafted hydrophobic molecules to carbon nanotubes, such as anthraquinone and adamantane, selected to interact with the hydrogenase surface. Such modified nanotubes were able to bind the NFSH and were used to manufacture gas diffusion electrodes for a H2/air BFC operating at room temperature and pressure. In another study, Hardt et al. [86] designed a redox-active film composed of an FFH embedded into a 2,2′-viologen-modified hydrogel to perform bidirectional and reversible hydrogen conversion. This film was used to manufacture a bioanode for a H2/O2 BFC achieving a 1.16 V open-circuit voltage. The same film was also usable as a highly efficient biocathode for hydrogen evolution. The authors proposed a detailed kinetic model to explain the observed catalytic activity. Viologen-based polymers have been used extensively in BFC designs, and their properties can also be exploited to protect the sensitive hydrogenases from inactivation caused by oxygen, as proposed by Szczesny et al. [87]. In their recent work, the authors proposed a noteworthy design for a dual-gas breathing H2/air BFC, relying on a bilirubin oxidase-based biocathode and a NFSH-based bioanode. The bioanode is composed of two layers of a redox-active polymer which serves both as support and as protection for the NFSH from O2. The cell achieved an open-circuit voltage equal to 1.13 V. In a more recent study, the same authors were also able to design a H2/O2 BFC, comprising a hydrogenase/redox polymer-based bioanode integrating a highly active but O2-sensitive FFH, protected by the redox polymer [88]. An additional concern with hydrogenases is their thermal stability, which limits the temperature range in which BFCs can operate. To address this problem, Wang et al. [89] designed a novel bioanode using a thermostable and highly active NFH isolated in Pyrococcus furiosus. They first developed a thermophilic archaea host system to express the hydrogenase, which was then immobilized on a nanotube-modified carbon felt electrode. This bioanode was integrated in a H2/air BFC with a Pt/C cathode, successfully operating in the temperature range 40–80 °C with high efficiency. The practical feasibility of biohybrid devices such as BPECs and BFCs has been clearly demonstrated by these technically impressive designs. Integration of natural hydrogenases and other biologically derived electrocatalysts with electrochemical devices has been accomplished using a variety of different strategies. Despite the incredible advances reported so far, a number of difficulties remain to be addressed. Many of the materials used to make these devices are delicate and quite sophisticated, inflating the cost of their wide adoption. The overall efficiency of the setup is also of paramount importance in this regard. Lastly, biohybrid devices also pose specific problems arising from the chemical fragility of biomolecules, which hamper long-term operation, as required for continuous energy/fuel generation. Solving these issues is a key aspect of the coming transition toward sustainable energy production, which will undoubtedly be the focus of much scientific inquiry in the foreseeable future.

Despite being conceptually very simple, the hydrogen evolution reaction (HER) is a deceivingly difficult process. The reaction is only kinetically viable when its overpotential is lowered by a suitable electrocatalyst (i.e., a hydrogen-evolving catalyst or HEC), serving to decrease the activation barriers involved in the various chemical steps. The best solid-state HEC is platinum [90], which displays several favorable properties, such as an almost-zero hydrogen adsorption free energy, durability, and high electrical conductivity. However, it is a rare and expensive material, and this makes its wide adoption unfeasible. Therefore, large-scale production of hydrogen, as required by a full shift toward renewable energy sources, is associated with the discovery of artificial electrocatalysts based on abundant transition metals such as Fe, Co, Ni, Cu, Mo, and W [91,92]. Numerous inorganic-based materials have been shown to be active as HECs, particularly sulfides, selenides, carbides, nitrides, and phosphides. In addition to heterogeneous catalysts, a variety of small-molecule metal complexes have been developed for homogeneous catalysis, aimed at reproducing the exceeding catalytic performances of hydrogenases into simple, robust, and easily tunable systems. This section presents an overview of molecular catalysts conceived as electrocatalysts for HER. Starting from the smallest molecular scaffolds, organic and biomolecular architectures with increasing levels of complexity have been constructed around the catalytic metal site (Figure 3), introducing secondary sphere and long-range interactions that have a substantial impact on redox and catalytic properties. A summary of the catalytic parameters for electrochemical H2 production by selected catalysts is reported in Table 1.

Inspired by hydrogenases, a plethora of molecular complexes have been reported as structural and functional mimics of the active sites of these enzymes. Studies on these mimics have afforded further insight into the behavior and properties of their natural counterparts, as thoroughly described in the literature [21,95,96,97,98,99]. Indeed, deep investigation of some hydrogenase mimics enabled the elucidation of a number of details regarding the catalytic mechanism through isotope experiments, as described in a recent study [100]. The majority of model complexes are inspired by the diiron dithiolate carbonyl cluster of [FeFe]-hydrogenases, and well-established synthetic strategies are currently available [101,102,103]. The first generation of such structural mimics dates back to 1999, when limited structural information on hydrogenases active sites was available. Complexes bearing azadithiolate (ADT) and oxadithiolate ligands were reported in [104,105], until the nature of the bridgehead atom in the propanedithiolate (pdt) moiety was identified as nitrogen. Since then, substantial work has been focused on preparing model compounds resembling natural [FeFe]-hydrogenases in terms of primary coordination sphere and oxidation state of the diiron center [21]. However, most of these synthetic complexes usually catalyzed proton reduction at high overpotential. Strategies to overcome this major drawback involved either substitution of the carbonyl ligands or modification of the ADT bridge [106,107,108,109]. Abiotic phosphine or carbene ligands [110,111] are found in compounds with redox properties that are the closest to the thermodynamic potential for proton reduction. These ligands mimic the electronic properties of cyanide and increase the electron density of the diiron center, thus favoring protonation [112]. Another effective method for modulating the redox properties of diiron model systems is the inclusion of a properly substituted aromatic dithiolate bridgehead [113,114]. In particular, the use of a benezenedithiolate (bdt) ligand bearing electron-withdrawing substituents allowed lowering the overpotential for catalysis but compromised the catalytic performances in terms of turnover frequency. In addition to hydrogenase active site mimics, a variety of molecular HECs have been developed through a carefully designed electronic environment provided by the interaction between metal and ligands [115] (Figure 4). Examples of this are complexes with non-innocent S-donating ligands [116,117], such as dithiolene, thiosemicarbazone, and ene-1,2-dithiolate, which display some unique properties due to their unusual electronic structure (Figure 4a). Both experimental and computational studies revealed that the presence of redox-active ligands in these complexes is strongly linked to the catalytic activity [118]. Recently, Niu et al. [119] synthesized new dithiolene complexes of nickel able to evolve hydrogen in both aqueous and nonaqueous media. The compound with the highest performance ([BzPy]2[Ni(tdas)2]; BzPy = benzyl pyridinium, tdas = 1,2,5-thiadiazole-3,4-dithiolate) reached a turnover frequency (TOF) of 5.36 min−1 at a potential equal to −0.99 V (vs. SHE) in acetonitrile solutions (Table 1, entry 1). The same compound showed a higher TOF of 9.25 min−1 at −0.49 V (vs. SHE) in neutral buffer solutions (Table 1, entry 2). The electronic structure of the metal complex can also be rationally modified by imposing structural constraints on ligand geometry, as proposed by Drosou et al. [120]. Alkylation of two dithiolene sulfur atoms with a propyl moiety afforded compounds ([Ni(pbdt)]; pbdt = 2,2′-(propane-1,3-diylbis(sulfanediyl))dibenzene thiol) with increased metal basicity, able to reduce protons from TFA in DMF (Table 1, entry 3), albeit at comparable or higher overpotentials than their non-alkylated counterpart. Novel compounds based on bis(1,2,5-thiadiazole-3,4-dithiolate) ligands were prepared by Yin et al. [121] which were active as HECs in both DMF (Table 1, entry 4) and neutral buffer (Table 1, entry 5). Kamatsos et al. [122] investigated three heteroleptic compounds bearing different diamine ligands and 2-aminobenzenethiolate (2-amnt, Figure 4b and Table 1, entry 6), active as electrocatalysts for proton reduction in DMF with TFA as a proton donor. Albeit reaching remarkably high TOF (6120 min−1), the highest catalytic current for these complexes was observed at very high overpotentials. Replacement of nitrogen with an oxygen atom in the 2-amnt ligand yielded analogous heteroleptic diamine oxathiolate nickel complexes (Figure 4c and Table 1, entry 7). This enabled catalysis at lower overpotentials, but significantly decreased both TON and TOF [123]. In another recent study, Chen et al. [124] synthetized five novel nickel complexes bearing o-methyldithiophosphate and aminodiphosphine monosulfide ligands (Figure 4d), which showed HEC activity in acetonitrile using TFA as a proton source. The electrocatalytic studies on these complexes showed that the catalytic performances depend on the nature of the substituent in the aminodiphosphine monosulfide ligand. In particular, the presence of a 4-CH3C6H4 group leads to the highest TOF and the lowest overpotential (Table 1, entry 8) due to the electron-donating effect. Complexes including nitrogen and oxygen-based ligands have also been extensively studied as HECs, and several classes have emerged as promising, including diamine dioxime, polypyridine, corrole, porphyrin, tetraimine, diglyoxime, and cyclam complexes, in most cases containing cobalt [125]. Cobalt diimine dioxime compounds (Figure 4e) have long being known as HECs [126] and several water-soluble analogues have also been reported. In particular, Peters and coworkers evaluated the effect of the bridging groups within the macrocycle and found that the proton bridged compound ([Co(DO)(DOH)pn(OH2)2]2+, Table 1, entry 9) displayed the highest TON (23) among the series, catalyzing proton reduction at notably low overpotential (260 mV) [127]. A recent study by Sun et al. [128] granted further insight into the mechanism. The authors synthesized and studied two cobalt diimine dioxime complexes and carried out a detailed quantitative analysis of their catalytic activity, which confirmed a previously proposed catalytic mechanism. A nickel salicylketoxime complex (Figure 4f) was instead reported by Wang et al. [129] to be active as an electrocatalyst in neutral aqueous buffer, achieving a TON of 1428 over 1 h of electrolysis (Table 1, entry 10). Additionally, the first Cu(I) quinoxaline complex (Figure 4g) active as an HEC was reported by Drosou et al. [155], and its activity was investigated by electrochemical and DFT methods, which also provided insight into the reaction mechanism. Polypyridine complexes (Figure 4h) have also been abundantly explored as potential HECs, as summarized in some excellent reviews [156,157]. Novel 5- and 6-coordinated cobalt complexes with polypyridyl ligands were designed by Liu et al. [130], and their properties were studied both experimentally and theoretically. Both complexes were found to be active HECs when using triethylammonium tetrafluoroborate (Et3NH+BF4−) as the proton source in DMF, with the most active complex ([Co(ppq)(MeCN)(TfO)]+; Table 1, entry 11) reaching a remarkably high TOF of 336,000 min−1. Voltammetry experiments also allowed demonstrating the dissociation of a phosphine axial ligand from one of the two complexes, prior to hydrogen evolution. An unusual cobalt complex was proposed by Kumar Padhi et al. [131], in which terpyridine molecules were conjugated to ferrocene ([Co(Fc-tpy)2]2+; Table 1, entry 12). The activity of the complex was tested in DMF/H2O solutions with acetic acid as the proton source and compared to the parent terpyridine compound lacking the ferrocene moieties. The complex conjugated to ferrocene showed a 200 mV lowered overpotential and a doubled TOF compared to the simple terpyridine cobalt derivative. These results were attributed to the electron-donating character of ferrocene. Another interesting example of sophisticated ligands was provided by Cui et al. [132], who prepared two bistriazolylpyridine cobalt complexes (Figure 4i and Table 1, entries 13–14) through click chemistry and performed a number of spectroelectrochemical experiments on them, also involving the use of cobaltocene as reductant. Both resulted active HECs in acetonitrile with acetic acid, and their redox and catalytic behaviors were shown to be influenced by substitution on the triazole moiety. The synthesis of copper(II) phenolate selenoether complexes (Figure 4j) was recently accomplished by Upadhyay et al. [133], who showed that these compounds were air-stable and catalytically active as electrocatalysts. The presence of electron-withdrawing groups on the ligands (e.g., formyl group) was also demonstrated to yield the highest TOF values (Table 1, entry 15). Lastly, some binuclear cobalt complexes (Figure 4k and Table 1, entry 16), recently reported by To et al. [134], were also found to be highly active HECs in DMF with acetic acid as the proton source. As mentioned, numerous porphyrin-based complexes have also been employed as homogeneous HECs (Figure 4l) [158,159,160,161]. The different mechanisms possible for hydrogen evolution by porphyrins were compared and discussed in a review on the subject [162]. Chou et al. [135] recently prepared TFP (5,10,15,20-tetrakis(4-fluoro-2,6-dimethylphenyl)porphyrin) complexes of Cu(II), Zn(II), and Ni(II). Among them, only the Cu(II) derivative exhibited proton reduction activity in acetonitrile with benzoic acid as the proton source (Table 1, entry 17). Novel corrole complexes of Co, Cu, Mo, Mn, and Ni have also been described in very recent work (Figure 4m) [136,137,138,163]. Gross and coworkers reported a nickel corrole derivative displaying HER activity in acetonitrile with low overpotential (360 mV, Table 1, entry 18) [136]. A manganese corrole catalyst was shown to catalyze proton reduction in both DMF and mixed aqueous/acetonitrile solution (Table 1, entries 19–20), reaching the highest TOF in the latter condition, but suffering from a high overpotential (1.027 V) [137]. To the best of our knowledge, the most active corrole-based catalysts are the copper complexes developed by Sudhakar and Panda [138], who evaluated the effect of chlorine substitution on electrocatalytic hydrogen evolution. The octafluoro-corrole derivative (CuCl8, Table 1, entry 21) showed an impressively high TOF (1.47 × 107 min−1) for proton reduction in acetonitrile with low overpotential (360 mV). Overall, an impressive amount of research effort has been directed toward the development of small-molecule complexes mimicking hydrogenase function. Despite the achievement of notable results, there is still ample margin for improvement, particularly regarding the overpotentials, TOF, costs, and toxicity of the proposed compounds. The catalytic activity and overpotentials of molecular HECs depend on a number of factors, including the identity of the metal and its ligands, electronic structure of the complex, and operational conditions. Minimization of the overpotential down to a nearly vanishing value is a crucial aspect to obtain synthetic analogues of natural hydrogenases. However, optimizing the electronic structure and substitution of metal complexes to achieve minimum overpotential is a challenging task, since hard ligands favor protonation of a low-valent metal atom, but hinder proton reduction [164]. Moreover, non-innocent ligands might lower the overpotential by allowing storage of reducing equivalents, as proposed for dithiolene complexes [116]. The presence of internal bases has also been shown to decrease the overpotential by facilitating protonation/deprotonation steps in the catalytic cycle [21]. Due to this intricate interplay of factors, it is difficult to delineate a unique structure–function relationship for all the hydrogen-evolving catalysts, although correlations between substitution and catalytic activity have been reported for selected systems [165].

Small-molecule complexes offer the advantage of being easily produced through synthesis. Although active as HECs, the lack of the protein environment hampers them in reproducing the efficiency of the natural counterparts, such as the capability of performing well in water. Indeed, the unique ability of hydrogenases to reversibly catalyze proton reduction is certainly connected to their highly evolved protein structures [166]. Despite the apparent simplicity of HER, this reaction involves electron transfer (ET) and proton transfer (PT) steps, which are usually interconnected as proton-coupled electron transfer (PCET) processes [167,168]. Incorporating acid–base groups into molecular electrocatalysts facilitates this mechanism by positioning the proton donor and acceptor close to each other, thus resembling natural catalysts. Following nature’s footpath, the insertion of proton relays near the active site of synthetic complexes has been actively pursued. Several mononuclear nickel bis diphosphine catalysts, incorporating a pendant amine group in the ligand backbone ([Ni(PR’2NR2)2]2+ complexes), have been developed to investigate the role of the base in catalysis [169,170]. Interestingly, the catalytic activity for complexes containing properly positioned pendant amines occurs at much higher rates and/or lower overpotentials, compared to analogous catalysts lacking the base group (Table 1, entry 22) [139]. Such improvement has been attributed not only to a facilitated proton transfer from the solution to the catalytic site, but also to a decrease in the activation barrier for H–H bond formation (or cleavage) [139,171,172,173]. Substitution of Ni(II) with Co(II) afforded complexes able to catalyze electrochemical H2 production only when a pendant amine was incorporated in the diphosphine ligand, further confirming the essential role of the proton relay in the catalytic reaction (Table 1, entry 23) [140]. This approach was fruitfully applied also to structural mimics of the [FeFe]-hydrogenase active site. Dey and coworkers recently developed a series of azadithiolate (ADT) diiron hexacarbonyl complexes carrying an ortho-substituted aryl moiety bound to the bridging nitrogen of ADT [141]. Ortho-substitution of the aryl group favors a tetrahedral geometry at the nitrogen atom, which disfavors electron delocalization to the aromatic ring, in contrast with the nearly planar geometry observed in other N-aryl ADT complexes. This leads to an increased electron density at the well-oriented, pyramidal nitrogen, enabling its protonation at pH 5.5, whereas unsubstituted N-aryl derivatives show proton reduction activity only in strongly acidic conditions (Table 1, entry 24). Thus, despite not affecting the formal redox potential of the diiron center, the presence of an active proton shuttle allows HER to occur at very low overpotential (180 mV). In follow-up work, such systems were shown to promote both proton reduction and H2 oxidation as natural enzymes, and they were implemented in the construction of a fuel cell [93]. Along these lines, Rauchfuss and coworkers developed, for the first time, a synthetic model containing the three functional components of the [FeFe]-hydrogenase active site: the reactive diiron center, an amine as the proton shuttle, and a ferrocene derivative as a one-electron redox partner [174]. Other relevant examples of this strategy are the “hangman porphyrins”. These complexes typically contain a xanthene ring that places a pendant carboxylate or amine group over the metal ion [175,176,177]. The boosting effect of introducing an acid–base functionality was first demonstrated by Nocera and coworkers, who reported the lower overpotentials and improved catalytic activities of hangman porphyrins compared to the non-hangman analogues. These systems also provided insights into the mechanism of HER catalyzed by cobalt, nickel, or iron porphyrins, mainly based on electrochemical and spectroscopic data coupled to DFT analysis [162,178,179]. Starting from this pioneering work by Nocera and coworkers, the hangman strategy has been subsequently expanded, by appending increasingly complex molecular frameworks to porphyrin-like catalysts [180,181,182]. As a notable example, Cao and coworkers designed a hydrogen-bonded water network by positioning a crown ether over a cobalt corrole catalyst, which significantly enhanced the catalytic performances by facilitated proton delivery (Table 1, entry 25) [142]. However, hangman porphyrins, as the majority of molecular complexes discussed above, have been typically studied in organic solvents due to their scarce solubility in water. This represents an important drawback for their potential application in fuel cells. Several alternative strategies have been proposed for the construction of secondary sphere interactions around the metal-containing catalysts, concurrently enabling good performance in water. For example, these strategies foresee the inclusion of molecular hydrogenase catalysts into cyclodextrins [183,184], micelles [143], and metal–organic supramolecular cages [144,185]. Encapsulation of a sulfonated derivative of ADT diiron hexacarbonyl into β-cyclodextrin provided water solubility to the insoluble complex, but disfavored electrochemical proton reduction, as demonstrated by increased overpotentials. This phenomenon was ascribed to the steric constraints imposed by the cyclodextrin cavity, hampering structural rearrangements during the catalytic cycle or inhibiting proton delivery to the active site [183,184]. Conversely, dispersion of benzenedithiolate diiron hexacarbonyl into micelles formed in aqueous sodium dodecyl sulfate (SDS) provided efficient electrochemical H2 evolution (Table 1, entry 26), occurring at an overpotential of <500 mV (pH 3) [143], which is competitive with the values reported for most molecular catalysts. The encapsulation of a pyridine-substituted diiron hydrogenase model into a synthetic Zn-porphyrin cage (Fe4(Zn-L)6) significantly decreased the catalytic overpotential (by 150 mV) at the expense of a slightly lower catalytic rate with respect to the same complex in freely diffusing conditions (Table 1, entry 27) [144].

As described above, with the aim of tuning and controlling the activity of artificial H2 evolving electrocatalysts, more and more sophisticated and diverse environments have been continuously exploited. As a step further in the development of artificial hydrogenases closely resembling the natural counterparts, the field has broadened to include more elaborate systems based on the construction of a peptide- or protein-based architecture to host the catalytic center. Small peptides binding hydrogenase model complexes lay in the middle ground between natural metalloenzymes and small-molecule complexes, and they provide the chance to carefully modulate secondary sphere interactions while also allowing activity in water. Early work in this field was focused on anchoring dithiolate bridged diiron hexacarbonyl clusters to peptide scaffolds, mimicking the active site of [FeFe]-hydrogenases. Dutton and coworkers reported the first peptide-based structural model of hydrogenases consisting of a de novo designed helical peptide capable of binding the diiron cluster through cysteine residues [186]. Along these lines, other groups have developed [FeFe]-hydrogenase models endowed with redox activity, exploiting either natural or artificial dithiol containing peptides to bind the diiron core [187,188,189]. Notably, Jones and coworkers reported an alternative strategy to link a dithiolate bridged diiron core to a peptide scaffold, which does not involve thiolates as anchoring sites. In this approach, an unnatural phosphine amino acid was incorporated into the peptide sequence, for binding the diiron cluster and replacing a terminal carbonyl [190]. The resulting phosphine-bound diiron pentacarbonyl complex displayed improved electrocatalytic properties, which were ascribed to the higher electron donating properties of the phosphine compared to CO, acting as a mimic of CN− found in hydrogenases [191]. Shaw and coworkers introduced an amino-acid-based outer coordination sphere by further modification of the diphosphine ligand with a mono- or dipeptide in the previously described [Ni(PR’2NR2)2]2+ complexes [145,192,193]. In particular, simple incorporation of a diglycine ethyl esther (DPE = dipeptide esther) into the ligand framework imparted water solubility to this class of molecular catalysts while boosting their catalytic activity both in proton reduction and in H2 oxidation (Table 1, entry 28) [145]. The carboxylic groups of glycines may facilitate proton delivery from the solvent to the metal and vice versa, creating a minimal proton channel [192]. Despite countless mimetics of the [FeFe]-hydrogenase diiron cluster being reported, functional models of [NiFe]-hydrogenases are rare [194,195]. The absent or little activity displayed by structural mimics of these enzymes underlines that the protein environment plays a fundamental role in dictating the properties of the metal cluster. Dutta et al. reported a synthetic methodology to construct minimal peptide-based model of [NiFe]-hydrogenases using a fragment derived from the enzyme nickel superoxide dismutase (NiSOD) [196]. Although the synthetic approach is promising, functional peptide-based models of this class of enzymes are still missing. The effect that a protein matrix can exert on catalyst robustness and overpotential was demonstrated by Bren and coworkers, who screened cobalt porphyrins covalently bound to peptide scaffolds with different levels of complexity (Figure 5). Among them, the first and simplest was Co-MP11, which consisted of the cobalt-substituted cytochrome c-derived heme-peptide complex, known as microperoxidase-11 (Figure 5a and Table 1, entry 29). Co-MP11 performed ~25,000 turnovers in electrocatalytic hydrogen evolution at an overpotential of ~850 mV, but showed a significant drop in the activity after only 15 min of electrolysis [146]. Conversely, the cobalt-substituted mutant of Hydrogenobacter thermophilus cytochrome c552 (Ht-CoM61A, Figure 5b, Table 1, entry 30) was more robust, lasting for 6 h of electrolysis and displaying a 10-fold enhanced turnover number (~270,000 TON). However, almost no improvement in the catalytic overpotential of Ht-CoM61A (830 mV) was observed with respect to Co-MP11 [147]. This drawback was recently addressed thanks to the joint efforts of the Bren and Lombardi labs, with the synthetic cobalt deuteroporphyrin-containing miniprotein Co-MC6*a (Figure 5c and Table 1, entry 31) [148,197]. Co-MC6*a promotes proton reduction from nearly neutral aqueous solutions (pH 6.5) with the highest catalytic current at an overpotential of 680 mV, which is a significantly lower value compared to both Co-MP11 and Ht-CoM61A. Furthermore, investigations of Co-MC6*a revealed that enhancing peptide folding via the addition of TFE (2,2,2-trfifluoroethanol) as a cosolvent lowered overpotential to 520 mV, while increasing catalyst longevity up to 1.3 h with respect to Co-MP11 [148]. Remarkably, this miniature enzyme performed more than 230,000 turnovers, a value comparable to Ht-CoM61A. Artificial hydrogenases have been also produced by incorporating small-molecule catalysts into protein scaffolds, with the aim of imparting enzyme-like features to synthetic complexes. This approach, usually referred to as the host–guest strategy [198,199,200], takes advantage of stable and evolutionary selected protein architectures to create a tunable protein environment around a catalytic metal site. Cobalt bis-glyoxime catalysts have been incorporated into several natural proteins replacing the native metal cofactor. Artero and coworkers reported two biohybrid hydrogenase catalysts via reconstitution of apo sperm-whale myoglobin (SwMb) with two different cobaloxime moieties (Table 1, entry 32) [149]. The electrocatalytic activity of the cobalt centers was retained in both systems, which performed proton reduction in aqueous conditions with low overpotential (~200 mV). However, the operational stability of these biohybrids remained limited to few turnovers (TON = 3.2 at pH 7). Slight improvements in terms of TON (up to 12) were observed upon insertion of the same cobaloxime catalysts into the heme-binding pocket of heme oxygenase (HO) as the host protein (Table 1, entry 33) [150]. In both cases, the complexes alone showed better performance compared to when embedded into the proteins. This finding was ascribed to a bimolecular mechanism involving two molecules, which is inhibited when the catalysts are housed within the protein scaffold. Ghirlanda and Fillol developed cobalt-based biohybrid systems exploiting the biotin–streptavidin technology, which relies on the exceedingly high binding affinity of the protein streptavidin (SA) toward biotin [151]. Biotinylated aminopyridine cobalt complexes encapsulated into SA retained their electrocatalytic activity, displaying a decreased overpotential compared to the lone complexes (Table 1, entry 34). The simple substitution of iron with nickel into Desulfovibrio desulfuricans rubredoxin led to an efficient hydrogenase model (NiRd), showing an electrocatalytic behavior closely resembling that of Desulfovibrio vulgaris [NiFe] hydrogenase, with H2 production rates that are comparable to those of the native enzyme (Table 1, entry 35) [152]. A library of NiRd mutants was generated using site-directed mutagenesis to evaluate the effect of mutations in the secondary coordination sphere. Interestingly, the introduction of a carboxylate residue, mimicking that found in the [NiFe] hydrogenase, considerably increased the turnover frequency, likely resulting from a proton transfer pathway into the active site, without significantly affecting the overpotential (Table 1, entry 36) [153]. Recently, a rational metalloprotein design approach was applied by Chakraborty and coworkers to repurpose a copper storage protein (Csp1) into a nickel-binding protein (NBP) able to promote hydrogen evolution [154]. A tetrathiolate nickel-binding site was engineered into Csp1, inspired by the active site of [NiFe] hydrogenases. Ni-NBP performed 210 turnovers over 1 h of bulk electrolysis, displaying an overpotential which is still far from those of native hydrogenases, but comparable to those of most molecular catalysts (Table 1, entry 37).

Photocatalytic hydrogen evolution is the most attractive process, as, in principle, it could exploit sunlight as the energy source. At the same time, it is also the most challenging, since it requires the presence of three components working in synergy: a photosensitizer, a catalyst, and an electron donor [201]. A suitable light-harvesting unit should be characterized by a strong absorption and a long-lived excited state to facilitate charge separation. Furthermore, the redox potentials of the photosensitizer and catalyst should be well matched, allowing the desired charge transfer among them to occur. Lastly, as H2 evolution is a reductive process, an electron reservoir is required to bring the photosensitizer back to the initial state. Sacrificial electron donors are typically used to this purpose, but coupling proton reduction with an oxidative process (such as water oxidation in photosynthesis) is highly desirable. Several reagents have been used as sacrificial electron donors in photocatalytic hydrogen evolution, with the most used being ascorbic acid (AscOH), alcohols, and aliphatic amines [202,203]. The choice of the electron donor is as important as that of the photosensitizer, since it must be oxidized to an inert species which does not interfere with catalysis. However, the use of a sacrificial reagent represents a limitation to the scale-up of the process, because it requires large quantities or repeated addition of the reagent, which may not be affordable from a sustainability and/or an economic point of view. One possibility to address this problem consists of using electron donors which can be reversibly regenerated upon light irradiation, as reported by Girault and coworkers [204,205]. In a homogeneous catalytic system, the energy transfer is limited by diffusion, which brings the active components close enough for the process to occur. Possible strategies to optimize the efficiency of the photoelectron transfer consist of creating covalent links between the catalyst and the photosensitizer or of gathering the active components into a micro- or nanometer-sized reactor. This problem is not encountered with heterogenous catalysts, as light harvesting and catalysis occur in the frame of a photoactive material. Since 1970s, there has been growing interest in the fabrication of highly efficient semiconductor-based photocatalysts, which include TiO2 [206], ZnO [207], Fe2O3 [208], and CdS [209]. Heterogeneous systems offer the advantage of being highly stable and long-lived; however, unlike molecular catalysts, they have limited light-harvesting abilities and tunability. Further, these materials suffer from other common drawbacks, including wide bandgaps, conduction band potentials lower than the proton reduction potential, toxicity, and high cost. The need to overcome the inherent limitations of inorganic materials has prompted great interest towards the study of organic polymers as photocatalysts [210]. The first and possibly the most studied photocatalytic polymeric semiconductor is graphitic carbon nitride (g-C3N4), due to its facile synthesis, appealing electronic band structure, and high physicochemical stability [211]. In the past few years, a new class of highly porous organic polymers emerged as prominent candidates for photocatalytic hydrogen evolution, namely, covalent organic frameworks (COFs) [212,213,214]. COFs are made up by molecular building blocks; hence, they grant the possibility of modulating the different functions, fundamental to the photocatalytic process: light harvesting, charge separation, charge transport, and catalysis. Furthermore, COFs are characterized by nanometer-sized structural pores entailing a high interaction surface area and providing accessibility to sensitizers, sacrificial components, and cocatalysts throughout the material [215]. A comprehensive analysis of heterogeneous systems for hydrogen production is beyond the scope of this review; interested readers can refer to the literature [216,217]. In the remainder of this section, homogeneous photocatalytic systems involving molecular and biomolecular components are described. A summary of the catalytic parameters from selected photocatalytic H2 production systems, described below, is reported in Table 2.

The choice of the most appropriate photosensitizer is of paramount importance in light-driven catalysis. Ruthenium, iridium, rhenium, and platinum polypyridine complexes are usually preferred (Figure 6a–d), since they display intense metal-to-ligand charge transfer absorption bands in the visible region, populating a long-lived excited-state, which provides enough time for the electron transfer [246,247]. However, their high cost, toxicity, and difficult recyclability hamper their use on a large scale, thus prompting research toward noble-metal free photosensitizers. In contrast, organic dyes such as fluorescein and eosin Y (Figure 6e,f) show higher absorption intensity of light and fluorescence quantum efficiency compared to metal complexes, but have generally lower chemical stability and shorter excited-state lifetime [248,249]. In order to overcome these issues, the use of semiconductor quantum dots (QDs) as photosensitizers has gained much interest (Figure 6g) [250,251]. QDs are nanocrystals with a diameter of few nanometers, characterized by a unique electronic structure which is intermediate between that of bulk semiconductors and discrete atoms. Due to their small dimensions, the charge carriers (electron and holes) experience a quantum confinement effect; thus, they can occupy only discrete energy levels. As a result, QDs display strong absorption and fluorescence at specific wavelengths. [252,253]. These wavelengths can be precisely tuned by altering the QD composition and size. In particular, the energy gap between the conduction and valence bands increases as the size of QDs decreases, providing extended stabilization to the charge-separated state [254]. The highest photoluminescence quantum yields have been reported for cadmium containing nanocrystals, such as CdS, CdSe, and CdTe; thus, research on environmentally friendly and highly efficient QDs is growing [255,256]. Nature developed an astonishing light-harvesting unit, namely, photosystem I (PSI, Figure 6h). It is the light-harvesting machinery found in bacterial reaction centers and works in synergy with photosystem II (PSII) to achieve photosynthesis [201,257]. PSI is provided with a molecular antenna made of multiple carotenoid and chlorophyll molecules on the protein periphery. Light excitation initiates a cascade of rapid electron-transfer steps between the cofactors, starting with the primary electron donor (P700) and terminating with a [4Fe-4S] cluster (FB). The quantum yield of PSI approaches 1; thus, almost every absorbed photon is converted to the charge-separated state (P700+ FB−), which is characterized by a long lifetime (~60 ms) [258]. The impressive photosensitizing properties of PSI have attracted much interest in switching its photocatalytic activity from NADH synthesis to hydrogen evolution, taking advantage of its ability to provide electrons to partners other than natural ones [259]. Biohybrid systems have been constructed by combining molecular or biomolecular catalysts with PSI.

Over the last few decades, great efforts have been focused on the design of photocatalytic systems coupling a [FeFe]-hydrogenase model complex with a variety of light-harvesting units [218,260,261,262,263]. Unlike native [FeFe]-hydrogenases, their mimics undergo gradual degradation and loss of CO under light irradiation over a few hours, resulting in relatively low TON [22]. Substitution of one carbonyl with a phosphine ligand results in improved photostability, but also increases the electron density at the diiron site [264,265]. This effect, which is beneficial for electrocatalysis, is detrimental for the photoinduced process, because it disfavors the charge transfer from the photosensitizer. The use of weakly donating phosphines such as P(pyr)3 (tris(N-pyrrolyl)phosphine) has been exploited for the construction of functional photocatalysts, performing up to 466 turnovers over 8 h of irradiation upon matching with a iridium(III) photosensitizer (Table 2, entry 1) [218]. Other strategies aimed at increasing the stability of diiron carbonyl complexes consist of replacing the aliphatic dithiolate bridged ligand with an aromatic one, such as benezenedithiolate (bdt) [266]. Due to their structural rigidity, aromatic bridgeheads enhance catalyst robustness and offer the possibility of easily tuning redox properties by introducing substituents on the conjugated system [125]. The incorporation of electron-withdrawing substituents such as chlorine atoms on the aromatic bridging ligands was proven to facilitate the electron transfer from the photosensitizer, but compromised complex stability by weakening the Fe–S bonds (Table 2, entry 2) [114,219]. This problem was addressed by Hou and coworkers, who exploited a dinaphthalene-derived dithiolate ligand (μ-BNT) for the construction of an efficient photo-activated catalyst [220]. The photocatalytic system composed of the dinaphthalene diiron complex, eosin Y as a light-harvester, and triethylamine (TEA) as a sacrificial reagent was capable of performing 404 turnovers under 4 h of visible-light irradiation in a mixed water/acetonitrile solution (Table 2, entry 3), reporting the highest TON among the aromatic bridged diiron carbonyl complexes. Another subset of [FeFe]-hydrogenase mimics was developed by incorporating hydrophilic groups in the dithiolate bridge, increasing water solubility of the complexes. In this context, Wu and coworkers developed two diiron model complex bearing sulfonate-substituted thiol ligands, using either a bridged dithiolate (Fe2S2-SO3Na) or two separate thiolate ligands (Fe2S2–2SO3Na) [221]. The two analogues displayed similar catalytic behavior, which was found to be strongly dependent on the nature of the photosensitizer. In detail, TON values of 114 and 178 were reported in the presence of [Ru(bpy)3]2+ (Table 2, entries 4–5), while values of 18,800 and 26,500 were found in the presence of CdSe QDs (Table 2, entries 6–7) for the bridged and unbridged complex, respectively. This comparative study confirmed the importance of an efficient electron transfer, from the photosensitizer to the catalytic unit, in photocatalytic hydrogen-evolving systems, highlighting the superior performances of QDs with respect to ruthenium(II) polypyridyl complexes [221]. One possible strategy to enhance the stability and performances of a homogeneous photocatalytic system consists of constructing supramolecular assemblies, where the distance between the catalyst and the light-harvesting unit can be precisely controlled. Covalently linked dyads were constructed associating diiron carbonyl complexes with several kinds of photosensitizers, where the organic dithiolate bridge serves as a connection between the components [267,268,269,270,271] (Figure 7). Rigid and conjugated organic molecular structures are preferred as spacers between the catalytic and the photosensitizer units, in order to favor the charge transfer process. However, moderate activities have been reported for covalent dyads containing zinc porphyrins [268,269] (Figure 7a), rhenium complexes, or ruthenium complexes as photosensitizers. A new class of organic, silicon-containing photosensitizer was recently reported, displaying the advantage of fast electron transfer through the Si–CH2–S linker [222,272] (Figure 7b). The system composed of 1-silafluorene linked to diiron dithiolate hexacarbonyl complex performed 539 turnovers over 7 h of light irradiation under optimized reaction conditions (Table 2, entry 8) [222]. In addition to model complexes of hydrogenases, a variety of synthetic molecular catalysts have been studied for photocatalytic hydrogen production. Among them, Ni-based catalysts have demonstrated excellent catalytic activity and received much attention among researchers. A series of mononuclear nickel(II) complexes with diphosphine ligands, developed for electrochemical proton reduction, were also tested for light-driven catalysis. Remarkably, the [Ni(PR’NRPR’)2]2+ complex, reported by DuBois and coworkers, was tested by Holland and Eisenberg et al. in the presence of different photosensitizers, displaying long-term stability and performing 2700 turnovers when paired to [Ru(bpy)3]2+ (Table 2, entry 9) [223]. In addition to the phosphine ligand, complexes with dithiolate or aminothiolate ligands were also developed, taking inspiration from the ligand composition of [NiFe]-hydrogenase active sites. Pyridinethiolate complexes such as [Ni(X-pyS)3]− (pyS = pyridine-2-thiolate; X = 5-H, 5-Cl, 5-CF3, 6-CH3) and [Ni(4,4′-Y-2,2′-bpy)(pyS)2] (Y = H, CH3, OCH3) displayed impressive photocatalytic properties and performed up to 7300 turnovers in the presence of fluorescein as photosensitizer (Table 2, entry 10) [224]. Interestingly, in the absence of substituents on the pyS ligand, the nickel complex assumed a square planar geometry, while distorted octahedral geometries were found in the other cases. Mechanistic studies showed that a key step in the catalytic cycle of these catalysts is the reversible dechelation and protonation of a pyridine nitrogen, which enables the formation and coupling of a metal hydride with the N-bound proton to produce H2 [224]. On the basis of these findings, nickel(II) bischelate complexes with a square-planar coordination have more recently gained attention [123]. The [Ni(bpy)(mp)] (mp = 2-hydroxythiophenol) complex displayed notable durability, being capable of promoting H2 evolution under light irradiation for more than 100 h, reaching exceedingly high turnover numbers (Table 2, entry 11). Cobalt-based catalysts have been also privileged candidates for hydrogen production, due to their relatively low cost and high tolerance to O2 [125,149,273]. A wide variety of cobalt complexes have been tested for light-driven hydrogen evolution, mostly involving macrocyclic or pseudomacrocyclic square-planar ligand arrays. Among them, difluoroboryl-bridged cobaloximes such as [Co(dmgBF2)2] (dmg2− = dimethylglyoximato) have been the most studied, owing to their stability against hydrolysis in acidic conditions [274]. Encouraging results were obtained by Artero and coworkers, which explored the use of cyclometallated iridium diamine ([Ir(ppy)2(phen)]+; ppy = 2-phenylpyridine, phen = 1,10-phenantroline) or tricarbonylrhenium diamine ([ReBr(CO)3(phen)]) as photosensitizers achieving 165 and 273 TON, respectively (Table 2, entries 12–13) [225]. H2 production was improved by enhancing the long-term stability of the rhenium photosensitizer upon substitution of the labile bromide ligand with thiocyanate [226]. With the highly stable photosensitizer [Re(NCS)(bpy)(CO)3], the stability of cobaloxime catalyst became performance-limiting, and up to 6000 turnovers with respect to the rhenium complex were performed (Table 2, entry 14). Eisenberg et al. also focused on photocatalytic activity of cobaloxime complexes, with a platinum terpyridyl acetylide ([Pt(ttpy)(C≡CPh)]+, ttpy = 4′-p-tolylterpyridine) [227] or eosin Y as photosensitizers [275], achieving up to 2400 TONs with the Pt complex (Table 2, entry 15). Cobalt porphyrins have been predominantly studied for electrochemical proton reduction, and only limited reports are available for photochemical H2 evolution with this class of catalysts (Table 2, entries 16–17) [228,229]. A notable example was reported by Hung et al., who developed a sulfonate-substituted tetraphenyl porphyrin (CoTPPS) capable of performing ~6400 turnovers in neutral water and in the presence of [Ru(bpy)3]2+ as photosensitizer [229].

Combining hydrogenase model complexes with peptides or protein scaffolds has proven to be a valuable strategy to modulate the surrounding environment of the catalyst, providing additional stabilization and/or altering the charge transfer mechanism. Feng and coworkers constructed protein-based nanoreactors by incorporation of [FeFe]-hydrogenase mimics into protein-based supramolecular structures [230,231]. First, the authors achieved catalyst encapsulation into horse spleen apoferritin ([FeFe]-HSF) [230]. The negatively charged inner surface of apo-HSF, due to the presence of multiple glutamate residues, was exploited to cage up to 312 diiron carbonyl moieties. The photocatalytic hydrogen evolution activity was studied in mild acidic aqueous solutions, using [Ru(bpy)3]2+ as photosensitizer and ascorbate as sacrificial electron donor (Table 2, entry 18). Interestingly, even though penetration of the ruthenium complex into the protein interior was negligible, photocatalytic activity per single diiron catalyst was improved ~8-fold upon inclusion into apo-HSF. This finding suggests that the electron transfer process does not need the catalyst and the photosensitizer to be in close proximity, but it may occur through a multistep electron tunneling pathway across the protein matrix. Similar results were reported in a subsequent work, where the diiron dithiolate complex was entrapped into self-assembling protein nanogels formed by ovalbumin (Table 2, entry 19) [231]. Remarkably, the composite nanogel ([FeFe]-OVA) showed a 15-fold improvement in photocatalytic activity compared to the freely diffusing catalyst, suggesting an enhanced electron transfer efficiency mediated by the protein framework. Artificial hydrogenases obtained by protein design and engineering have also been studied for light-driven catalysis. Hayashi and coworkers anchored the diiron carbonyl cluster to cytochrome c, exploiting the cysteine residues in the CXXC heme-binding motif [232]. Substitution of the native cofactor with the synthetic diiron complex endowed the resulting conjugate ([FeFe]-Cyt c) with photocatalytic H2 production activity in the presence of [Ru(bpy)3]2+ as a photosensitizer and ascorbate as a sacrificial electron donor (Table 2, entry 20) [232]. In a follow-up study, aimed at improving the efficiency of electron transfer, the authors selected an 18-residue fragment (Pep18) of the cytochrome c sequence as a peptide matrix to covalently bind both the diiron carbonyl core and the photoactive ruthenium complex (Table 2, entry 21) [189]. In particular, the native heme coordinating histidine was chosen as the anchoring site for [Ru(bpy)(tpy)]2+ (bpy = 2,2′-bipyridine, tpy = 2,2′:6′,2″-terpyridine), and the two cysteines in the CXXCH sequence acted as ligands for the diiron complex. The intramolecular photochemical system ([FeFe][Ru]-Pep18) promoted H2 evolution in aqueous solution at pH 8.5 in the presence of excess ascorbate as sacrificial reagent, albeit displaying a limited turnover number (TON ~9). Addition of the imidazole-bound Ru complex ([Ru(bpy)(tpy)(im)]2+) as an external photosensitizer to the [FeFe]-Pep18 conjugate, did not lead to H2 evolution, evidencing that efficient electron transfer might occur only within the peptide framework. Further expanding the range of protein scaffolds, the Hayashi group covalently attached the diiron cluster into a rigid β-barrel protein, i.e., nitrobindin (NB), containing a suitable cavity to host the metal complex [94]. Inclusion of the [FeFe]-hydrogenase active site mimic within the barrel yielded an artificial metalloenzyme ([FeFe]-NB) capable of performing up to 130 turnovers over 6 h of activity in the presence of [Ru(bpy)3]2+ as photosensitizer (Table 2, entry 22). The catalytic complex without the protein matrix displayed similar total hydrogen production but required only 2 h to reach the plateau. The decreased initial reaction rate was attributed to a lowered accessibility of the ruthenium complex to the diiron site when embedded in the protein matrix. Biomimetic diiron hexacarbonyl complexes have also been encapsulated into SA using the biotin–SA technology (Figure 8a and Table 2, entry 23) [233]. Incorporation in the protein scaffold provided the diiron catalyst ([FeFe]biot-SA) with prolonged stability and enhanced activity compared to the isolated complex. The conjugation of synthetic complexes to proteins has not been limited to structural mimics of hydrogenases; it has also been applied with other classes of molecular catalysts. Ward and coworkers incorporated a biotinylated cobalt pentapyridyl complex ([CoBr(appy)-Biot]) into several SA mutants (Figure 8b,c and Table 2, entries 24–25) to evaluate the effect of the amino acids close to the catalytic center on H2 evolution [234]. The presence of lysine residues near the metal site resulted not only in increased turnover number with respect to the wildtype protein, but also in a higher reaction rate, suggesting a beneficial proton shuttling effect. Ghirlanda and coworkers inserted cobalt protoporphyrin IX (CoPPIX) into myoglobin, replacing the native heme cofactor. When coupled with Ru(bpy)32+ in the presence of ascorbate, the beneficial role of the protein scaffold was demonstrated by a threefold increased TON (518 at pH 7) of Co-Mb with respect to the free porphyrin (Table 2, entry 26) [235]. In a subsequent work, the same authors explored the effect of replacing heme with CoPPIX in Cytochrome b562 [236]. Cofactor substitution afforded a photoactive enzyme (Co-cyt b562), whose activity could be modulated through mutagenesis in the primary coordination sphere (Table 2, entries 27–29). In particular, substitution of the axial methionine with an alanine or an aspartate residue led to a ~2.5-fold improvement in total hydrogen production with respect to the wildtype protein. Such enhancement was attributed to a synergic effect of the increased basicity at the metal center, which shifted the redox potential toward more negative values and facilitated proton transfer, with the aspartate acting as a proton relay. Along these lines, Mahy and coworkers recently developed an artificial hydrogenase via insertion of maleimide functionalized CoPPIX into an artificial protein, named αRep (Table 2, entry 30) [237]. This protein has been shown to adopt a closed shell-like conformation hosting a large hydrophobic cavity, which has been previously exploited for the anchoring of catalytic metal complexes. Theoretical studies suggested that covalent anchoring of CoPPIX to αRep resulted in the complete encapsulation of the cofactor into the protein interior, with a glutamine sidechain acting as an axial ligand to the cobalt ion. The Co-αRep biohybrid showed hydrogen evolution photoinduced activity in the presence of [Ru(bpy)3]2+ as the photosensitizer and ascorbate as the electron donor, performing up to 163 TONs over 25 h of irradiation. As observed by Ghirlanda and coworkers, incorporation of CoPPIX into the protein scaffold significantly enhanced catalytic performance as compared to the lone porphyrin in the same experimental conditions. A substantial contribution to this field was granted by Bren and coworkers, who showcased the extraordinary catalytic potential of peptide–cobalt–porphyrin conjugates in light-driven hydrogen production [238,239]. On the basis of the promising results obtained in electrochemical catalysis, both Co-MP11 and Co-MC6*a were screened for photocatalytic activity in the presence of [Ru(bpy)3]2+ as a photosensitizer. In the photochemical system, Co-MP11 showed extended longevity, performing up to 2390 turnovers over 20 h of catalysis (Table 2, entry 31) [238]. This finding is in contrast with electrochemical studies, where substantial loss of activity was found after only 15 min of electrolysis. The short lifetime of Co-MP11 in electrocatalysis was attributed to porphyrin degradation by direct interaction with the electrode surface. Remarkably, Co-MC6*a outperformed Co-MP11 in terms of both catalyst durability (40 h) and turnover number (10,400), placing among the most efficient photocatalysts reported so far (Table 2, entry 32) [239]. In terms of catalyst longevity, Co-porphyrin-based artificial hydrogenases can be ordered as follows: Co-MC6*a > Co-αRep > Co-MP11 > Co-cyt b562 ≈ Co-Mb > CoP > CoTPPS. It can be noted that all the biohybrid systems display higher longevity with respect to the molecular metalloporphyrins CoP and CoTPPS, possibly resulting from a protective effect exerted by the protein matrix against light-induced catalyst degradation. This effect is much more evident when Co-porphyrin is covalently linked to the protein, as observed for Co-MP11, Co-αRep, and Co-MC6*a. Additionally, the higher durability of Co-MP11 compared to catalysts possessing a proper protein scaffold was attributed to tight histidine coordination in the covalently bound peptide [238]. Lastly, Co-MC6*a showed the highest longevity among Co-porphyrin biohybrids, despite its reduced dimensions. This was granted by the distal peptide chain which is designed to strongly interact with the porphyrin, thus acting as a shield against catalyst degradation [276]. An alternative approach for constructing nature-inspired photocatalytic systems has involved the use of the photosynthetic protein PSI in place of molecular photosensitizers [258]. Utschig and coworkers applied this strategy to develop biohybrid systems by combining cobalt or nickel containing molecular catalysts with PSI. Self-assembly of the cobaloxime catalyst Co(dmgH)2pyCl (dmgH = dimethylglyoximato, py = pyridine) with PSI provided the biohybrid system containing 2–4 cobalt catalyst units per PSI (Table 2, entry 33) [240]. This cobaloxime–PSI hybrid was capable of promoting hydrogen evolution with one of the highest rates (TOF = 170 min−1) reported for artificial photocatalysts, performing up to 5200 turnovers. However, significant loss of catalytic activity was observed after relatively short times (1.5 h), which was ascribed to cobaloxime dissociation from the protein. Subsequently, the DuBois nickel diphosphine complex was assembled with PSI either via direct interaction or delivery by apo flavodoxin (apoFld, Table 2, entries 34–35) [241]. The [Ni(P2PhN2Ph)2]-PSI biohybrid displayed a lower turnover frequency (44 min−1) and turnover number (1870) compared to the cobaloxime–PSI biohybrid. However, when the association to PSI was mediated by apoFld, the biohybrid reached 2865 turnovers in aqueous solution at nearly neutral pH, a value that is comparable to that displayed by the parent nickel molecular complex in acetonitrile in strongly acidic conditions [172]. Despite PSI-based biohybrids rapidly producing hydrogen, the large size of PSI (350 kDa) and the presence of multiple spectroscopically overlapping [4Fe-4S] cofactors hamper a thorough spectroscopic characterization and monitoring of the catalytic events. In subsequent work, the same authors designed a covalently assembled three component system by linking a cobaloxime catalyst and a ruthenium photosensitizer to Spinacia oleracea ferredoxin (Fd, Figure 9 and Table 2, entry 36) [242]. The small (10.5 kDa) electron transfer protein served as a scaffold to anchor the catalyst and the photosensitizer while also facilitating the charge transfer between them. Whereas the Ru–Fd–Co biohybrid performed up to 320 turnovers, no significant activity was observed when native Fd was replaced with apo-Fd, suggesting that the photoinduced electron transfer proceeds through the [2Fe2S] cofactor. Shaafat and coworkers designed an even smaller biohybrid system by covalently assembling a ruthenium photosensitizer ([RuII(2,2′-bipyridine)2(5,6-epoxy-5,6-dihydro- [1,10]-phenanthroline)]2+) to the previously described nickel-substituted rubredoxin (NiRd, Figure 10a and Table 2, entry 37) [243]. Attachment of the ruthenium complex to NiRd was provided through a native free cysteine residue, close to the metal-binding site. The effective distance between the nickel center and the photosensitizer was modulated by developing a series of mutants, featuring a free cysteine at different positions in the sequence. All the RuNiRd variants were capable of producing hydrogen upon light-irradiation, displaying enhanced activity with respect to the bimolecular Ru-complex/NiRd system. Furthermore, the activity of the biohybrid constructs was strongly dependent on the Ru–Ni distance calculated on the basis of rubredoxin crystal structure, suggesting an intramolecular electron-transfer pathway. Inspired by the nickel center of [NiFe]-hydrogenases, Chackraborty and coworkers recently developed artificial metalloenzymes by engineering a tetrathiolate nickel center (NiS4) into de novo designed coiled coils [244,245]. In their first study, self-assembly of two de novo peptides, each bearing a CXXC motif, in the presence of nickel(II) yielded a two-stranded coiled coil (2SCC) enclosing the NiS4 site (Figure 10b and Table 2, entry 38) [244]. The artificial enzyme was able to produce H2 under photocatalytic conditions with a bell-shaped pH dependence, reaching its maximum activity at pH 5.6. Spectroscopic studied revealed that a precise acidity is required for activity, suggesting that the protonation state of Cys is crucial for H2 production. In the latest study, a peptide sequence was designed to favor a four-stranded coiled coil (4SCC), and the presence of a rubredoxin-like CXXC motif in each peptide strand was conceived to provide a dual nickel-binding site (Table 2, entry 39) [245]. Interestingly, nickel binding induced an oligomerization shift, stabilizing the formation of dimers, while the apo form was predominantly trimeric. This behavior indicates that nickel induced dissociation of the assembly by imposing its coordination preferences, which were stabilized in the 2SCC, as found in previous work. Notably, the newly designed Ni-4SCC peptide displayed a slightly improved activity compared to the 2SCC analogue. Even though these de novo metalloenzymes display relatively low turnover numbers, they represent a significant expansion of the de novo design approach for the development of artificial metalloproteins able to catalyze energy-relevant processes.

In the last few decades, research interest toward hydrogen production has greatly exploded, prompted by the urgency to face the current energy crisis. Hydrogen has a huge potential as a clean fuel, due to its incredibly high energy density and the possibility to be combusted without generating CO2. However, the development of cost-effective and environmentally friendly methods for hydrogen production is far from trivial, and much effort has been devoted to the elaboration of enzyme-based devices. These typically involve the use of hydrogenases either for light-driven hydrogen production (photoelectrochemical cells) or its conversion into electricity, upon matching of hydrogen oxidation with a reductive chemical reaction (biofuel cells). Despite hydrogenases showing impressive catalytic properties, their application in large-scale processes is currently hampered by difficulties with their expression and their limited tolerance to ambient oxygen. Extensive research has been conducted to elucidate the catalytic mechanism of hydrogenases, with the aim of replicating their function into artificial and tailorable systems. Since the late 1990s, a plethora of molecular catalysts have been proposed as hydrogenase mimics. While some have focused on reproducing the structure of the [FeFe]-hydrogenase active site, a wide variety of HECs have also been developed, combining redox non-innocent ligands with nickel or cobalt ions. Considerable improvements in catalytic properties have been accomplished through the incorporation of secondary coordination sphere interactions and, particularly, proton shuttle functionalities. Even though remarkably high catalytic performances have been achieved in some cases, the majority of these complexes only function in organic solvents and in the presence of strong acids, further highlighting the essential role of a bioinspired ligand framework in imparting enzyme-like properties. Incorporation of catalytic centers into natural or designed protein scaffolds conceivably represents the most promising strategy to obtain artificial hydrogenases capable of paralleling their natural counterparts. Indeed, protein design methods allow precisely controlling the dielectric properties and the interactions within the residues in the active site, enabling a fine modulation of the metal redox properties. Furthermore, protein scaffolds can also serve as a platform to bind together catalytic and light-harvesting units, as reported in prominent examples discussed within the review. In summary, biohybrid and bioinspired catalysts are promising candidates for sustainable hydrogen and energy production; nevertheless, some important challenges still need to be addressed. In particular, most artificial biomolecular catalysts are still fragile and display limited durability under operative conditions, especially in light-driven catalysis. One possible approach to overcome this problem is the immobilization of catalysts onto nanomaterials, which could extend their durability, allowing for easy catalyst recycling. Furthermore, as the majority of HECs reported so far contain transition metal ions such as nickel and cobalt, evaluation of their toxicity is a crucial aspect to consider, especially when scaling up the process to the industrial scale [277]. Indeed, the toxicity of metal complexes greatly depends on the nature of the ligands; thus, it is not trivial to predict [278]. A step further in the development of sustainable methods for hydrogen production would consist of the construction of completely protein or peptide-based materials, opening the way for the assembly of biosynthetic nanoreactors for energy-related catalysis.