Matthew G. Quesne

Quantum Mechanics/Molecular Mechanics Study on the Oxygen Binding and Substrate Hydroxylation Step in AlkB Repair Enzymes

AlkB repair enzymes are important nonheme iron enzymes that catalyse the demethylation of alkylated DNA bases in humans, which is a vital reaction in the body that heals externally damaged DNA bases. Its mechanism is currently controversial and in order to resolve the catalytic mechanism of these enzymes, a quantum mechanics/molecular mechanics (QM/MM) study was performed on the demethylation of the N1-methyladenine fragment by AlkB repair enzymes. Firstly, the initial modelling identified the oxygen binding site of the enzyme. Secondly, the oxygen activation mechanism was investigated and a novel pathway was found, whereby the catalytically active iron(IV)–oxo intermediate in the catalytic cycle undergoes an initial isomerisation assisted by an Arg residue in the substrate binding pocket, which then brings the oxo group in close contact with the methyl group of the alkylated DNA base. This enables a subsequent rate-determining hydrogen-atom abstraction on competitive σ-and π-pathways on a quintet spin-state surface. These findings give evidence of different locations of the oxygen and substrate binding channels in the enzyme and the origin of the separation of the oxygen-bound intermediates in the catalytic cycle from substrate. Our studies are compared with small model complexes and the effect of protein and environment on the kinetics and mechanism is explained.

Drug Metabolism by Cytochrome P450 Enzymes: What Distinguishes the Pathways Leading to Substrate Hydroxylation Over Desaturation?

Cytochrome P450 enzymes are highly versatile biological catalysts in our body that react with a broad range of substrates. Key functions in the liver include the metabolism of drugs and xenobiotics. One particular metabolic pathway that is poorly understood relates to the P450 activation of aliphatic groups leading to either hydroxylation or desaturation pathways. A DFT and QM/MM study has been carried out on the factors that determine the regioselectivity of aliphatic hydroxylation over desaturation of compounds by P450 isozymes. The calculations establish multistate reactivity patterns, whereby the product distributions differ on each of the spin-state surfaces; hence spin-selective product formation was found. The electronic and thermochemical factors that determine the bifurcation pathways were analysed and a model that predicts the regioselectivity of aliphatic hydroxylation over desaturation pathways was established from valence bond and molecular orbital theories. Thus, the difference in energy of the OH versus the OC bond formed and the π-conjugation energy determines the degree of desaturation products. In addition, environmental effects of the substrate binding pocket that affect the regioselectivities were identified. These studies imply that bioengineering P450 isozymes for desaturation reactions will have to include modifications in the substrate binding pocket to restrict the hydroxylation rebound reaction.

Secondary coordination sphere influence on the reactivity of nonheme iron(II) complexes: an experimental and DFT approach.

The new biomimetic ligands N4Py2Ph (1) and N4Py2Ph,amide (2) were synthesized and yield the iron(II) complexes [FeII(N4Py2Ph)(NCCH3)](BF4)2 (3) and [FeII(N4Py2Ph,amide)](BF4)2 (5). Controlled orientation of the Ph substituents in 3 leads to facile triplet spin reactivity for a putative FeIV(O) intermediate, resulting in rapid arene hydroxylation. Addition of a peripheral amide substituent within hydrogen-bond distance of the iron first coordination sphere leads to stabilization of a high-spin FeIIIOOR species which decays without arene hydroxylation. These results provide new insights regarding the impact of secondary coordination sphere effects at nonheme iron centers.

Quantum Mechanics/Molecular Mechanics Modeling of Enzymatic Processes: Caveats and Breakthroughs

Nature has developed large groups of enzymatic catalysts with the aim to transfer substrates into useful products, which enables biosystems to perform all their natural functions. As such, all biochemical processes in our body (we drink, we eat, we breath, we sleep, etc.) are governed by enzymes. One of the problems associated with research on biocatalysts is that they react so fast that details of their reaction mechanisms cannot be obtained with experimental work. In recent years, major advances in computational hardware and software have been made and now large (bio)chemical systems can be studied using accurate computational techniques. One such technique is the quantum mechanics/molecular mechanics (QM/MM) technique, which has gained major momentum in recent years. Unfortunately, it is not a black-box method that is easily applied, but requires careful set-up procedures. In this work we give an overview on the technical difficulties and caveats of QM/MM and discuss work-protocols developed in our groups for running successful QM/MM calculations.

Origin of the Regioselective Fatty-Acid Hydroxylation versus Decarboxylation by a Cytochrome P450 Peroxygenase: What Drives the Reaction to Biofuel Production?

The cytochromes P450 are heme-based mono-oxygenases or peroxygenases involved in vital reaction processes for human health. A recently described P450 per-oxygenase, OleTJE , converts long-chain fatty acids to terminal olefins and as such may have biotechnological relevance in biodiesel production. However, the reaction produces significant amounts of α- and β-hydroxylation by-products, and their origin are poorly understood. Herein, we elucidate through a QM/MM study on the bifurcation pathways how the three possible products are generated and show how the enzyme can be further engineered for optimum desaturase activity. The studies showed that the polarity and the solvent accessibility of the substrate in the binding pocket destabilize the OH-rebound pathways and kinetically enable a thermodynamically otherwise unfavorable decarboxylation reaction. The origins of the bifurcation pathways are analyzed with valence-bond models that highlight the differences in reaction mechanism.

Oxygen-atom transfer reactivity of axially ligated Mn(V)-oxo complexes: Evidence for enhanced electrophilic and nucleophilic pathways

Addition of anionic donors to the manganese(V)–oxo corrolazine complex MnV(O)(TBP8Cz) has a dramatic influence on oxygen-atom transfer (OAT) reactivity with thioether substrates. The six-coordinate anionic [MnV(O)(TBP8Cz)(X)]− complexes (X = F–, N3–, OCN–) exhibit a ∼5 cm–1 downshift of the Mn–O vibrational mode relative to the parent MnV(O)(TBP8Cz) complex as seen by resonance Raman spectroscopy. Product analysis shows that the oxidation of thioether substrates gives sulfoxide product, consistent with single OAT. A wide range of OAT reactivity is seen for the different axial ligands, with the following trend determined from a comparison of their second-order rate constants for sulfoxidation: five-coordinate ≈ thiocyanate ≈ nitrate < cyanate < azide < fluoride ≪ cyanide. This trend correlates with DFT calculations on the binding of the axial donors to the parent MnV(O)(TBP8Cz) complex. A Hammett study was performed with p-X-C6H4SCH3 derivatives and [MnV(O)(TBP8Cz)(X)]− (X = CN– or F–) as the oxidant, and unusual “V-shaped” Hammett plots were obtained. These results are rationalized based upon a change in mechanism that hinges on the ability of the [MnV(O)(TBP8Cz)(X)]− complexes to function as either an electrophilic or weak nucleophilic oxidant depending upon the nature of the para-X substituents. For comparison, the one-electron-oxidized cationic MnV(O)(TBP8Cz•+) complex yielded a linear Hammett relationship for all substrates (ρ = −1.40), consistent with a straightforward electrophilic mechanism. This study provides new, fundamental insights regarding the influence of axial donors on high-valent MnV(O) porphyrinoid complexes.

Catalytic Mechanism of Cofactor-Free Dioxygenases and How They Circumvent Spin-Forbidden Oxygenation of Their Substrates

Dioxygenases catalyze a diverse range of biological reactions by incorporating molecular oxygen into organic substrates. Typically, they use transition metals or organic cofactors for catalysis. Bacterial 1-H-3-hydroxy-4-oxoquinaldine-2,4-dioxygenase (HOD) catalyzes the spin-forbidden transfer of dioxygen to its N-heteroaromatic substrate in the absence of any cofactor. We combined kinetics, spectroscopic and computational approaches to establish a novel reaction mechanism. The present work gives insight into the rate limiting steps in the reaction mechanism, the effect of first-coordination sphere amino acids as well as electron-donating/electron-withdrawing substituents on the substrate. We highlight the role of active site residues Ser101/Trp160/His251 and their involvement in the reaction mechanism. The work shows, for the first time, that the reaction is initiated by triplet dioxygen and its binding to deprotonated substrate and only thereafter a spin state crossing to the singlet spin state occurs. As revealed by steady- and transient-state kinetics the oxygen-dependent steps are rate-limiting, whereas Trp160 and His251 are essential residues for catalysis and contribute to substrate positioning and activation, respectively. Computational modeling further confirms the experimental observations and rationalizes the electron transfer pathways, and the effect of substrate and substrate binding pocket residues. Finally, we make a direct comparison with iron-based dioxygenases and explain the mechanistic and electronic differences with cofactor-free dioxygenases. Our multidisciplinary study confirms that the oxygenation reaction can take place in absence of any cofactor by a unique mechanism in which the specially designed fit-for-purpose active-site architecture modulates substrate reactivity toward oxygen.

Differences and Comparisons of the Properties and Reactivities of Iron(III)–hydroperoxo Complexes with Saturated Coordination Sphere

Heme and nonheme monoxygenases and dioxygenases catalyze important oxygen atom transfer reactions to substrates in the body. It is now well established that the cytochrome P450 enzymes react through the formation of a high-valent iron(IV)–oxo heme cation radical. Its precursor in the catalytic cycle, the iron(III)–hydroperoxo complex, was tested for catalytic activity and found to be a sluggish oxidant of hydroxylation, epoxidation and sulfoxidation reactions. In a recent twist of events, evidence has emerged of several nonheme iron(III)–hydroperoxo complexes that appear to react with substrates via oxygen atom transfer processes. Although it was not clear from these studies whether the iron(III)–hydroperoxo reacted directly with substrates or that an initial O–O bond cleavage preceded the reaction. Clearly, the catalytic activity of heme and nonheme iron(III)–hydroperoxo complexes is substantially different, but the origins of this are still poorly understood and warrant a detailed analysis. In this work, an extensive computational analysis of aromatic hydroxylation by biomimetic nonheme and heme iron systems is presented, starting from an iron(III)–hydroperoxo complex with pentadentate ligand system (L52). Direct C–O bond formation by an iron(III)–hydroperoxo complex is investigated, as well as the initial heterolytic and homolytic bond cleavage of the hydroperoxo group. The calculations show that [(L52)FeIII(OOH)]2+ should be able to initiate an aromatic hydroxylation process, although a low-energy homolytic cleavage pathway is only slightly higher in energy. A detailed valence bond and thermochemical analysis rationalizes the differences in chemical reactivity of heme and nonheme iron(III)–hydroperoxo and show that the main reason for this particular nonheme complex to be reactive comes from the fact that they homolytically split the O–O bond, whereas a heterolytic O–O bond breaking in heme iron(III)–hydroperoxo is found.

Direct Observation of a Nonheme Iron(IV)–Oxo Complex That Mediates Aromatic C–F Hydroxylation

The synthesis of a pentadentate ligand with strategically designed fluorinated arene groups in the second coordination sphere of a nonheme iron center is reported. The oxidatively resistant fluorine substituents allow for the trapping and characterization of an FeIV(O) complex at −20 °C. Upon warming of the FeIV(O) complex, an unprecedented arene C–F hydroxylation reaction occurs. Computational studies support the finding that substrate orientation is a critical factor in the observed reactivity. This work not only gives rare direct evidence for the participation of an FeIV(O) species in arene hydroxylation but also provides the first example of a high-valent iron–oxo complex that mediates aromatic C–F hydroxylation.

Identification and Spectroscopic Characterization of Nonheme Iron (III) Hypochlorite Intermediates

FeIII–hypohalite complexes have been implicated in a wide range of important enzyme-catalyzed halogenation reactions including the biosynthesis of natural products and antibiotics and post-translational modification of proteins. The absence of spectroscopic data on such species precludes their identification. Herein, we report the generation and spectroscopic characterization of nonheme FeIII–hypohalite intermediates of possible relevance to iron halogenases. We show that FeIII-OCl polypyridylamine complexes can be sufficiently stable at room temperature to be characterized by UV/Vis absorption, resonance Raman and EPR spectroscopies, and cryo-ESIMS. DFT methods rationalize the pathways to the formation of the FeIII-OCl, and ultimately FeIV=O, species and provide indirect evidence for a short-lived FeII-OCl intermediate. The species observed and the pathways involved offer insight into and, importantly, a spectroscopic database for the investigation of iron halogenases.