Experimental and DFT Investigations Reveal the Influence of the Outer Coordination Sphere on the Vibrational Spectra of Nickel-Substituted Rubredoxin, a Model Hydrogenase Enzyme

Targeted Modulation of Photocatalytic Hydrogen Evolution Activity by Nickel-Substituted Rubredoxin through Functionalized Ruthenium Phototriggers

Light-driven hydrogen evolution is a promising means of sustainable energy production to meet global energy demand. This study investigates the photocatalytic hydrogen evolution activity of nickel-substituted rubredoxin (NiRd), an artificial hydrogenase mimic, covalently attached to a ruthenium phototrigger (RuNiRd). By systematically modifying the para-substituents on Ru(II) polypyridyl complexes, we sought to optimize the intramolecular electron transfer processes within the RuNiRd system. A series of electron-donating and electron-withdrawing groups were introduced to tune the photophysical, photochemical, and electrochemical properties of the ruthenium complexes. Our findings reveal that electron-donating substituents can increase the hydrogen evolution capabilities of the artificial enzyme to a point; however, the complexes with the most electron-donating substituents suffer from short lifetimes and inefficient reductive quenching, rendering them inactive. The present work highlights the intricate balance required between driving force, lifetime, and quenching efficiency for effective light-driven catalysis, providing valuable insights into the design of artificial enzyme-photosensitizer constructs.

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

Electrocatalytic Hydrogen Evolution using a Nickel-based Calixpyrrole Complex: Controlling the Secondary Coordination Sphere on an Electrode Surface

Incorporating design elements from homogeneous catalysts to construct well defined active sites on electrode surfaces is a promising approach for developing next generation electrocatalysts for energy conversion reactions. Furthermore, if functionalities that control the electrode microenvironment could be integrated into these active sites it would be particularly appealing. In this context, a square planar nickel calixpyrrole complex, Ni(DPMDA) (DPMDA=2,2'-((diphenylmethylene)bis(1H-pyrrole-5,2-diyl))bis(methaneylylidene))bis(azaneylylidene))dianiline) with pendant amine groups is reported that forms a heterogeneous hydrogen evolution catalyst using anilinium tetrafluoroborate as the proton source. The supported Ni(DPMDA) catalyst was surprisingly stable and displayed fast reaction kinetics with turnover frequencies (TOF) up to 25,900 s-1 or 366,000 s-1  cm-2 . Kinetic isotope effect (KIE) studies revealed a KIE of 5.7, and this data, combined with Tafel slope analysis, suggested that a proton-coupled electron transfer (PCET) process involving the pendant amine groups was rate-limiting. While evidence of an outer-sphere reduction of the Ni(DPMDA) catalyst was observed, it is hypothesized that the control over the secondary coordination sphere provided by the pendant amines facilitated such high TOFs and enabled the PCET mechanism. The results reported herein provide insight into heterogeneous catalyst design and approaches for controlling the secondary coordination sphere on electrode surfaces.

Rubredoxin Protein Scaffolds Sourced from Diverse Environmental Niches as an Artificial Hydrogenase Platform

Nickel-substituted rubredoxin (NiRd) from Desulfovibrio desulfuricans has previously been shown to act as both a structural and functional mimic of the [NiFe] hydrogenase. However, improvements both in turnover frequency and overpotential are needed to rival the native [NiFe] hydrogenase enzymes. Characterization of a library of NiRd mutants with variations in the secondary coordination sphere suggested that protein dynamics played a substantial role in modulating activity. In this work, rubredoxin scaffolds were selected from diverse organisms to study the effects of distal sequence variation on catalytic activity. It was found that though electrochemical catalytic activity was only slightly impacted across the series, the Rd sequence from a psychrophilic organism exhibited substantially higher levels of solution-phase hydrogen production. Additionally, Eyring analyses suggest that catalytic activation properties relate to the growth temperature of the parent organism, implying that the general correlation between the parent organism environment and catalytic activity often seen in naturally occurring enzymes may also be observed in artificial enzymes. Selecting protein scaffolds from hosts that inhabit diverse environments, particularly low-temperature environments, represents an alternative approach for engineering artificial metalloenzymes.

Protein-based models offer mechanistic insight into complex nickel metalloenzymes

There are ten nickel enzymes found across biological systems, each with a distinct active site and reactivity that spans reductive, oxidative, and redox-neutral processes. We focus on the reductive enzymes, which catalyze reactions that are highly germane to the modern-day climate crisis: [NiFe] hydrogenase, carbon monoxide dehydrogenase, acetyl coenzyme A synthase, and methyl coenzyme M reductase. The current mechanistic understanding of each enzyme system is reviewed along with existing knowledge gaps, which are addressed through the development of protein-derived models, as described here. This opinion is intended to highlight the advantages of using robust protein scaffolds for modeling multiscale contributions to reactivity and inspire the development of novel artificial metalloenzymes for other small molecule transformations.

The underappreciated influence of ancillary halide on metal–ligand proton tautomerism

Syntheses of Vaska-type complexes [IrP₂X(CO)] (P = phosphine, X = halide) with all four common halides (fluoride, chloride, bromide, and iodide) was attempted using a protic and hemilabile imidazolyl di-tert-butyl phosphine ligand. In the solid-state, all four complexes were found to be ionic with the halides in the outer-sphere, and the fourth coordination site of the square plane occupied by the imidazole arm of the ligand. In solution, however, the chloride complex was found to be in equilibrium with an octahedral Ir^III–H species at room temperature. For the bromide and iodide analogs, the corresponding Ir^III–H species were also observed but only after heating the solutions. The neutral Ir^I Vaska's analogs for X = Cl, Br, and I were obtained upon addition of excess halide salt, albeit heating was required for X = Br and I. The Ir^III–H species are proposed to originate from tautomerization of minor amounts of the electron rich neutral Vaska analog (halide inner-sphere and phosphines monodentate) that are in equilibrium with the ionic species. Heating is required for the larger anions of bromide and iodide to overcome a kinetic barrier associated with their movement to an inner-sphere position prior to tautormerization. For the fluoride analog, the Ir^III–H was not observed, attributable to strong hydrogen bonding interactions of the imidazolyl proton with the fluoride anion.

Ligand protonated Ir^I bisphosphine carbonyl complexes isolated as halide salts equilibrate with their neutral Ir^III–H congeners in solution. The equilibrium constant and energy barrier to interconversion are dependent on the identity of the halide.

Exploring Structure and Function of Redox Intermediates in [NiFe]-Hydrogenases by an Advanced Experimental Approach for Solvated, Lyophilized and Crystallized Metalloenzymes

To study metalloenzymes in detail, we developed a new experimental setup allowing the controlled preparation of catalytic intermediates for characterization by various spectroscopic techniques. The in situ monitoring of redox transitions by infrared spectroscopy in enzyme lyophilizate, crystals, and solution during gas exchange in a wide temperature range can be accomplished as well. Two O2 -tolerant [NiFe]-hydrogenases were investigated as model systems. First, we utilized our platform to prepare highly concentrated hydrogenase lyophilizate in a paramagnetic state harboring a bridging hydride. This procedure proved beneficial for 57 Fe nuclear resonance vibrational spectroscopy and revealed, in combination with density functional theory calculations, the vibrational fingerprint of this catalytic intermediate. The same in situ IR setup, combined with resonance Raman spectroscopy, provided detailed insights into the redox chemistry of enzyme crystals, underlining the general necessity to complement X-ray crystallographic data with spectroscopic analyses.

Metal/Ligand Proton Tautomerism Facilitates Dinuclear H2 Reductive Elimination

Using the doubly protic bis-pyrazole-pyridine ligand (N(NN^H)₂), we have synthesized an octahedral Ir^III-H [HIr(κ³-N(NN^H)(NN^-))(CO)(^tBuPy)]⁺ ([1-MH]⁺) from an Ir^I starting material. This hydride was generated by adding sufficient electron density to the metal center such that it became the thermodynamically preferred site of protonation. It was observed via UV-vis spectroscopy that [1-MH]⁺ establishes a [^tBuPy] dependent equilibrium with a ligand protonated square-planar Ir^I [Ir(N(NN^H)₂)(CO)]⁺ ([2-LH]⁺). This example of metal/ligand proton tautomerism is unusual in that the position of the equilibrium can be controlled by the concentration of exogeneous ligand (i.e., ^tBuPy). This equilibrium was shown to be key to the reactivity of the Ir^III-H; 2 equiv of [1-MH]⁺ release H₂, converting to the Ir^II dimer [[Ir(N(NN^-)(NN^H))(CO)(^tBuPy)]₂]²⁺ ([7]²⁺) under mild conditions (observable at room temperature). Mechanistic evidence is presented to support that this dinuclear reductive elimination occurs by tautomerization of the metal hydride [1-MH]⁺ to a ligand protonated species [1-LH]⁺, from which ligand dissociation is facile, generating [2-LH]⁺. Subsequent reaction of [2-LH]⁺ with [1-MH]⁺ allows for production of H₂ and the Ir^II dimer [7]²⁺. The tautomerization between the metal-hydride and the ligand protonated species provides a low energy pathway for ligand dissociation, opening the needed coordination site. The ability to control the interconversion between a metal-hydride and a ligand-protonated congener using an exogeneous ligand introduces a new strategy for catalyst design with proton responsive ligands.

Light-driven catalysis with engineered enzymes and biomimetic systems

Efforts to drive catalytic reactions with light, inspired by natural processes like photosynthesis, have a long history and have seen significant recent growth. Successfully engineering systems using biomolecular and bioinspired catalysts to carry out light-driven chemical reactions capitalizes on advantages offered from the fields of biocatalysis and photocatalysis. In particular, driving reactions under mild conditions and in water, in which enzymes are operative, using sunlight as a renewable energy source yield environmentally friendly systems. Furthermore, using enzymes and bioinspired systems can take advantage of the high efficiency and specificity of biocatalysts. There are many challenges to overcome to fully capitalize on the potential of light-driven biocatalysis. In this mini-review, we discuss examples of enzymes and engineered biomolecular catalysts that are activated via electron transfer from a photosensitizer in a photocatalytic system. We place an emphasis on selected forefront chemical reactions of high interest, including CH oxidation, proton reduction, water oxidation, CO2 reduction, and N2 reduction.

H2 Evolution from a Thiolate-Bound Ni(III) Hydride

Terminal Ni^III hydrides are proposed intermediates in proton reduction catalyzed by both molecular electrocatalysts and metalloenzymes, but well-defined examples of paramagnetic nickel hydride complexes are largely limited to bridging hydrides. Herein, we report the synthesis of an S = 1/2, terminally bound thiolate-Ni^III-H complex. This species and its terminal hydride ligand in particular have been thoroughly characterized by vibrational and EPR techniques, including pulse EPR studies. Corresponding DFT calculations suggest appreciable spin leakage onto the thiolate ligand. The hyperfine coupling to the terminal hydride ligand of the thiolate-Ni^III-H species is comparable to that of the hydride ligand proposed for the Ni-C hydrogenase intermediate (Ni^III-H-Fe^II). Upon warming, the featured thiolate-Ni^III-H species undergoes bimolecular reductive elimination of H₂. Associated kinetic studies are discussed and compared with a structurally related Fe^III-H species that has also recently been reported to undergo bimolecular H-H coupling.

On the Nature of C(sp3)-C(sp2) Bond Formation in Nickel-Catalyzed Tertiary Radical Cross-Couplings: A Case Study of Ni/Photoredox Catalytic Cross-Coupling of Alkyl Radicals and Aryl Halides

The merger of photoredox and nickel catalysis has enabled the construction of quaternary centers. However, the mechanism, role of the ligand, and effect of the spin state for this transformation and related Ni-catalyzed cross-couplings involving tertiary alkyl radicals in combination with bipyridine and diketonate ligands remain unknown. Several mechanisms have been proposed, all invoking a key Ni(III) species prior to undergoing irreversible inner-sphere reductive elimination. In this work, we have used open-shell dispersion-corrected DFT calculations, quasi-classical dynamics calculations, and experiments to study in detail the mechanism of carbon-carbon bond formation in Ni bipyridine- and diketonate-based catalytic systems. These calculations revealed that access to high spin states (e.g., triplet spin state tetrahedral Ni(II) species) is critical for effective radical cross-coupling of tertiary alkyl radicals. Further, these calculations revealed a disparate mechanism for the C-C bond formation. Specifically, contrary to the neutral Ni-bipyridyl system, diketonate ligands lead directly to the corresponding tertiary radical cross-coupling products via an outer-sphere reductive elimination step via triplet spin state from the Ni(III) intermediates. Implications to related Ni-catalyzed radical cross-couplings and the design of new transformations are discussed.

The oxygen reactivity of an artificial hydrogenase designed in a reengineered copper storage protein

The O₂ reactivity of an artificial biomolecular hydrogenase, the nickel binding protein (NBP) is investigated.

Quantitative prediction of electronic absorption spectra of copper(II)-bioligand systems: Validation and applications

The visible region of the electronic absorption spectra of Cu(II) complexes was studied by time-dependent density functional theory (TD-DFT). The performance of twelve functionals in the prediction of absorption maxima (λ_max) was tested on eleven compounds with different geometry, donors and charge. The ranking of the functionals for λ_max was determined in terms of mean absolute percent deviation (MAPD) and standard deviation (SD) and it is as follows: BHandHLYP > M06 ≫ CAM-B3LYP ≫ MPW1PW91 ~ B1LYP ~ BLYP > HSE06 ~ B3LYP > B3P86 ~ ω-B97x-D ≫ TPSSh ≫ M06-2X (MAPD) and BHandHLYP > M06 ~ HSE06 > ω-B97x-D ~ CAM-B3LYP ~ MPW1PW91 > B1LYP ~ B3LYP > B3P86 > BLYP ≫ TPSSh ≫ M06-2X (SD). With BHandHLYP functional the MAPD is 3.1% and SD is 2.3%, while with M06 the MAPD is 3.7% and SD is 3.7%. The protocol validated in the first step of the study was applied to: i) calculate the number of transitions in the spectra and relate them to the geometry of Cu(II) species; ii) determine the coordination of axial water(s); iii) predict the electronic spectra of the systems where Cu(II) is bound to human serum albumin (HSA) and to the regions 94-97 and 108-112 of prion protein (PrP). The results indicate that the proposed computational protocol allows a successful prediction of the electronic spectra of Cu(II) species and to relate an experimental spectrum to a specific structure.

Redesign of a Copper Storage Protein into an Artificial Hydrogenase

We report the construction of an artificial hydrogenase (ArH) by reengineering a Cu storage protein (Cspl) into a Ni-binding protein (NBP) employing rational metalloprotein design. The hypothesis driven design approach involved deleting existing Cu sites of Csp1 and identification of a target tetrathiolate Ni binding site within the protein scaffold followed by repacking the hydrophobic core. Guided by modeling, the NBP was expressed and purified in high purity. NBP is a well-folded and stable construct displaying native-like unfolding behavior. Spectroscopic and computational studies indicated that the NBP bound nickel in a distorted square planar geometry that validated the design. Ni(II)-NBP is active for photo-induced H₂ evolution following a reductive quenching mechanism. Ni(II)-NBP catalyzed H⁺ reduction to H₂ gas electrochemically as well. Analysis of the catalytic voltammograms established a proton-coupled electron transfer (PCET) mechanism. Electrolysis studies confirmed H₂ evolution with quantitative Faradaic yields. Our studies demonstrate an important scope of rational metalloprotein design that allows imparting functions into protein scaffolds that have natively not evolved to possess the same function of the target metalloprotein constructs.

Understanding and Modulating Metalloenzymes with Unnatural Amino Acids, Non-Native Metal Ions, and Non-Native Metallocofactors

Metalloproteins set the gold standard for performing important functions, including catalyzing demanding reactions under mild conditions. Designing artificial metalloenzymes (ArMs) to catalyze abiological reactions has been a major endeavor for many years, but most ArM activities are far below those of native enzymes, making them unsuitable for most pratical applications. A critical step to advance the field is to fundamentally understand what it takes to not only confer but also fine-tune ArM activities so they match those of native enzymes. Indeed, only once we can freely modulate ArM activity to rival (or surpass!) natural enzymes can the potential of ArMs be fully realized. A key to unlocking ArM potential is the observation that one metal primary coordination sphere can display a range of functions and levels of activity, leading to the realization that secondary coordination sphere (SCS) interactions are critically important. However, SCS interactions are numerous, long-range, and weak, making them very difficult to reproduce in ArMs. Furthermore, natural enzymes are tied to a small set of biologically available functional moieties from canonical amino acids and physiologically available metal ions and metallocofactors, severely limiting the chemical space available to probe and tune ArMs. In this Account, we summarize the use of unnatural amino acids (UAAs) and non-native metal ions and metallocofactors by our group and our collaborators to probe and modulate ArM functions. We incorporated isostructural UAAs in a type 1 copper (T1Cu) protein azurin to provide conclusive evidence that axial ligand hydrophobicity is a major determinant of T1Cu redunction potential ( E°'). Closely related work from other groups are also discussed. We also probed the role of protein backbone interactions that cannot be altered by standard mutagenesis by replacing the peptide bond with an ester linkage. We used insight gained from these studies to tune the E°' of azurin across the entire physiological range, the broadest range ever achieved in a single metalloprotein. Introducing UAA analogues of Tyr into ArM models of heme-copper oxidase (HCO) revealed a linear relationship between p K_a, E°', and activity. We also substituted non-native hemes and non-native metal ions for their native equivalents in these models to resolve several issues that were intractable in native HCOs and the closely related nitric oxide reductases, such as their roles in modulating substrate affinity, electron transfer rate, and activity. We incorporated abiological cofactors such as ferrocene and Mn(salen) into azurin and myoglobin, respectively, to stabilize these inorganic and organometallic compounds in water, confer abiological functions, tune their E°' and activity through SCS interactions, and show that the approach to metallocofactor anchoring and orientation can tune enantioselectivity and alter function. Replacing Cu in azurin with non-native Fe or Ni can impart novel activities, such as superoxide reduction and C-C bond formation. While progress was made, we have identified only a small fraction of the interactions that can be generally applied to ArMs to fine-tune their functions. Because SCS interactions are subtle and heavily interconnected, it has been difficult to characterize their effects quantitatively. It is vital to develop spectroscopic and computational techniques to detect and quantify their effects in both resting states and catalytic intermediates.

Driving Protein Conformational Changes with Light: Photoinduced Structural Rearrangement in a Heterobimetallic Oxidase

The heterobimetallic R2lox protein binds both manganese and iron ions in a site-selective fashion and activates oxygen, ultimately performing C-H bond oxidation to generate a tyrosine-valine cross-link near the active site. In this work, we demonstrate that, following assembly, R2lox undergoes photoinduced changes to the active site geometry and metal coordination motif. Through spectroscopic, structural, and mass spectrometric characterization, the photoconverted species is found to consist of a tyrosinate-bound iron center following light-induced decarboxylation of a coordinating glutamate residue and cleavage of the tyrosine-valine cross-link. This process occurs with high quantum efficiencies (Φ = 3%) using violet and near-ultraviolet light, suggesting that the photodecarboxylation is initiated via ligand-to-metal charge transfer excitation. Site-directed mutagenesis and structural analysis suggest that the cross-linked tyrosine-162 is the coordinating residue. One primary product is observed following irradiation, indicating potential use of this class of proteins, which contains a putative substrate channel, for controlled photoinduced decarboxylation processes, with relevance for in vivo functionality of R2lox as well as application in environmental remediation.

Spectroelectrochemical investigations of nickel cyclam indicate different reaction mechanisms for electrocatalytic CO2 and H+ reduction

Characterization of a Ni^III species during reductive catalysis by [Ni(cyclam)]²⁺ implicates an ECCE mechanism for hydrogen production in aqueous solution.

Light-Driven Hydrogen Evolution by Nickel-Substituted Rubredoxin

An enzymatic system for light-driven hydrogen generation has been developed through covalent attachment of a ruthenium chromophore to nickel-substituted rubredoxin (NiRd). The photoinduced activity of the hybrid enzyme is significantly greater than that of a two-component system and is strongly dependent on the position of the ruthenium phototrigger relative to the active site, indicating a role for intramolecular electron transfer in catalysis. Steady-state and time-resolved emission spectra reveal a pathway for rapid, direct quenching of the ruthenium excited state by nickel, but low overall turnover numbers suggest initial electron transfer is not the rate-limiting step. This approach is ideally suited for detailed mechanistic investigations of catalysis by NiRd and other molecular systems, with implications for generation of solar fuels.