Iron porphyrin carbenes constitute a new frontier of species with considerable synthetic potential. Exquisitely engineered myoglobin and cytochrome P450 enzymes can generate these complexes and facilitate the transformations they mediate. The current work harnesses density functional theoretical methods to provide insight into the electronic structure, formation, and N-H insertion reactivity of an iron porphyrin carbene, [Fe(Por)(SCH3)(CHCO2Et)](-), a model of a complex believed to exist in an experimentally studied artificial metalloenzyme. The ground state electronic structure of the terminal form of this complex is an open-shell singlet, with two antiferromagnetically coupled electrons residing on the iron center and carbene ligand. As we shall reveal, the bonding properties of [Fe(Por)(SCH3)(CHCO2Et)](-) are remarkably analogous to those of ferric heme superoxide complexes. The carbene forms by dinitrogen loss from ethyl diazoacetate. This reaction occurs preferentially through an open-shell singlet transition state: iron donates electron density to weaken the C-N bond undergoing cleavage. Once formed, the iron porphyrin carbene accomplishes N-H insertion via nucleophilic attack. The resulting ylide then rearranges, using an internal carbonyl base, to form an enol that leads to the product. The findings rationalize experimentally observed reactivity trends reported in artificial metalloenzymes employing iron porphyrin carbenes. Furthermore, these results suggest a possible expansion of enzymatic substrate scope, to include aliphatic amines. Thus, this work, among the first several computational explorations of these species, contributes insights and predictions to the surging interest in iron porphyrin carbenes and their synthetic potential.

This Account outlines interplay of theory and experiment in the quest to identify the reactive-spin-state in chemical reactions that possess a few spin-dependent routes. Metalloenzymes and synthetic models have forged in recent decades an area of increasing appeal, in which oxometal species bring about functionalization of hydrocarbons under mild conditions and via intriguing mechanisms that provide a glimpse of Nature's designs to harness these reactions. Prominent among these are oxoiron(IV) complexes, which are potent H-abstractors. One of the key properties of oxoirons is the presence of close-lying spin-states, which can mediate H-abstractions. As such, these complexes form a fascinating chapter of spin-state chemistry, in which chemical reactivity involves spin-state interchange, so-called two-state reactivity (TSR) and multistate reactivity (MSR). TSR and MSR pose mechanistic challenges. How can one determine the structure of the reactive transition state (TS) and its spin state for these mechanisms? Calculations can do it for us, but the challenge is to find experimental probes. There are, however, no clear kinetic signatures for the reactive-spin-state in such reactions. This is the paucity that our group has been trying to fill for sometime. Hence, it is timely to demonstrate how theory joins experiment in realizing this quest. This Account uses a set of the H-abstraction reactions of 24 synthetic oxoiron(IV) complexes and 11 hydrocarbons, together undergoing H-abstraction reactions with TSR/MSR options, which provide experimentally determined kinetic isotope effect (KIE) data. For this set, we demonstrate that comparing KIE results with calculated tunneling-augmented KIE (KIE) data leads to a clear identification of the reactive spin-state during H-abstraction reactions. In addition, generating KIE data for a reaction of interest, and comparing these to KIE values, provides the mechanistic chemist with a powerful capability to identify the reactive-TS in terms of not only its spin state but also its geometry and ligand-sphere constitution. Since tunneling "cuts through" barriers, it serves as a chemical selectivity factor. Thus, we show that in a family of oxoirons reacting with one hydrocarbon, the tunneling efficiency increases as the ligands become better electron donors. This generates counterintuitive-reactivity patterns, like antielectrophilic reactivity, and induces spin-state reactivity reversals because of differing steric demands of the corresponding TS species, etc. Finally, for the same series, the Account reaches intuitive understanding of tunneling trends. It is shown that the increase of ligand's donicity results in electrostatic narrowing of the barrier, while the decrease of donicity and increase of bond-order asymmetry in the TS (inter alia due to Bell-Evans-Polanyi effects) broadens the barrier. Predictions are made that usage of powerful electron-donating ligands may train H-abstractors to activate the strongest C-H bond in a molecule. The concepts developed here are likely to be ap...

We present here a versatile computational code named “elecTric fIeld generaTion And maNipulation (TITAN),” capable of generating various types of external electric fields, as well as quantifying the local (or intrinsic) electric fields present in proteins and other biological systems according to Coulomb's Law. The generated electric fields can be coupled with quantum mechanics (QM), molecular mechanics (MM), QM/MM, and molecular dynamics calculations in most available software packages. The capabilities of the TITAN code are illustrated throughout the text with the help of examples. We end by presenting an application, in which the effects of the local electric field on the hydrogen transfer reaction in cytochrome P450 OleTJE enzyme and the modifications induced by the application of an oriented external electric field are examined. We find that the protein matrix in P450 OleTJE acts as a moderate catalyst and that orienting an external electric field along the Fe─O bond of compound I has the biggest impact on the reaction barrier. The induced catalysis/inhibition correlates with the calculated spin density on the O‐atom. © 2019 Wiley Periodicals, Inc.

This contribution follows the recent remarkable catalysis observed by Groves et al. in hydrogen‐abstraction reactions by a) an oxoferryl porphyrin radical‐cation complex [Por⋅+FeIV(O)Lax] and b) a hydroxoiron porphyrazine ferric complex [PyPzFeIII(OH)Lax], both of which involve positively charged substituents on the outer circumference of the respective macrocyclic ligands. These charge‐coronated complexes are analogues of the biologically important Compound I (Cpd I) and synthetic hydroxoferric species, respectively. We demonstrate that the observed enhancement of the H‐abstraction catalysis for these systems is a purely electrostatic effect, elicited by the local charges embedded on the peripheries of the respective macrocyclic ligands. Our findings provide new insights into how electrostatics can be employed to tune the catalytic activity of metalloenzymes and can thus contribute to the future design of new and highly efficient hydrogen‐abstraction catalysts.

An H/D kinetic isotope effect (KIE) of 80 is found at -20 °C for the oxidation of 9,10-dihydroanthracene by [Fe(O)(TMCS)], a complex supported by the tetramethylcyclam (TMC) macrocycle with a tethered thiolate. This KIE value approaches that previously predicted by DFT calculations. Other [Fe(O)(TMC)(anion)] complexes exhibit values of 20, suggesting that the thiolate ligand of [Fe(O)(TMCS)] plays a unique role in facilitating tunneling. Calculations show that tunneling is most enhanced (a) when the bond asymmetry between C-H bond breaking and O-H bond formation in the transition state is minimized, and (b) when the electrostatic interactions in the O---H---C moiety are maximal. These two factors-which peak for the best electron donor, the thiolate ligand-afford a slim and narrow barrier through which the H-atom can tunnel most effectively.

Iron(III) complexes [FeL(B)] (1-5) of a tetradentate trianionic phenolate-based ligand (L) and modified dipyridophenazine bases (B), namely, dipyrido-6,7,8,9-tetrahydrophenazine (dpqC in 1), dipyrido[3,2-a:2',3'-c]phenazine-2-carboxylic acid (dppzc in 2), dipyrido[3,2-a:2',3'-c]phenazine-11-sulfonic acid (dppzs in 3), 7-aminodipyrido[3,2-a:2',3'-c]phenazine (dppza in 4) and benzo[i]dipyrido[3,2-a:2',3'-c]phenazine (dppn in 5), have been synthesized and their photocytotoxic properties studied along with their dipyridophenazine analogue (6). The complexes have a five electron paramagnetic iron(III) center, and the Fe(III)/Fe(II) redox couple appears at about -0.69 V versus SCE in DMF-0.1 M TBAP. The physicochemical data also suggest that the complexes possess similar structural features as that of its parent complex [FeL(dppz)] with FeO3N3 coordination in a distorted octahedral geometry. The DNA-complex and protein-complex interaction studies have revealed that the complexes interact favorably with the biomolecules, the degree of which depends on the nature of the substituents present on the dipyridophenazine ring. Photocleavage of pUC19 DNA by the complexes has been studied using visible light of 476, 530, and 647 nm wavelengths. Mechanistic investigations with inhibitors show formation of HO(•) radicals via a photoredox pathway. Photocytotoxicity study of the complexes in HeLa cells has shown that the dppn complex (5) is highly active in causing cell death in visible light with sub micromolar IC(50) value. The effect of substitutions and the planarity of the phenazine moiety on the cellular uptake are quantified by determining the total cellular iron content using the inductively coupled plasma-optical emission spectrometry (ICP-OES) technique. The cellular uptake increases marginally with an increase in the hydrophobicity of the dipyridophenazine ligands whereas complex 3 with dppzs shows very high uptake. Insights into the cell death mechanism by the dppn complex 5, obtained through DAPI nuclear staining in HeLa cells, reveal a rapid programmed cell death mechanism following photoactivation of complex 5 with visible light. The effect of substituent on the DNA photocleavage activity of the complexes has been rationalized from the theoretical studies.

Bispropargyl sulfones equipped with aromatic rings of dissimilar nature were synthesized. Under basic conditions, these sulfones isomerized to the bisallenic sulfones, creating a competitive scenario between two alternate Garratt-Braverman (GB) cyclization pathways. The observed product distribution ruled out the involvement of any ionic intermediate and supported the diradical mechanism with greater involvement of the electron-rich aromatic ring via the more nucleophilic radical. DFT-based calculations supported the diradical mechanism along with the observed selectivity.

A competitive scenario between Myers-Saito (MS) and Garratt-Braverman (GB) cyclization has been created in a molecule. High-level computations indicate a preference for GB over MS cyclization. The activation energies for the rate-determining steps of the GB and MS cyclizations were found to be the same (24.4 kcal/mol) at the B3LYP/6-31G* level of theory; thus, from the kinetic point of view, both reactions are feasible. However, the main biradical intermediate GB2 of the GB reaction is 6.2 kcal/mol lower in energy than the biradical MS2, which is the main intermediate of MS reaction, so GB cyclization is thermodynamically favored over MS cyclization. To verify the prediction by computational techniques, bisenediynyl sulfones 1-4 and bisenediynyl sulfoxide 17 were synthesized. Under basic conditions, these molecules isomerized to a system possessing both the ene-yne-allene and the bisallenic sulfone. The isolation of only one product, identified as the corresponding naphthalene- or benzene-fused sulfone 8-11, indicated the occurrence of GB cyclization as the sole reaction pathway. No product corresponding to the MS cyclization pathway could be isolated. Though the theoretical prediction showed a preference for the GB pathway over the MS pathway, the exclusive preference for GB over MS cyclization is very striking. Further analysis showed that the intramolecular self-quenching nature of the GB pathway may play an important role in the complete preference for this reaction. Apart from the mechanistic studies, these sulfones showed DNA cleavage activity that had an inverse relation with the reactivity order. Our findings are important for the design of artificial DNA-cleaving agents.