The Golden Age of Transfer Hydrogenation

Artificial photosynthesis: molecular systems for catalytic water oxidation.

Photosynthesis converts light energy into biologically useful chemical energy vital to life on Earth. The initial reaction of photosynthesis takes place in photosystem II (PSII), a 700-kilodalton homodimeric membrane protein complex which catalyses photo-oxidation of water into dioxygen through an S-state cycle of the oxygen evolving complex (OEC). The structure of PSII has been solved by X-ray diffraction (XRD) at 1.9-angstrom (A) resolution, which revealed that the OEC is a Mn4CaO5-cluster coordinated by a well-defined protein environment1. However, extended X-ray absorption fine structure (EXAFS) studies showed that the manganese cations in the OEC are easily reduced by X-ray irradiation2, and slight differences were found in the Mn–Mn distances between the results of XRD1, EXAFS3–7 and theoretical studies8–14. Here we report a ‘radiation-damage-free’ structure of PSII from Thermosynechococcus vulcanus in the S1 state at a resolution of 1.95 A using femtosecond X-ray pulses of the SPring-8 angstrom compact free-electron laser (SACLA) and a huge number of large, highly isomorphous PSII crystals. Compared with the structure from XRD, the OEC in the X-ray free electron laser structure has Mn–Mn distances that are shorter by 0.1–0.2 A. The valences of each manganese atom were tentatively assigned as Mn1D(III), Mn2C(IV), Mn3B(IV) and Mn4A(III), based on the average Mn–ligand distances and analysis of the Jahn–Teller axis on Mn(III). One of the oxo-bridged oxygens, O5, has significantly longer Mn–O distances in contrast to the other oxo-oxygen atoms, suggesting that it is a hydroxide ion instead of a normal oxygen dianion and therefore may serve as one of the substrate oxygen atoms. These findings provide a structural basis for the mechanism of oxygen evolution, and we expect that this structure will provide a blueprint for design of artificial catalysts for water oxidation.

https://ousar.lib.okayama-u.ac.jp/files/public/5/53637/20160528122001540368/Nature-Shen manuscript for deposit(151002).pdf

Native structure of photosystem II at 1.95 Å resolution viewed by femtosecond X-ray pulses.

Dialkylbiaryl phosphines are a valuable class of ligand for Pd-catalyzed amination reactions and have been applied in a range of contexts. This perspective attempts to aid the reader in the selection of the best choice of reaction conditions and ligand of this class for the most commonly encountered and practically important substrate combinations.

Dialkylbiaryl phosphines in Pd-catalyzed amination: a user's guide

Because of its fundamental importance in many branches of science, hydrogen bonding is a subject of intense contemporary research interest. The physical and chemical properties of hydrogen bonds in the ground state have been widely studied both experimentally and theoretically by chemists, physicists, and biologists. However, hydrogen bonding in the electronic excited state, which plays an important role in many photophysical processes and photochemical reactions, has scarcely been investigated.Upon electronic excitation of hydrogen-bonded systems by light, the hydrogen donor and acceptor molecules must reorganize in the electronic excited state because of the significant charge distribution difference between the different electronic states. The electronic excited-state hydrogen-bonding dynamics, which are predominantly determined by the vibrational motions of the hydrogen donor and acceptor groups, generally occur on ultrafast time scales of hundreds of femtoseconds. As a result, state-of-the-art femtos...

/pdf/hydrogen-bonding-in-the-electronic-excited-state-1wj02r16hb.pdf

Hydrogen Bonding in the Electronic Excited State

ing a hydrogen atom from substrates with relatively weak X−H bonds (Figure 56). Such a situation was found in the catalytic cycle of ACCO, where cleavage of the N−H bond is facilitated by the ferrous ion capable of (partly) reducing the resulting N-based radical to an anion. In the catalytic cycle of HEPD and in the reaction of HppE with an enantiomer of the native substrate, a C−H bond at the carbon hosting a deprotonated alcohol group is severed by the superoxide. In these cases transfer of a hydrogen atom to the superoxide is coupled to one-electron reduction of Fe(III) to Fe(II), which yields aldehyde or ketone products. In a similar way a thioaldehyde is produced during the initial steps of the IPNS catalytic cycle. Figure 53. Reaction catalyzed by NDO. Figure 54. X-ray structure of the non-heme iron cofactor in the NDO active site (PDB 1O7G). Figure 55. Reaction mechanism for the dihydroxylation reaction catalyzed by NDO suggested on the basis of the results of a DFT study. Relative energy values in kilocalories per mole are given beneath the structures shown and above the arrows for the transition states connecting them. Chemical Reviews Review dx.doi.org/10.1021/cr400388t | Chem. Rev. 2014, 114, 3601−3658 3636 The most often encountered reaction of the Fe(III)−O2 species is an electrophilic attack on an electron-rich (co)substrate that yields an Fe(II) intermediate with a peroxide bridge between the ferrous ion and an organic molecule (Figure 57). In all examples shown in the figure the reactive spin state is a quintet with high-spin Fe(III) antiferromagnetically coupled to the superoxide radical. This particular electronic structure allows for a smooth two-electron reaction, leading to a ferrous intermediate featuring a high-spin Fe(II) ion, i.e., also lying on a quintet PES. When the substrates are easily one-electron-oxidized, as, for example, carotenoids or catecholates, already the ternary enzyme−dioxygen−organic substrate complex may contain an organic radical along with the superoxide anion stabilized by coordination to the metal (Figure 58). With proper (antiferromagnetic) alignment of the spins of the unpaired electrons on the two radicals, direct coupling between them proceeds with a straightforward formation of a new C−O bond. With the O2 ligand already protonated, several different reaction scenarios are usually plausible, and one of them, arguably the simplest, involves a transfer of the HOO group from the metal ion to the organic substrate (Figure 59). In the reaction of HEPD, it is the distal, i.e., originally protonated, oxygen atom that attacks the aldehyde group and the proton is transferred to the second (proximal) oxygen atom with the help of a phosphonic group of the substrate. In this reaction HOO acts as a nucleophile. In the case of HGD, the HOO ligand has a partial radical character, and hence, in its attack on the aromatic ring, it behaves as an electrophilic reagent. Moreover, in the HGD and ACO cases, it is the proximal oxygen atom that is directly transferred from the metal ion to the organic radical, i.e., the proton remains on the same oxygen atom throughout the reaction. When a catalytic reaction involves a homolytic O−O bond cleavage, one end of the peroxo group is in contact with the metal ion, and it is reduced to HO− or RO− when the O−O bond breaks (Figure 60). The electron required for the reduction is provided usually by the ferrous ion; in intradiol dioxygenases, which host Fe(III) in the active site, a tyrosinate ligand can serve as a reductant. Heterolytic cleavage of the O−O bond typically yields highvalent iron(IV)−oxo species, and such a reaction requires Fe(II), a deprotonated proximal oxygen atom, and usually also that the distal oxygen atom has a chance to develop a second covalent bond when the O−O gets broken (Figure 61). The bonding partner for the distal oxygen can be a hydrogen (PDHs, IPNS, ACCO) or a (co)substrate’s carbon (αKAOs, ACO, Dke1) atom. Heterolysis to Fe(IV)O usually proceeds on the quintet PES. Finally, heterolytic O−O bond cleavage may proceed without changing the oxidation state of the metal but instead with coupled two-electron oxidation of an organic substrate (Figure Figure 56. Examples of Fe(II)/Fe(III)−O2 reactions involving hydrogen atom abstraction. Only the most relevant fragment of the substrate is shown. Figure 57. Examples of Fe(III)−O2 reactions involving electrophilic attack and two-electron oxidation of an organic (co)substrate. Figure 58. Examples of Fe(II)/Fe(III)−O2 reactions involving radical coupling between the superoxide and an organic radical. Only the most relevant fragment of the substrate is shown. Chemical Reviews Review dx.doi.org/10.1021/cr400388t | Chem. Rev. 2014, 114, 3601−3658 3637 62). In both cases depicted in the figure, a hydroperoxo group is bound to high-spin Fe(III), and when theO−Obond cleaves, the OH group remains on iron and the other oxygen atom forms two covalent bonds with the organic substrate. This kind of heterolytic O−O cleavage requires that the donor orbital of the organic substrate, i.e., the one that provides two electrons for reduction of the peroxo group, has a good overlap with the O−O σ* orbital. 5.2. Dinuclear Non-Heme Iron Enzymes 5.2.1. Methane Monooxygenase. MMO is an enzyme which inserts one oxygen from O2 into methane to form methanol. The active site of the soluble form of MMO is shown in Figure 63. It contains an iron dimer complex linked by oxygen-derived ligands and has four glutamates and two histidines. Due to the presence of a well-resolved X-ray structure, and the technical importance of the reaction catalyzed, this was actually the first redox-active enzymemechanism that was treated with the cluster model using modern DFT functionals in 1997. Several groups were active at an early stage, and the most essential parts of the reaction mechanism were determined more than a decade ago. A comprehensive review of this development was written by Friesner et al. In short, the active species (compound Q) has a diamond core structure with two bridging oxo groups and is in an Fe2(IV,IV) state, one of themost oxidized species in nature. One of the oxo groups of compound Q activates nethane by an abstraction of one hydrogen atom. The TS is linear in C···H···O to form an Fe2(III,IV) state and amethyl radical. The loss of entropy at the TS is a large part of the barrier. In the second step, the methyl radical recombines with the bridging hydroxide formed in the first step. The TS for this step was first located by Basch et al. This rebound mechanism was criticized by interpretations of radical clock experiments, which appeared to show that there could not be a sufficiently long-lived alkyl radical to be consistent with a two-step mechanism. A concerted mechanism with simultaneous cleavage of the C−H bond and formation of the C−O bond was therefore suggested. Several explanations were suggested to resolve this apparent discrepancy between experiment and theory. In one of them, it was concluded that the radical clock probe molecules Figure 59. Examples of HOO transfer reactions. Only the most relevant fragment of the substrate is shown. Figure 60. Examples of O−O bond homolysis facilitated by oneelectron reduction of the proximal oxygen atom. Only the most relevant fragment of the substrate is shown. Figure 61. Examples of heterolytic O−O bond cleavage yielding iron(IV)−oxo species. Figure 62. Examples of heterolytic O−Obond cleavage proceeding with direct two-electron oxidation of the organic molecule. Chemical Reviews Review dx.doi.org/10.1021/cr400388t | Chem. Rev. 2014, 114, 3601−3658 3638 were chemically so different from methane that different mechanisms were likely. The probes are much easier to ionize and could form cations instead of radicals. Another suggestion was a so-called two-channel mechanism of the dynamics involving a bound radical intermediate. A third possibility could be something analogous to the two-state reactivity mechanism suggested for P450 to explain similar discrepancies in that case, but this has not been tested. In summary, a concerted mechanism, as suggested by the experiments, has never been found in DFT modeling calculations for methane hydroxylation by MMO, at least when reasonable models have been tried, and the suggestion has therefore not survived. In the initial phase of the MMO studies, there were problems converging to proper electronic states. These problems were solved by Friesner et al., who were able to obtain the correct antiferromagnetic coupling of compound Q with two high-spin irons. A state of key importance is also the first intermediate after Q, with an electronic structure characterized as Fe2(III,IV)− O•, discussed further below in connection with mixedMn−Fe dimers. It was found that already at the TS for hydrogen abstraction the iron dimer is in an Fe2(III,IV) state as it is in the product of this step. The bridging oxygen radical would then act as a hydrogen atom abstractor. At the TS, the spins are divided between a bridging oxo ligand and the methyl, while the iron spins stay essentially constant from the Fe2(III,IV)−O state to the product. The O−O bond cleavage to reach compound Q is also a significant step in the catalytic cycle. Again, several groups were involved in studying this step at an early stage. There was essential agreement among these studies on the mechanism. First, a peroxide (compound P) is formed between the two irons in an Fe2(III,III) state. Several different structures of P are nearly degenerate. In the TS for the O−O cleavage, the oxygens are symmetric, but only one of the irons is redox-active. In the final Figure 63. Active site of methane monooxygenase. Figure 64. R1 and R2 proteins in RNR. Chemical Reviews Review dx.doi.org/10.1021/cr400388t | Chem. Rev. 2014, 114, 3601−3658 3639 stage, the other iron also changes its oxidation state to IV, and compound Q is formed. In the study by Friesner et al., 

Quantum chemical studies of mechanisms for metalloenzymes.

Electrocatalytic water oxidation using the oxidatively robust 2,7-[bis(2-pyridylmethyl)aminomethyl]-1,8-naphthyridine ligand (BPMAN)-based dinuclear copper(II) complex, [Cu2(BPMAN)(μ-OH)](3+), has been investigated. This catalyst exhibits high reactivity and stability towards water oxidation in neutral aqueous solutions. DFT calculations suggest that the O-O bond formation takes place by an intramolecular direct coupling mechanism rather than by a nucleophilic attack of water on the high-oxidation-state Cu(IV)=O moiety.

https://onlinelibrary.wiley.com/doi/pdf/10.1002/anie.201411625

Electrocatalytic water oxidation by a dinuclear copper complex in a neutral aqueous solution.

We report systematic quantum mechanics-only (QM-only) and QM/molecular mechanics (MM) calculations on an enzyme-catalyzed reaction to assess the convergence behavior of QM-only and QM/MM energies with respect to the size of the chosen QM region. The QM and MM parts are described by density functional theory (typically B3LYP/def2-SVP) and the CHARMM force field, respectively. Extending our previous work on acetylene hydratase with QM regions up to 157 atoms (Liao and Thiel, J. Chem. Theory Comput. 2012, 8, 3793), we performed QM/MM geometry optimizations with a QM region M4 composed of 408 atoms, as well as further QM/MM single-point calculations with even larger QM regions up to 657 atoms. A charge deletion analysis was conducted for the previously used QM/MM model (M3a, with a QM region of 157 atoms) to identify all MM residues with strong electrostatic contributions to the reaction energetics (typically more than 2 kcal/mol), which were then included in M4. QM/MM calculations with this large QM region M4 lead to the same overall mechanism as the previous QM/MM calculations with M3a, but there are some variations in the relative energies of the stationary points, with a mean absolute deviation (MAD) of 2.7 kcal/mol. The energies of the two relevant transition states are close to each other at all levels applied (typically within 2 kcal/mol), with the first (second) one being rate-limiting in the QM/MM calculations with M3a (M4). QM-only gas-phase calculations give a very similar energy profile for QM region M4 (MAD of 1.7 kcal/mol), contrary to the situation for M3a where we had previously found significant discrepancies between the QM-only and QM/MM results (MAD of 7.9 kcal/mol). Extension of the QM region beyond M4 up to M7 (657 atoms) leads to only rather small variations in the relative energies from single-point QM-only and QM/MM calculations (MAD typically about 1–2 kcal/mol). In the case of acetylene hydratase, a model with 408 QM atoms thus seems sufficient to achieve convergence in the computed relative energies to within 1–2 kcal/mol.Copyright © 2013 Wiley Periodicals, Inc.

Convergence in the QM-only and QM/MM modeling of enzymatic reactions: A case study for acetylene hydratase.

Acetylene hydratase is a tungsten-dependent enzyme that catalyzes the nonredox hydration of acetylene to acetaldehyde. Density functional theory calculations are used to elucidate the reaction mechanism of this enzyme with a large model of the active site devised on the basis of the native X-ray crystal structure. Based on the calculations, we propose a new mechanism in which the acetylene substrate first displaces the W-coordinated water molecule, and then undergoes a nucleophilic attack by the water molecule assisted by an ionized Asp13 residue at the active site. This is followed by proton transfer from Asp13 to the newly formed vinyl anion intermediate. In the subsequent isomerization, Asp13 shuttles a proton from the hydroxyl group of the vinyl alcohol to the α-carbon. Asp13 is thus a key player in the mechanism, but also W is directly involved in the reaction by binding and activating acetylene and providing electrostatic stabilization to the transition states and intermediates. Several other mechanisms are also considered but the energetic barriers are found to be very high, ruling out these possibilities.

https://www.pnas.org/content/pnas/107/52/22523.full.pdf

Mechanism of tungsten-dependent acetylene hydratase from quantum chemical calculations

We report a comparison of QM-only and QM/MM approaches for the modeling of enzymatic reactions For this purpose, we present a QM/MM case study on the formation of vinyl alcohol in the catalytic cycle of tungsten-dependent acetylene hydratase Three different QM regions ranging from 32 to 157 atoms are designed for the reinvestigation of the previously suggested one-water attack mechanism The QM/MM calculations with the minimal QM region M1 (32 atoms) yield a two-step reaction profile, with an initial nucleophilic attack followed by the protonation of the formed vinyl anion intermediate, as previously proposed on the basis of QM-only calculations on cluster model M2 (116 atoms); however, the overall QM/MM barrier with M1 is much too high, mainly due to an overestimate of the QM/MM electrostatic repulsions QM/MM calculations with QM region M2 (116 atoms) fail to reproduce the published QM-only results, giving a one-step profile with a very high barrier This is traced back to the strong electrostatic influence of the two neighboring diphosphate groups that were neglected in the QM-only work but are present at the QM/MM level These diphosphate groups and other electrostatically important nearby residues are included in QM region M3 (157 atoms) QM/MM calculations with M3 recover the two-step mechanism and yield a reasonable overall barrier of 167 kcal/mol at the B3LYP/MM level They thus lead to a similar overall mechanistic scenario as the previous QM-only calculations, but there are also some important variations Most notably, the initial nucleophilic attack becomes rate limiting at the QM/MM level A modified two-water attack mechanism is also considered but is found to be less favorable than the previously proposed one-water attack mechanism Detailed residue interaction analyses and comparisons between QM/MM results with electronic and mechanical embedding and QM-only results without and with continuum solvation show that the protein environment plays a key role in determining the mechanistic preferences in acetylene hydratase The combined use of QM-only and QM/MM methods provides a powerful approach for the modeling of enzyme catalysis

/pdf/comparison-of-qm-only-and-qm-mm-models-for-the-mechanism-of-wbinvt266a.pdf

Rong-Zhen Liao

Papers

Quantum chemical studies of mechanisms for metalloenzymes.

Electrocatalytic water oxidation by a dinuclear copper complex in a neutral aqueous solution.

Convergence in the QM-only and QM/MM modeling of enzymatic reactions: A case study for acetylene hydratase.

Mechanism of tungsten-dependent acetylene hydratase from quantum chemical calculations

Comparison of QM-Only and QM/MM Models for the Mechanism of Tungsten-dependent Acetylene Hydratase