An Active Site Tyr Residue Guides the Regioselectivity of Lysine Hydroxylation by Nonheme Iron Lysine-4-hydroxylase Enzymes through Proton-Coupled Electron Transfer.

Lysine dioxygenase (KDO) is an important enzyme in human physiology involved in bioprocesses that trigger collagen cross-linking and blood pressure control. There are several KDOs in nature; however, little is known about the factors that govern the regio- and stereoselectivity of these enzymes. To understand how KDOs can selectively hydroxylate their substrate, we did a comprehensive computational study into the mechanisms and features of 4-lysine dioxygenase. In particular, we selected a snapshot from the MD simulation on KDO5 and created large QM cluster models (A, B, and C) containing 297, 312, and 407 atoms, respectively. The largest model predicts regioselectivity that matches experimental observation with rate-determining hydrogen atom abstraction from the C4-H position, followed by fast OH rebound to form 4-hydroxylysine products. The calculations show that in model C, the dipole moment is positioned along the C4-H bond of the substrate and, therefore, the electrostatic and electric field perturbations of the protein assist the enzyme in creating C4-H hydroxylation selectivity. Furthermore, an active site Tyr233 residue is identified that reacts through proton-coupled electron transfer akin to the axial Trp residue in cytochrome c peroxidase. Thus, upon formation of the iron(IV)-oxo species in the catalytic cycle, the Tyr233 phenol loses a proton to the nearby Asp179 residue, while at the same time, an electron is transferred to the iron to create an iron(III)-oxo active species. This charged tyrosyl residue directs the dipole moment along the C4-H bond of the substrate and guides the selectivity to the C4-hydroxylation of the substrate.

Axial Ligation Impedes Proton-Coupled Electron-Transfer Reactivity of a Synthetic Compound-I Analogue.

The nature of the axial ligand in high-valent iron-oxo heme enzyme intermediates and related synthetic catalysts is a critical structural element for controlling proton-coupled electron-transfer (PCET) reactivity of these species. Herein, we describe the generation and characterization of three new 6-coordinate, iron(IV)-oxo porphyrinoid-π-cation-radical complexes and report their PCET reactivity together with a previously published 5-coordinate analogue, FeIV(O)(TBP8Cz+•) (TBP8Cz = octakis(p-tert-butylphenyl)corrolazinato3-) (2) (Cho, K. A high-valent iron-oxo corrolazine activates C-H bonds via hydrogen-atom transfer. J. Am. Chem. Soc. 2012, 134, 7392-7399). The new complexes FeIV(O)(TBP8Cz+•)(L) (L = 1-methyl imidazole (1-MeIm) (4a), 4-dimethylaminopyridine (DMAP) (4b), cyanide (CN-)(4c)) can be generated from either oxidation of the ferric precursors or by addition of L to the Compound-I (Cpd-I) analogue at low temperatures. These complexes were characterized by UV-vis, electron paramagnetic resonance (EPR), and Mössbauer spectroscopies, and cryospray ionization mass spectrometry (CSI-MS). Kinetic studies using 4-OMe-TEMPOH as a test substrate indicate that coordination of a sixth axial ligand dramatically lowers the PCET reactivity of the Cpd-I analogue (rates up to 7000 times slower). Extensive density functional theory (DFT) calculations together with the experimental data show that the trend in reactivity with the axial ligands does not correlate with the thermodynamic driving force for these reactions or the calculated strengths of the O-H bonds being formed in the FeIV(O-H) products, pointing to non-Bell-Evans-Polanyi behavior. However, the PCET reactivity does follow a trend with the bracketed reduction potential of Cpd-I analogues and calculated electron affinities. The combined data suggest a concerted mechanism (a concerted proton electron transfer (CPET)) and an asynchronous movement of the electron/proton pair in the transition state.

QM/MM Study Into the Mechanism of Oxidative C=C Double Bond Cleavage by Lignostilbene-α,ß-Dioxygenase.

The enzymatic biosynthesis of fragrance molecules from lignin fragments is an important reaction in biotechnology for the sustainable production of fine chemicals. In this work we investigated the biosynthesis of vanillin from lignostilbene by a nonheme iron dioxygenase using QM/MM and tested several suggested proposals via either an epoxide or dioxetane intermediate. Binding of dioxygen to the active site of the protein results in the formation of an iron(II)-superoxo species with lignostilbene cation radical. The dioxygenase mechanism starts with electrophilic attack of the terminal oxygen atom of the superoxo group on the central C=C bond of lignostilbene, and the second-coordination sphere effects in the substrate binding pocket guide the reaction towards dioxetane formation. The computed mechanism is rationalized with thermochemical cycles and valence bond schemes that explain the electron transfer processes during the reaction mechanism. Particularly, the polarity of the protein and the local electric field and dipole moments enable a facile electron transfer and an exergonic dioxetane formation pathway.

Computational Study Into the Oxidative Ring-Closure Mechanism During the Biosynthesis of Deoxypodophyllotoxin.

The nonheme iron dioxygenase deoxypodophyllotoxin synthase performs an oxidative ring-closure reaction as part of natural product synthesis in plants. How the enzyme enables the oxidative ring-closure reaction of (-)-yatein and avoids substrate hydroxylation remains unknown. To gain insight into the reaction mechanism and understand the details of the pathways leading to products and by-products we performed a comprehensive computational study. The work shows that substrate is bound tightly into the substrate binding pocket with the C7'-H bond closest to the iron(IV)-oxo species. The reaction proceeds through a radical mechanism starting with hydrogen atom abstraction from the C7'-H position followed by ring-closure and a final hydrogen transfer to form iron(II)-water and deoxypodophyllotoxin. Alternative mechanisms including substrate hydroxylation and an electron transfer pathway were explored but found to be higher in energy. The mechanism is guided by electrostatic perturbations of charged residues in the second-coordination sphere that prevent alternative pathways.

What is the Origin of the Regioselective C3-Hydroxylation of L-Arg by the Nonheme Iron Enzyme Capreomycin C?

The nonheme iron dioxygenase capreomycin C (CmnC) hydroxylates a free L-arginine amino acid regio- and stereospecifically at the C3-position as part of the capreomycin antibiotics biosynthesis. Little is known on its structure, catalytic cycle and substrate specificity and, therefore, a comprehensive computational study was performed. A large QM cluster model of CmnC was created of 297 atoms and the mechanisms for C3-H, C4-H and C5-H hydroxylation and C3-C4 desaturation were investigated. All low-energy pathways correspond to radical reaction mechanisms with an initial hydrogen atom abstraction followed by OH rebound to form alcohol product complexes. The work is compared to alternative L-Arg hydroxylating nonheme iron dioxygenases and the differences in active site polarity are compared. We show that a tight hydrogen bonding network in the substrate binding pocket positions the substrate in an ideal orientation for C3-H activation, whereby the polar groups in the substrate binding pocket induce an electric field effect that guides the selectivity.

Tutorial Review on the Set-Up and Running of Quantum Mechanical Cluster Models for Enzymatic Reaction Mechanisms.

Enzymes turnover substrates into products with amazing efficiency and selectivity and as such have great potential for use in biotechnology and pharmaceutical applications. However, details of their catalytic cycles and the origins surrounding the regio- and chemoselectivity of enzymatic reaction processes remain unknown, which makes the engineering of enzymes and their use in biotechnology challenging. Computational modelling can assist experimental work in the field and establish the factors that influence the reaction rates and the product distributions. A popular approach in modelling is the use of quantum mechanical cluster models of enzymes that take the first- and second coordination sphere of the enzyme active site into consideration. These QM cluster models are widely applied but often the results obtained are dependent on model choice and model selection. Herein, we show that QM cluster models can give highly accurate results that reproduce experimental product distributions and free energies of activation within several kcal mol-1, regarded that large cluster models with >300 atoms are used that include key hydrogen bonding interactions and charged residues. In this tutorial review, we give general guidelines on the set-up and applications of the QM cluster method and discuss its accuracy and reproducibility. Finally, several representative QM cluster model examples on metal-containing enzymes are presented, which highlight the strength of the approach.

CO2 Reduction by an Iron(I) Porphyrinate System: Effect of Hydrogen Bonding on the Second Coordination Sphere.

Transforming CO2 into valuable materials is an important reaction in catalysis, especially because CO2 concentrations in the atmosphere have been growing steadily due to extensive fossil fuel usage. From an environmental perspective, reduction of CO2 to valuable materials should be catalyzed by an environmentally benign catalyst and avoid the use of heavy transition-metal ions. In this work, we present a computational study into a novel iron(I) porphyrin catalyst for CO2 reduction, namely, with a tetraphenylporphyrin ligand and analogues. In particular, we investigated iron(I) tetraphenylporphyrin with one of the meso-phenyl groups substituted with o-urea, p-urea, or o-2-amide groups. These substituents can provide hydrogen-bonding interactions in the second coordination sphere with bound ligands and assist with proton relay. Furthermore, our studies investigated bicarbonate and phenol as stabilizers and proton donors in the reaction mechanism. Potential energy landscapes for double protonation of iron(I) porphyrinate with bound CO2 are reported. The work shows that the bicarbonate bridges the urea/amide groups to the CO2 and iron center and provides a tight bonding pattern with strong hydrogen-bonding interactions that facilitates easy proton delivery and reduction of CO2. Specifically, bicarbonate provides a low-energy proton shuttle mechanism to form CO and water efficiently. Furthermore, the o-urea group locks bicarbonate and CO2 in a tight orientation and helps with ideal proton transfer, while there is more mobility and lesser stability with an o-amide group in that position instead. Our calculations show that the o-urea group leads to reduction in proton-transfer barriers, in line with experimental observation. We then applied electric-field-effect calculations to estimate the environmental effects on the two proton-transfer steps in the reaction. These calculations describe the perturbations that enhance the driving forces for the proton-transfer steps and have been used to make predictions about how the catalysts can be further engineered for more enhanced CO2 reduction processes.

Enhanced Reactivity through Equatorial Sulfur Coordination in Nonheme Iron(IV)-Oxo Complexes: Insights from Experiment and Theory.

Sulfur ligation in metalloenzymes often gives the active site unique properties, whether it is the axial cysteinate ligand in the cytochrome P450s or the equatorial sulfur/thiol ligation in nonheme iron enzymes. To understand sulfur ligation to iron complexes and how it affects the structural, spectroscopic, and intrinsic properties of the active species and the catalysis of substrates, we pursued a systematic study and compared sulfur with amine-ligated iron(IV)-oxo complexes. We synthesized and characterized a biomimetic N4S-ligated iron(IV)-oxo complex and compared the obtained results with an analogous N5-ligated iron(IV)-oxo complex. Our work shows that the amine for sulfur replacement in the equatorial ligand framework leads to a rate enhancement for oxygen atom and hydrogen atom transfer reactions. Moreover, the sulfur-ligated iron(IV)-oxo complex reacts through a different reaction mechanism as compared to the N5-ligated iron(IV)-oxo complex, where the former reacts through hydride transfer with the latter reacting via radical pathways. We show that the reactivity differences are caused by a dramatic change in redox potential between the two complexes. Our studies highlight the importance of implementing a sulfur ligand into the equatorial ligand framework of nonheme iron(IV)-oxo complexes and how it affects the physicochemical properties of the oxidant and its reactivity.

Machine learning-aided engineering of a cytochrome P450 for optimal bioconversion of lignin fragments.

Using machine learning, molecular dynamics simulations, and density functional theory calculations we gain insight into the selectivity patterns of substrate activation by the cytochromes P450. In nature, the reactions catalyzed by the P450s lead to the biodegradation of xenobiotics, but recent work has shown that fungi utilize P450s for the activation of lignin fragments, such as monomer and dimer units. These fragments often are the building blocks of valuable materials, including drug molecules and fragrances, hence a highly selective biocatalyst that can produce these compounds in good yield with high selectivity would be an important step in biotechnology. In this work a detailed computational study is reported on two reaction channels of two P450 isozymes, namely the O-deethylation of guaethol by CYP255A and the O-demethylation versus aromatic hydroxylation of p-anisic acid by CYP199A4. The studies show that the second-coordination sphere plays a major role in substrate binding and positioning, heme access, and in the selectivity patterns. Moreover, the local environment affects the kinetics of the reaction through lowering or raising barrier heights. Furthermore, we predict a site-selective mutation for highly specific reaction channels for CYP199A4.

Mechanism of the Oxidative Ring-Closure Reaction during Gliotoxin Biosynthesis by Cytochrome P450 GliF.

During gliotoxin biosynthesis in fungi, the cytochrome P450 GliF enzyme catalyzes an unusual C-N ring-closure step while also an aromatic ring is hydroxylated in the same reaction cycle, which may have relevance to drug synthesis reactions in biotechnology. However, as the details of the reaction mechanism are still controversial, no applications have been developed yet. To resolve the mechanism of gliotoxin biosynthesis and gain insight into the steps leading to ring-closure, we ran a combination of molecular dynamics and density functional theory calculations on the structure and reactivity of P450 GliF and tested a range of possible reaction mechanisms, pathways and models. The calculations show that, rather than hydrogen atom transfer from the substrate to Compound I, an initial proton transfer transition state is followed by a fast electron transfer en route to the radical intermediate, and hence a non-synchronous hydrogen atom abstraction takes place. The radical intermediate then reacts by OH rebound to the aromatic ring to form a biradical in the substrate that, through ring-closure between the radical centers, gives gliotoxin products. Interestingly, the structure and energetics of the reaction mechanisms appear little affected by the addition of polar groups to the model and hence we predict that the reaction can be catalyzed by other P450 isozymes that also bind the same substrate. Alternative pathways, such as a pathway starting with an electrophilic attack on the arene to form an epoxide, are high in energy and are ruled out.

A Cytochrome P450 TxtE Model System with Mechanistic and Theoretical Evidence for a Heme Peroxynitrite Active Species.

The cytochrome P450 homolog, TxtE, efficiently catalyzes the direct and regioselective aromatic nitration of the indolyl moiety of L-tryptophan to 4-nitro-L-tryptophan, using nitric oxide (NO) and dioxygen (O2) as co-substrates. Pathways for such direct and selective nitration of heteroaromatic motifs present platforms for engineering new nitration biocatalysts for pharmacologically beneficial targets, among a medley of other pivotal industrial applications. Precise mechanistic details concerning this pathway are only weakly understood, albeit a heme iron(III)-peroxynitrite active species has been postulated. To shed light on this unique reaction landscape, we investigated the indole nitration pathway of a series of biomimetic ferric heme superoxide mimics, [(Por)FeIII(O2 -⋅)], in the presence of NO. Therein, our model systems gave rise to three distinct nitroindole products, including 4-nitroindole, the product analogous to that obtained with TxtE. Moreover, 15N and 18O isotope labeling studies, along with meticulously designed control experiments lend credence to a heme peroxynitrite active nitrating agent, drawing close similarities to the tryptophan nitration mechanism of TxtE. All organic and inorganic reaction components have been fully characterized using spectroscopic methods. Theoretical investigation into several mechanistic possibilities deem a unique indolyl radical based reaction pathway as the most energetically favorable, products of which, are in excellent agreement with experimental findings.

Underlying Role of Hydrophobic Environments in Tuning Metal Elements for Efficient Enzyme Catalysis.

The catalytic functions of metalloenzymes are often strongly correlated with metal elements in the active sites. However, dioxygen-activating nonheme quercetin dioxygenases (QueD) are found with various first-row transition-metal ions when metal swapping inactivates their innate catalytic activity. To unveil the molecular basis of this seemingly promiscuous yet metal-specific enzyme, we transformed manganese-dependent QueD into a nickel-dependent enzyme by sequence- and structure-based directed evolution. Although the net effect of acquired mutations was primarily to rearrange hydrophobic residues in the active site pocket, biochemical, kinetic, X-ray crystallographic, spectroscopic, and computational studies suggest that these modifications in the secondary coordination spheres can adjust the electronic structure of the enzyme-substrate complex to counteract the effects induced by the metal substitution. These results explicitly demonstrate that such noncovalent interactions encrypt metal specificity in a finely modulated manner, revealing the underestimated chemical power of the hydrophobic sequence network in enzyme catalysis.

Local Charge Distributions, Electric Dipole Moments, and Local Electric Fields Influence Reactivity Patterns and Guide Regioselectivities in α-Ketoglutarate-Dependent Non-heme Iron Dioxygenases.

Non-heme iron dioxygenases catalyze vital processes for human health related to the biosynthesis of essential products and the biodegradation of toxic metabolites. Often the natural product biosyntheses by these non-heme iron dioxygenases is highly regio- and chemoselective, which are commonly assigned to tight substrate-binding and positioning. However, recent high-level computational modeling has shown that substrate-binding and positioning is only part of the story and long-range electrostatic interactions can play a major additional role.In this Account, we review and summarize computational viewpoints on the high regio- and chemoselectivity of α-ketoglutarate-dependent non-heme iron dioxygenases and how external perturbations affect the catalysis. In particular, studies from our groups have shown that often a regioselectivity in enzymes can be accomplished by stabilization of the rate-determining transition state for the reaction through external charges, electric dipole moments, or local electric field effects. Furthermore, bond dissociation energies in molecules are shown to be influenced by an electric field effect, and through targeting a specific bond in an electric field, this can lead to an unusually specific reaction. For instance, in the carbon-induced starvation protein, we studied two substrate-bound conformations and showed that regardless of what C-H bond of the substrate is closest to the iron(IV)-oxo oxidant, the lowest hydrogen atom abstraction barrier is always for the pro-S C2-H abstraction due to an induced dipole moment of the protein that weakens this bond. In another example of the hygromycin biosynthesis enzyme, an oxidative ring-closure reaction in the substrate forms an ortho-δ-ester ring. Calculations on this enzyme show that the selectivity is guided by a protonated lysine residue in the active site that, through its positive charge, triggers a low energy hydrogen atom abstraction barrier. A final set of examples in this Account discuss the viomycin biosynthesis enzyme and the 2-(trimethylammonio)ethylphosphonate dioxygenase (TmpA) enzyme. Both of these enzymes are shown to possess a significant local dipole moment and local electric field effect due to charged residues surrounding the substrate and oxidant binding pockets. The protein dipole moment and local electric field strength changes the C-H bond strengths of the substrate as compared to the gas-phase triggers the regioselectivity of substrate activation. In particular, we show that in the gas phase and in a protein environment C-H bond strengths are different due to local electric dipole moments and electric field strengths. These examples show that enzymes have an intricately designed structure that enables a chemical reaction under ambient conditions through the positioning of positively and negatively charged residues that influence and enhance reaction mechanisms. These computational insights create huge possibilities in bioengineering to apply local electric field and dipole moments in proteins to achieve an unusual selectivity and specificity and trigger a fit-for-purpose biocatalyst for unique biotransformations.

Caffeine Biodegradation by Cytochrome P450 1A2. What Determines the Product Distributions?

Caffeine is a natural compound found in plant seeds that after consumption by humans effects the central nervous system as well as the cardiovascular system. In general, the cytochrome P450 enzymes in the liver are involved in the biodegradation of caffeine, which gives paraxanthine, theobromine and theophylline products. There has been debate for many years why multiple products are obtained and how their distributions are determined. To this end we performed a high-level computational study using a combination of molecular dynamics and quantum mechanical approaches. A series of quantum chemical cluster models on the mechanism of caffeine activation by P450 model complexes give hydrogen atom abstraction barriers that predicts the correct ordering and statistical distribution of products. Our studies highlight that second-coordination sphere effects and thermochemical properties of the substrate determine the product distributions.

Mechanism of CO2 Reduction to Methanol with H2 on an Iron(II)-scorpionate Catalyst.

CO2 utilization is an important process in the chemical industry with great environmental power. In this work we show how CO2 and H2 can be reacted to form methanol on an iron(II) center and highlight the bottlenecks for the reaction and what structural features of the catalyst are essential for efficient turnover. The calculations predict the reactions to proceed through three successive reaction cycles that start with heterolytic cleavage of H2 followed by sequential hydride and proton transfer processes. The H2 splitting process is an endergonic process and hence high pressures will be needed to overcome this step and trigger the hydrogenation reaction. Moreover, H2 cleavage into a hydride and proton requires a metal to bind hydride and a nearby source to bind the proton, such as an amide or pyrazolyl group, which the scorpionate ligand used here facilitates. As such the computations highlight the non-innocence of the ligand scaffold through proton shuttle from H2 to substrate as an important step in the reaction mechanism.

Reactivity Differences of Trigonal Pyramidal Nonheme Iron(IV)-Oxo and Iron(III)-Oxo Complexes: Experiment and Theory.

High-valent metal-oxo species play critical roles in enzymatic catalysis yet their properties are still poorly understood. In this work we report a combined experimental and computational study into biomimetic iron(IV)-oxo and iron(III)-oxo complexes with tight second-coordination sphere environments that restrict substrate access. The work shows that the second-coordination sphere slows the hydrogen atom abstraction step from toluene dramatically and the kinetics is zeroth order in substrate. However, the iron(II)-hydroxo that is formed has a low reduction potential and hence cannot do OH rebound favorably. The tolyl radical in solution then reacts further with alternative reaction partners. By contrast, the iron(IV)-oxo species reacts predominantly through OH rebound to form alcohol products. Our studies show that the oxidation state of the metal influences reactivities and selectivities with substrate dramatically and that enzymes will likely need an iron(IV) center to catalyze C-H hydroxylation reactions.

Defluorination of Fluorophenols by a Nonheme Iron(IV)-Oxo Species: Observation of a New Intermediate Along the Reaction.

High-valent iron(IV)-oxo intermediates are versatile oxidants in the biotransformation of various substrates by metalloenzymes and catalyze essential reactions for human health as well as in the biodegradation of toxic organic pollutants in the environment. Herein, we report a biomimetic system that efficiently reacts with halophenols through defluorination reactions and characterize various short-lived intermediates along the reaction mechanism. We study the reactivity pattern of a nonheme iron(IV)-oxo species with a series of trihalophenols (X=F, Cl, Br). A combined experimental and computational study reveals that the oxidative dehalogenation of 2,4,6-trifluorophenol is initiated with an H-atom abstraction from the phenolic group by the iron(IV)-oxo species resulting in the formation of a phenolate radical and an iron(III)-hydroxo species. This iron(III)-hydroxo species forms an adduct with the oxidized substrate with λmax at 558 nm which subsequently decays to give quinones as products.

Computational Study on the Influence of Mo/V Centers on the Electronic Structure and Hydrazine Reduction Capability of [MFe3S4]3+/2+ Complexes.

[MFe3S4] cubanes have for some time been of interest for their ability to mimic the electronic and geometric structure of the active site of nitrogenase, the enzyme responsible for fixing N2 to NH3. Nitrogenase naturally occurs in three forms, with the major difference being that the metal ion present in the cofactor active site is either molybdenum (FeMoco), vanadium (FeVco), or iron. The molybdenum and vanadium versions of these cofactors are more closely studied, owing to their larger abundance and rate of catalysis. In this study, we compare free energy profiles and electronic properties of the Mo/V cubanes at various stages during the reduction of N2H4 to NH3. Our findings highlight the differences in how the complexes facilitate the reaction, in particular, vanadium's comparatively weaker ability to interact with the Fe/S network and stabilize reducing electrons prior to N-N bond cleavage, which may have implications when considering the lower efficiency of the vanadium-dependent nitrogenase.

Mechanism of Nitrogen Reduction to Ammonia in a Diiron Model of Nitrogenase.

Nitrogenase is a fascinating enzyme in biology that reduces dinitrogen from air to ammonia through stepwise reduction and protonation. Despite it being studied in detail by experimental and computational groups, there are still many unknown factors in the catalytic cycle of nitrogenase, especially related to the addition of protons and electrons and their order. A recent biomimetic study characterized a potential dinitrogen-bridged diiron cluster as a synthetic model of nitrogenase. Using strong acid and reductants, the dinitrogen was converted into ammonia molecules, but details of the mechanism remains unknown. In particular, it was unclear from the experimental studies whether the proton and electron transfer steps are sequential or alternating. Moreover, the work failed to establish what the function of the diiron core is and whether it split into mononuclear iron fragments during the reaction. To understand the structure and reactivity of the biomimetic dinitrogen-bridged diiron complex [(P2P'PhFeH)2(µ-N2)] with triphenylphosphine ligands, we performed a density functional theory study. Our computational methods were validated against experimental crystal structure coordinates, Mössbauer parameters, and vibrational frequencies and show excellent agreement. Subsequently, we investigated the alternating and consecutive addition of electrons and protons to the system. The calculations identify a number of possible reaction channels, namely, same-site protonation, alternating protonation, and complex dissociation into mononuclear iron centers. The calculations show that the overall mechanism is not a pure sequential set of electron and proton transfers but a mixture of alternating and consecutive steps. In particular, the first reaction steps will start with double proton transfer followed by an electron transfer, while thereafter, there is another proton transfer and a second electron transfer to give a complex whereby ammonia can split off with a low energetic barrier. The second channel starts with alternating protonation of the two nitrogen atoms, whereafter the initial double proton transfer, electrons and protons are added sequentially to form a hydrazine-bound complex. The latter split off ammonia spontaneously after further protonation. The various reaction channels are analyzed with valence bond and orbital diagrams. We anticipate the nitrogenase enzyme to operate with mixed alternating and consecutive protonation and electron transfer steps.

Biotransformation of Bisphenol by Human Cytochrome P450 2C9 Enzymes: A Density Functional Theory Study.

Bisphenol A (BPA, 2,2-bis-(4-hydroxyphenyl)propane) is used as a precursor in the synthesis of polycarbonate and epoxy plastics; however, its availability in the environment is causing toxicity as an endocrine-disrupting chemical. Metabolism of BPA and their analogues (substitutes) is generally performed by liver cytochrome P450 enzymes and often leads to a mixture of products, and some of those are toxic. To understand the product distributions of P450 activation of BPA, we have performed a computational study into the mechanisms and reactivities using large model structures of a human P450 isozyme (P450 2C9) with BPA bound. Density functional theory (DFT) calculations on mechanisms of BPA activation by a P450 compound I model were investigated, leading to a number of possible products. The substrate-binding pocket is tight, and as a consequence, aliphatic hydroxylation is not feasible as the methyl substituents of BPA cannot reach compound I well due to constraints of the substrate-binding pocket. Instead, we find low-energy pathways that are initiated with phenol hydrogen atom abstraction followed by OH rebound to the phenolic ortho- or para-position. The barriers of para-rebound are well lower in energy than those for ortho-rebound, and consequently, our P450 2C9 model predicts dominant hydroxycumyl alcohol products. The reactions proceed through two-state reactivity on competing doublet and quartet spin state surfaces. The calculations show fast and efficient substrate activation on a doublet spin state surface with a rate-determining electrophilic addition step, while the quartet spin state surface has multiple high-energy barriers that can also lead to various side products including C4-aromatic hydroxylation. This work shows that product formation is more feasible on the low spin state, while the physicochemical properties of the substrate govern barrier heights of the rate-determining step of the reaction. Finally, the importance of the second-coordination sphere is highlighted that determines the product distributions and guides the bifurcation pathways.