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

Theoretical Insight into the Fluorescence Spectral Tuning Mechanism: A Case Study of Flavin-Dependent Bacterial Luciferase

Bioluminescence of bacteria is widely applied in biological imaging, environmental toxicant detection, and many other situations. Understanding the spectral tuning mechanism not only helps explain the diversity of colors observed in nature but also provides principles for bioengineering new color variants for practical applications. In this study, time-dependent density functional theory (TD-DFT) and quantum mechanics and molecular mechanics (QM/MM) calculations have been employed to understand the fluorescence spectral tuning mechanism of bacterial luciferase with a focus on the electrostatic effect. The spectrum can be tuned by both a homogeneous dielectric environment and oriented external electric fields (OEEFs). Increasing the solvent polarity leads to a redshift of the fluorescence emission maximum, λF, accompanied by a substantial increase in density. In contrast, applying an OEEF along the long axis of the isoalloxazine ring (X-axis) leads to a significant red- or blue-shift in λF, depending on the direction of the OEEF, yet with much smaller changes in intensity. The effect of polar solvents is directionless, and the red-shifts can be attributed to the larger dipole moment of the S1 state compared with that of the S0 state. However, the effect of OEEFs directly correlates with the difference dipole moment between the S1 and S0 states, which is directional and is determined by the charge redistribution upon deexcitation. Moreover, the electrostatic effect of bacterial luciferase is in line with the presence of an internal electric field (IEF) pointing in the negative X direction. Finally, the key residues that contribute to this IEF and strategies for modulating the spectrum through site-directed point mutations are discussed.

Probing Ferryl Reactivity in a Nonheme Iron Oxygenase Using an Expanded Genetic Code

The ability to introduce noncanonical amino acids as axial ligands in heme enzymes has provided a powerful experimental tool for studying the structure and reactivity of their FeIV=O ("ferryl") intermediates. Here, we show that a similar approach can be used to perturb the conserved Fe coordination environment of 2-oxoglutarate (2OG) dependent oxygenases, a versatile class of enzymes that employ highly-reactive ferryl intermediates to mediate challenging C-H functionalizations. Replacement of one of the cis-disposed histidine ligands in the oxygenase VioC with a less electron donating N δ-methyl-histidine (MeHis) preserves both catalytic function and reaction selectivity. Significantly, the key ferryl intermediate responsible for C-H activation can be accumulated in both the wildtype and the modified protein. In contrast to heme enzymes, where metal-oxo reactivity is extremely sensitive to the nature of the proximal ligand, the rates of C-H activation and the observed large kinetic isotope effects are only minimally affected by axial ligand replacement in VioC. This study showcases a powerful tool for modulating the coordination sphere of nonheme iron enzymes that will enhance our understanding of the factors governing their divergent activities.

Unraveling the Mechanism of the Oxidative C-C Bond Coupling Reaction Catalyzed by Deoxypodophyllotoxin Synthase

Deoxypodophyllotoxin synthase (DPS), a nonheme Fe(II)/2-oxoglutarate (2OG)-dependent oxygenase, is a key enzyme that is involved in the construction of the fused-ring system in (-)-podophyllotoxin biosynthesis by catalyzing the C-C coupling reaction. However, the mechanistic details of DPS-catalyzed ring formation remain unclear. Herein, our quantum mechanics/molecular mechanics (QM/MM) calculations reveal a novel mechanism that involves the recycling of CO2 (a product of decarboxylation of 2OG) to prevent the formation of hydroxylated byproducts. Our results show that CO2 can react with the FeIII-OH species to generate an unusual FeIII-bicarbonate species. In this way, hydroxylation is avoided by consuming the OH group. Then, the C-C coupling followed by desaturation yields the final product, deoxypodophyllotoxin. This work highlights the crucial role of the CO2 molecule, generated in the crevice between the iron active site and the substrate, in controlling the reaction selectivity.

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.

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.

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.

Selective cascade activation of polycyclic aromatic hydrocarbons in human cells: Role of enzyme's intrinsic electric field

Polycyclic aromatic hydrocarbons (PAHs) are major environmental organic pollutants. Some metabolites of PAHs show greater toxicity to humans while the others do not. It is highly important to decipher PAHs' regioselective activation mechanism and identify the major metabolites to accurately evaluate their public health risk. Here, we have performed a thorough computational study of benzo[a]anthracene (BA) metabolized by P450 1A1 by employing molecular docking, molecular dynamics simulations, quantum chemical calculation, and quantum mechanics/molecular mechanics calculations. Our findings show that highly-reactive species such as 3,4-epoxide, 8,9-epoxide, 3,4-diol-1,2-epoxide, and 8,9-diol-10,11-epoxide were major metabolites, which can efficiently react with guanine and damage DNA with extremely low energy barrier, therefore, supports the regioselective metabolism of BA. The origin of this selective activation is mainly contributed to both the oxygen‑carbon distance and previously overlooked enzyme's intrinsic electric field. Consequently, based on the resolved activation selectivity of BA. We built a high-throughput strategy to efficiently predict the metabolites of other PAHs. The accuracy of the strategy is validated by studying 16 PAHs on the priority control list. Hopefully this will aid the accurate evaluation of public health risks associated with PAH emissions.

Catalytic divergencies in the mechanism of L-arginine hydroxylating nonheme iron enzymes

Many enzymes in nature utilize a free arginine (L-Arg) amino acid to initiate the biosynthesis of natural products. Examples include nitric oxide synthases, which generate NO from L-Arg for blood pressure control, and various arginine hydroxylases involved in antibiotic biosynthesis. Among the groups of arginine hydroxylases, several enzymes utilize a nonheme iron(II) active site and let L-Arg react with dioxygen and α-ketoglutarate to perform either C3-hydroxylation, C4-hydroxylation, C5-hydroxylation, or C4-C5-desaturation. How these seemingly similar enzymes can react with high specificity and selectivity to form different products remains unknown. Over the past few years, our groups have investigated the mechanisms of L-Arg-activating nonheme iron dioxygenases, including the viomycin biosynthesis enzyme VioC, the naphthyridinomycin biosynthesis enzyme NapI, and the streptothricin biosynthesis enzyme OrfP, using computational approaches and applied molecular dynamics, quantum mechanics on cluster models, and quantum mechanics/molecular mechanics (QM/MM) approaches. These studies not only highlight the differences in substrate and oxidant binding and positioning but also emphasize on electronic and electrostatic differences in the substrate-binding pockets of the enzymes. In particular, due to charge differences in the active site structures, there are changes in the local electric field and electric dipole moment orientations that either strengthen or weaken specific substrate C-H bonds. The local field effects, therefore, influence and guide reaction selectivity and specificity and give the enzymes their unique reactivity patterns. Computational work using either QM/MM or density functional theory (DFT) on cluster models can provide valuable insights into catalytic reaction mechanisms and produce accurate and reliable data that can be used to engineer proteins and synthetic catalysts to perform novel reaction pathways.

Mechanistic Insights into Amphoteric Reactivity of an Iron-Bispidine Complex

The reactivity of FeIII -alkylperoxido complexes has remained a riddle to inorganic chemists owing to their thermal instability and impotency towards organic substrates. These iron-oxygen adducts have been known as sluggish oxidants towards oxidative electrophilic and nucleophilic reactions. Herein, we report the synthesis and spectroscopic characterization of a relatively stable mononuclear high-spin FeIII -alkylperoxido complex supported by an engineered bispidine framework. Against the notion, this FeIII -alkylperoxido complex serves as a rare example of versatile reactivity in both electrophilic and nucleophilic reactions. Detailed mechanistic studies and computational calculations reveal a novel reaction mechanism, where a putative superoxido intermediate orchestrates the amphoteric property of the oxidant. The design of the backbone is pivotal to convey stability and reactivity to alkylperoxido and superoxido intermediates. Contrary to the well-known O-O bond cleavage that generates an FeIV -oxido species, the FeIII -alkylperoxido complex reported here undergoes O-C bond scission to generate a superoxido moiety that is responsible for the amphiphilic reactivity.

Monovalent metal ion binding promotes the first transesterification reaction in the spliceosome

Cleavage and formation of phosphodiester bonds in nucleic acids is accomplished by large cellular machineries composed of both protein and RNA. Long thought to rely on a two-metal-ion mechanism for catalysis, structure comparisons revealed many contain highly spatially conserved second-shell monovalent cations, whose precise function remains elusive. A recent high-resolution structure of the spliceosome, essential for pre-mRNA splicing in eukaryotes, revealed a potassium ion in the active site. Here, we employ biased quantum mechanics/ molecular mechanics molecular dynamics to elucidate the function of this monovalent ion in splicing. We discover that the K+ ion regulates the kinetics and thermodynamics of the first splicing step by rigidifying the active site and stabilizing the substrate in the pre- and post-catalytic state via formation of key hydrogen bonds. Our work supports a direct role for the K+ ion during catalysis and provides a mechanistic hypothesis likely shared by other nucleic acid processing enzymes.

The elusive reaction mechanism of Mn(II)-mediated benzylic oxidation of alkylarene by H2O2: a gem-diol mechanism or a dual hydrogen abstraction mechanism?

The direct oxygenation of alkylarenes at the benzylic position employing bioinspired nonheme catalysts has emerged as a promising strategy for the production of bioactive arene ketone scaffolds in drugs. However, the structure-activity relationship of the active species and the mechanism of these reactions remain elusive. Herein, the reaction mechanism of the Mn(II)-mediated benzylic oxygenation of phenylbutanoic acid (PBA) to 4-oxo-4-phenylbutyric acid (4-oxo-PBA) by H2O2 was investigated using density functional theory calculations. The calculated results demonstrated that the MnIII-OOH species (1) is a sluggish oxidant and needs to be converted to a high-valent manganese-oxo species (2). The conversion of PBA to 4-oxo-PBA by 2 occurs via the consecutive hydroxylation of PBA to 4-hydroxyl-4-phenylbutyric acid (4-OH-PBA) and the alcohol oxidation of 4-OH-PBA to 4-oxo-PBA. The hydroxylation of PBA proceeds via a novel hydride transfer/hydroxyl-rebound mechanism and the alcohol oxidation of 4-OH-PBA occurs via three pathways (gem-diol, dual hydrogen abstraction (DHA), and reversed-DHA pathways). The regio-selectivity of benzylic oxidations was caused by a strong π-π stacking interaction between the pyridine ring of the nonheme ligand and the phenyl ring of the substrate. These mechanistic findings enrich the knowledge of biomimetic alcohol oxidations and play a positive role in the rational design of new non-heme catalysts.

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.

Insights into Cytochrome P450 Enzyme Catalyzed Defluorination of Aromatic Fluorides

Density functional calculations establish a novel mechanism of aromatic defluorination by P450 Compound I. This is achieved via either an initial epoxide intermediate or through a 1,2-fluorine shift in an electrophilic intermediate, which highlights that the P450s can defluorinate fluoroarenes. However, in the absence of a proton donor a strong Fe-F bond can be obtained as shown from the calculations.

Preorganized Internal Electric Field Promotes a Double-Displacement Mechanism for the Adenine Excision Reaction by Adenine DNA Glycosylase

Adenine DNA glycosylase (MutY) is a monofunctional glycosylase, removing adenines (A) misinserted opposite 8-oxo-7,8-dihydroguanine (OG), a common product of oxidative damage to DNA. Through multiscale calculations, we decipher a detailed adenine excision mechanism of MutY that is consistent with all available experimental data, involving an initial protonation step and two nucleophilic displacement steps. During the first displacement step, N-glycosidic bond cleavage is accompanied by the attack of the carboxylate group of residue Asp144 at the anomeric carbon (C1'), forming a covalent glycosyl-enzyme intermediate to stabilize the fleeting oxocarbenium ion. After departure of the excised base, water nucleophiles can be recruited to displace Asp144, completing the catalytic cycle with retention of stereochemistry at the C1' position. The two displacement reactions are found to mostly involve the movement of the oxocarbenium ion, occurring with large charge reorganization and thus sensitive to the internal electric field (IEF) exerted by the polar protein environment. Intriguingly, we find that the negatively charged carboxylate group is a good nucleophile for the oxocarbenium ion, yet an unactivated water molecule is not, and that the electric field catalysis strategy is used by the enzyme to enable its unique double-displacement reaction mechanism. A strong IEF, pointing toward 5' direction of the substrate sugar ring, greatly facilitates the second displacement reaction at the expense of elevating the barrier of the first one, thereby allowing both reactions to occur. These findings not only increase our understanding of the strategies used by DNA glycosylases to repair DNA lesions, but also have important implications for how internal/external electric field can be applied to modulate chemical reactions.

Melatonin Activation by Human Cytochrome P450 Enzymes: A Comparison between Different Isozymes

Cytochrome P450 enzymes in the human body play a pivotal role in both the biosynthesis and the degradation of the hormone melatonin. Melatonin plays a key role in circadian rhythms in the body, but its concentration is also linked to mood fluctuations as well as emotional well-being. In the present study, we present a computational analysis of the binding and activation of melatonin by various P450 isozymes that are known to yield different products and product distributions. In particular, the P450 isozymes 1A1, 1A2, and 1B1 generally react with melatonin to provide dominant aromatic hydroxylation at the C6-position, whereas the P450 2C19 isozyme mostly provides O-demethylation products. To gain insight into the origin of these product distributions of the P450 isozymes, we performed a comprehensive computational study of P450 2C19 isozymes and compared our work with previous studies on alternative isozymes. The work covers molecular mechanics, molecular dynamics and quantum mechanics approaches. Our work highlights major differences in the size and shape of the substrate binding pocket amongst the different P450 isozymes. Consequently, substrate binding and positioning in the active site varies substantially within the P450 isozymes. Thus, in P450 2C19, the substrate is oriented with its methoxy group pointing towards the heme, and therefore reacts favorably through hydrogen atom abstraction, leading to the production of O-demethylation products. On the other hand, the substrate-binding pockets in P450 1A1, 1A2, and 1B1 are tighter, direct the methoxy group away from the heme, and consequently activate an alternative site and lead to aromatic hydroxylation instead.

Catalysis by KDM6 Histone Demethylases - A Synergy between the Non-Heme Iron(II) Center, Second Coordination Sphere, and Long-Range Interactions

KDM6A (UTX) and KDM6B (JMJD3) are human non-heme Fe(II) and 2-oxoglutarate (2OG) dependent JmjC oxygenases that catalyze the demethylation of trimethylated lysine 27 in the N-terminal tail of histone H3, a post-translational modification that regulates transcription. A Combined Quantum Mechanics/ Molecular Mechanics (QM/MM) and Molecular Dynamics (MD) study on the catalytic mechanism of KDM6A/B reveals that the transition state for the rate-limiting hydrogen atom transfer (HAT) reaction in KDM6A catalysis is stabilized by polar (Asn217) and aromatic (Trp369)/non-polar (Pro274) residues in contrast to KDM4, KDM6B and KDM7 demethylases where charged residues (Glu, Arg, Asp) are involved. KDM6A employs both σ- and π-electron transfer pathways for HAT, whereas KDM6B employs the σ-electron pathway. Differences in hydrogen bonding of the Fe-chelating Glu252(KDM6B) contribute to the lower energy barriers in KDM6B vs. KDM6A. The study reveals a dependence of the activation barrier of the rebound hydroxylation on the Fe-O-C angle in the transition state of KDM6A. Anti-correlation of the Zn-binding domain with the active site residues is a key factor distinguishing KDM6A/B from KDM7/4s. The results reveal the importance of communication between the Fe center, second coordination sphere, and long-range interactions in catalysis by KDMs and, by implication, other 2OG oxygenases.

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.

How Is Substrate Halogenation Triggered by the Vanadium Haloperoxidase from Curvularia inaequalis?

Vanadium haloperoxidases (VHPOs) are unique enzymes in biology that catalyze a challenging halogen transfer reaction and convert a strong aromatic C-H bond into C-X (X = Cl, Br, I) with the use of a vanadium cofactor and H₂O₂. The VHPO catalytic cycle starts with the conversion of hydrogen peroxide and halide (X = Cl, Br, I) into hypohalide on the vanadate cofactor, and the hypohalide subsequently reacts with a substrate. However, it is unclear whether the hypohalide is released from the enzyme or otherwise trapped within the enzyme structure for the halogenation of organic substrates. A substrate-binding pocket has never been identified for the VHPO enzyme, which questions the role of the protein in the overall reaction mechanism. Probing its role in the halogenation of small molecules will enable further engineering of the enzyme and expand its substrate scope and selectivity further for use in biotechnological applications as an environmentally benign alternative to current organic chemistry synthesis. Using a combined experimental and computational approach, we elucidate the role of the vanadium haloperoxidase protein in substrate halogenation. Activity studies show that binding of the substrate to the enzyme is essential for the reaction of the hypohalide with substrate. Stopped-flow measurements demonstrate that the rate-determining step is not dependent on substrate binding but partially on hypohalide formation. Using a combination of molecular mechanics (MM) and molecular dynamics (MD) simulations, the substrate binding area in the protein is identified and even though the selected substrates (methylphenylindole and 2-phenylindole) have limited hydrogen-bonding abilities, they are found to bind relatively strongly and remain stable in a binding tunnel. A subsequent analysis of the MD snapshots characterizes two small tunnels leading from the vanadate active site to the surface that could fit small molecules such as hypohalide, halide, and hydrogen peroxide. Density functional theory studies using electric field effects show that a polarized environment in a specific direction can substantially lower barriers for halogen transfer. A further analysis of the protein structure indeed shows a large dipole orientation in the substrate-binding pocket that could enable halogen transfer through an applied local electric field. These findings highlight the importance of the enzyme in catalyzing substrate halogenation by providing an optimal environment to lower the energy barrier for this challenging aromatic halide insertion reaction.

Regioselective oxidation of heterocyclic aromatic hydrocarbons catalyzed by cytochrome P450: A case study of carbazole

Recently there are increasing interests in accurately evaluating the health effects of heterocyclic PAHs. However, the activation mechanism and possible metabolites of heterocyclic PAHs catalyzed by human CYP1A1 is still elusive to a great extent. Here, leveraged to high level QM/MM calculations, the corresponding activation pathways of a representative heterocyclic PAHs, carbazole, were systematically explored. The first stage is electrophilic addition or hydrogen abstraction from N-H group. Electrophilic addition was evidenced to be more feasible and regioselectivity at C3 and C4 sites were identified. Correlations between energy barriers and key structural/electrostatic parameters reveal that O-Cα distance and Fe-O-Cα angle are the main origin for the catalytic regioselectivity. Electrophilic addition was determined as the rate-determining step and the subsequent possible reactions include epoxidation, NIH shift (the hydrogen migration from the site of hydroxylation to the adjacent carbon) and proton shuttle. The corresponding products are epoxides, ketones and hydroxylated carbazoles, respectively. The main metabolites (hydroxylated carbazoles) are estimated to be more toxic than carbazole. The regioselectivity of carbazole activated by CYP1A1 is different from the environmental processes (gas and aqueous phase). Collectively, these results will inform the in-depth understanding the metabolic processes of heterocyclic PAHs and aid the accurate evaluation of their health effects.

Phenylalanine Residues in the Active Site of CYP2E1 Participate in Determining the Binding Orientation and Metabolism-Dependent Genotoxicity of Aromatic Compounds

The composition of amino acids forming the active site of a CYP enzyme is impactful in its substrate selectivity. For CYP2E1, the role of PHE residues in the formation of effective binding orientations for its aromatic substrates remains unclear. In this study, molecular docking and molecular dynamics analysis were performed to reflect the interactions between PHEs in the active site of human CYP2E1 and various aromatic compounds known as its substrates. The results indicated that the orientation of 1-methylpyrene (1-MP) in the active site was highly determined by the presence of PHEs, PHE478 contributing to the binding free energy most significantly. Moreover, by building a random forest model the relationship between each of 19 molecular descriptors of polychlorinated biphenyl (PCB) compounds (from molecular docking, quantum mechanics, and physicochemical properties) and their human CYP2E1-dependent mutagenicityas established mostly in our lab, was investigated. The presence of PHEs did not appear to significantly modify the electronic or structural feature of each bound ligand (PCB), instead, the flexibility of the conformation of PHEs contributed substantially to the effective binding energy and orientation. It is supposed that PHE residues adjust their own conformation to permit a suitablly shaped cavity for holding the ligand and forming its orientation as favorable for a biochemical reaction. This study has provided some insights into the role of PHEs in guiding the interactive adaptation of the active site of human CYP2E1 for the binding and metabolism of aromatic substrates.

Brønsted Acids Promote Olefin Oxidations by Bioinspired Nonheme CoIII(PhIO)(OH) Complexes: A Role for Low-Barrier Hydrogen Bonds

Introduction of Brønsted acids into biomimetic nonheme reactions promotes the oxidative ability of metal-oxygen complexes significantly. However, the molecular machinery of the promoted effects is missing. Herein, a comprehensive investigation of styrene oxidation by a cobalt(III)-iodosylbenzene complex, [(TQA)Co^III(OIPh)(OH)]²⁺ (1, TQA = tris(2-quinolylmethyl)amine), in the presence and absence of triflic acid (HOTf) was performed using density functional theory calculations. Results revealed for the first time that there is a low-barrier hydrogen bond (LBHB) between HOTf and the hydroxyl ligand of 1, which forms two valence-resonance structures [(TQA)Co^III(OIPh)(HO^---HOTf)]²⁺ (1_LBHB) and [(TQA)Co^III(OIPh)(H₂O--OTf^-)]²⁺ (1'_LBHB). Due to the oxo-wall, these complexes (1_LBHB and 1'_LBHB) cannot convert to high-valent cobalt-oxyl species. Instead, styrene oxidation by these oxidants (1_LBHB and 1'_LBHB) shows novel spin-state selectivity, i.e., on the ground closed-shell singlet state, styrene is oxidized to an epoxide, whereas on the excited triplet and quintet states, an aldehyde product, phenylacetaldehyde, is formed. The preferred pathway is styrene oxidation by 1'_LBHB, which is initiated by a rate-limiting bond-formation-coupled electron transfer process with an energy barrier of 12.2 kcal mol^-1. The nascent PhIO-styrene-radical-cation intermediate undergoes an intramolecular rearrangement to produce an aldehyde. The halogen bond between the OH^-/H₂O ligand and the iodine of PhIO modulates the activity of the cobalt-iodosylarene complexes 1_LBHB and 1'_LBHB. These new mechanistic findings enrich our knowledge of nonheme chemistry and hypervalent iodine chemistry and will play a positive role in the rational design of new catalysts.

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.

Bio-Inspired Iron Pentadentate Complexes as Dioxygen Activators in the Oxidation of Cyclohexene and Limonene

The use of dioxygen as an oxidant in fine chemicals production is an emerging problem in chemistry for environmental and economical reasons. In acetonitrile, the [(N4Py)Fe^II]²⁺ complex, [N4Py-N,N-bis(2-pyridylmethyl)-N-(bis-2-pyridylmethyl)amine] in the presence of the substrate activates dioxygen for the oxygenation of cyclohexene and limonene. Cyclohexane is oxidized mainly to 2-cyclohexen-1-one, and 2-cyclohexen-1-ol, cyclohexene oxide is formed in much smaller amounts. Limonene gives as the main products limonene oxide, carvone, and carveol. Perillaldehyde and perillyl alcohol are also present in the products but to a lesser extent. The investigated system is twice as efficient as the [(bpy)₂Fe^II]²⁺/O₂/cyclohexene system and comparable to the [(bpy)₂Mn^II]²⁺/O₂/limonene system. Using cyclic voltammetry, it has been shown that, when the catalyst, dioxgen, and substrate are present simultaneously in the reaction mixture, the iron(IV) oxo adduct [(N4Py)Fe^IV=O]²⁺ is formed, which is the oxidative species. This observation is supported by DFT calculations.

Theoretical investigation on the elusive structure-activity relationship of bioinspired high-valence nickel-halogen complexes in oxidative fluorination reactions

Very recently, bioinspired high-valence metal-halogen complexes have been proven to be competent oxidants in the C-H bond activation and heteroatom dihalogenation reactions. However, the structure-activity relationship of such active species and the reaction mechanisms of oxidations mediated by these oxidants are still elusive. In this study, density functional theory (DFT) calculations were performed to systematically study the oxidizing ability of the high-valence Ni^III-X (X = F and Cl) complexes Et₄N[Ni^III(Cl/F)(L)], (1_Cl/F, Et = ethyl, L = N,N'-(2,6-dimethylphenyl)-2,6-pyridinedicarboxamide), such as the reaction mechanism of fluorination of 1,4-cyclohexadiene (CHD) by 1_F in the presence of AgF and the reaction mechanism of difluorination of triphenyl phosphine (PPh₃) by 1_F. All calculated results fit well with the experiments and present new mechanistic findings. The C-H bond activation by the high-valence nickel(III)-halogen complexes was found to proceed via a hydrogen-atom transfer (HAT) mechanism by analysis of the molecular orbitals of the transition states. C-H bond activation by 1_F takes a Ni-F-H angle of ca. 180°, whereas that by 1_Cl takes an angle of ca. 120° on the transition states. These results indicate that the exchange-enhanced reactivity is responsible for the dramatic oxidative difference between these two oxidants. The role of AgF in C-H fluorination of CHD by 1_F is proposed to act as a Lewis acid adduct, AgF-binding Ni(III)-fluorine complex 1_F-Ag-F, which acts both as an oxidant in C-H bond activation and as a fluorine donor in the fluorination step. A cooperative oxidation mechanism involving two 1_F oxidants was proposed for the difluorination of PPh₃ by 1_F. These theoretical findings will enrich the knowledge of high-valence metal-halogen chemistry and play a positive role in the rational design of new catalysts.

Computing the Differences between Asn-X and Gln-X Deamidation and Their Impact on Pharmaceutical and Physiological Proteins: A Theoretical Investigation Using Model Dipeptides

Protein deamidation is a degradation mechanism that significantly impacts both pharmaceutical and physiological proteins. Deamidation impacts two amino acids, Asn and Gln, where the net neutral residues are converted into their acidic forms. While there are multiple similarities between the reaction mechanisms of the two residues, the impact of Gln deamidation has been noted to be most significant on physiological proteins while Asn deamidation has been linked to both pharmaceutical and physiological proteins. For this purpose, we sought to analyze the thermochemical and kinetic properties of the different reactions of Gln deamidation relative to Asn deamidation. In this study, we mapped the deamidation of Gln-X dipeptides into Glu-X dipeptides using density functional theory (DFT). Full network mapping facilitated the prediction of reaction selectivity between the two primary pathways, as well as between the two products of Gln-X deamidation as a function of solvent dielectric. To achieve this analysis, we studied a total of 77 dipeptide reactions per solvent dielectric (308 total reactions). Modeled at a neutral pH and using quantum chemical and statistical thermodynamic methods, we computed the following values: enthalpy of reaction (ΔH_RXN), entropy (ΔS_RXN), Gibbs free energy of reaction (ΔG_RXN), activation energy (E_A), and the Arrhenius preexponential factor (log(A)) for each dipeptide. Additionally, using chemical reaction principles, we generated a database of computed rate coefficients for all possible N-terminus Gln-X deamidation reactions at a neutral pH, predicted the most likely deamidation reaction mechanism for each dipeptide reaction, analyzed our results against our prior study on Asn-X deamidation, and matched our results against qualitative trends previously noted by experimental literature.

The Performance of Different Water Models on the Structure and Function of Cytochrome P450 Enzymes

Modeling approaches and modern simulations to investigate the biomolecular structure and function rely on various methods. Since water molecules play a crucial role in all sorts of chemistry, the accurate modeling of water molecules is vital for such simulations. In cytochrome P450 (CYP450), in particular, water molecules play a key role in forming active oxidant that ultimately performs oxidation and metabolism. In the present study, we have highlighted the behavior of the three most widely used water models─TIP3P, SPC/E, and OPC─for three different CYP450 enzymes─CYP450_BM3, CYP450_OleT, and CYP450_BSβ─during MD simulations and QM/MM calculations. We studied the various properties, such as RMSD, RMSF, H-bond, water occupancy, and hydrogen atom transfer (HAT), using QM/MM calculations and compared them for all three water models. Our study shows that the stabilities of the enzyme complexes are well maintained in all three water models. However, the OPC water model performs well for the polar active sites, that is, in CYP450_OleT and CYP450_BSβ, while the TIP3P water model is superior for the hydrophobic site, such as CYP450_BM3.

Can Second Coordination Sphere and Long-Range Interactions Modulate Hydrogen Atom Transfer in a Non-Heme Fe(II)-Dependent Histone Demethylase?

Fe(II)-dependent oxygenases employ hydrogen atom transfer (HAT) to produce a myriad of products. Understanding how such enzymes use dynamic processes beyond the immediate vicinity of the active site to control the selectivity and efficiency of HAT is important for metalloenzyme engineering; however, obtaining such knowledge by experiments is challenging. This study develops a computational framework for identifying second coordination sphere (SCS) and especially long-range (LR) residues relevant for catalysis through dynamic cross-correlation analysis (DCCA) using the human histone demethylase PHF8 (KDM7B) as a model oxygenase. Furthermore, the study explores the mechanistic pathways of influence of the SCS and LR residues on the HAT reaction. To demonstrate the plausibility of the approach, we investigated the effect of a PHF8 F279S clinical mutation associated with X-linked intellectual disability, which has been experimentally shown to ablate PHF8-catalyzed demethylation. In agreement, the molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) studies showed a change in the H3_1-14K9me2 substrate orientation and an increased HAT barrier. We systematically analyzed the pathways by which the identified SCS and LR residues may influence HAT by exploring changes in H3K9me2 substrate orientation, interdomain correlated motions, HAT transition state stabilization, reaction energetics, electron transfer mechanism, and alterations in the intrinsic electric field of PHF8. Importantly, SCS and LR variations decrease key motions of α9-α12 of the JmjC domain toward the Fe(IV)-center that are associated with tighter binding of the H3_1-14K9me2 substrate. SCS and LR residues alter the intrinsic electric field of the enzyme along the reaction coordinate and change the individual energetic contributions of residues toward TS stabilization. The overall results suggest that DCCA can indeed identify non-active-site residues relevant for catalysis. The substitutions of such dynamically correlated residues might be used as a tool to tune HAT in non-heme Fe(II)- and 2OG-dependent enzymes.

Enzymatic Control over Reactive Intermediates Enables Direct Oxidation of Alkenes to Carbonyls by a P450 Iron-Oxo Species

The aerobic oxidation of alkenes to carbonyls is an important and challenging transformation in synthesis. Recently, a new P450-based enzyme (aMOx) has been evolved in the laboratory to directly oxidize styrenes to their corresponding aldehydes with high activity and selectivity. The enzyme utilizes a heme-based, high-valent iron-oxo species as a catalytic oxidant that normally epoxidizes alkenes, similar to other catalysts. How the evolved aMOx enzyme suppresses the commonly preferred epoxidation and catalyzes direct carbonyl formation is currently not well understood. Here, we combine computational modelling together with mechanistic experiments to study the reaction mechanism and unravel the molecular basis behind the selectivity achieved by aMOx. Our results describe that although both pathways are energetically accessible diverging from a common covalent radical intermediate, intrinsic dynamic effects determine the strong preference for epoxidation. We discovered that aMOx overrides these intrinsic preferences by controlling the accessible conformations of the covalent radical intermediate. This disfavors epoxidation and facilitates the formation of a carbocation intermediate that generates the aldehyde product through a fast 1,2-hydride migration. Electrostatic preorganization of the enzyme active site also contributes to the stabilization of the carbocation intermediate. Computations predicted that the hydride migration is stereoselective due to the enzymatic conformational control over the intermediate species. These predictions were corroborated by experiments using deuterated styrene substrates, which proved that the hydride migration is cis- and enantioselective. Our results demonstrate that directed evolution tailored a highly specific active site that imposes strong steric control over key fleeting biocatalytic intermediates, which is essential for accessing the carbonyl forming pathway and preventing competing epoxidation.

Oriented External Electric Fields Regurating the Reaction Mechanism of CH4 Oxidation Catalyzed by Fe(IV)-Oxo-Corrolazine: Insight from Density Functional Calculations

Methane is the simplest alkane and can be used as an alternative energy source for oil and coal, but the greenhouse effect caused by its leakage into the air is not negligible, and its conversion into liquid methanol not only facilitates transportation, but also contributes to carbon neutrality. In order to find an efficient method for converting methane to methanol, CH4 oxidation catalyzed by Fe(IV)-Oxo-corrolazine (Fe(IV)-Oxo-Cz) and its reaction mechanism regulation by oriented external electric fields (OEEFs) are systematically studied by density functional calculations. The calculations show that Fe(IV)-Oxo-Cz can abstract one H atom from CH4 to form the intermediate with OH group connecting on the corrolazine ring, with the energy barrier of 25.44 kcal mol-1. And then the product methanol is formed through the following rebound reaction. Moreover, the energy barrier can be reduced to 20.72 kcal mol-1 through a two-state reaction pathway. Furthermore, the effect of OEEFs on the reaction is investigated. We found that OEEFs can effectively regulate the reaction by adjusting the stability of the reactant and the transition state through the interaction of electric field-molecular dipole moment. When the electric field is negative, the energy barrier of the reaction decreases with the increase of electric intensity. Moreover, the OEEF aligned along the intrinsic Fe‒O reaction axis can effectively regulate the ability of forming the OH on the corrolazine ring by adjusting the charges of O and H atoms. When the electric field intensity is -0.010 a.u., the OH can be directly rebounded to the CH3· before it is connecting on the corrolazine ring, thus forming the product directly from the transition state without passing through the intermediate with only an energy barrier of 17.34 kcal mol-1, which greatly improves the selectivity of the reaction.

Enhancing regioselectivity of sucrose phosphorylase by loop engineering for glycosylation of L-ascorbic acid

Sucrose phosphorylase (SPase) has a remarkable capacity to synthesize numerous glucosides from abundantly available sucrose under mild conditions but suffers from specificity and regioselectivity issues. In this study, a loop engineering strategy was introduced to enhance the regioselectivity and substrate specificity of SPase for the efficient synthesis of 2-O-α-D-glucopyranosyl-L-ascorbic acid (AA-2G) via L-ascorbic acid (L-AA). P134, L341, and L343 were identified as "hotspots" for modulating the flexibility of loops, which significantly influenced the H-bonding network of L-AA in the active site, as well as the entrance of the substrate channel, thereby altering the regioselectivity and substrate specificity. Finally, the mutant L341V/L343F, with near-perfect control of the selectivity synthesis of the 2-OH group of L-AA (> 99%), was obtained. The AA-2G production by the mutant reached 244 g L^-1 in a whole-cell biotransformation system, and the conversion rate of L-AA reached 64%, which is the highest level reported to date. Our work also provides a successful loop engineering case for modulating the regioselectivity and specificity of sucrose phosphorylase. KEY POINTS: • "Hotspots" were identified in the flexible loops of sucrose phosphorylase. • Mutants exhibited improved regioselectivity and specificity against L-ascorbic acid. • Synthesized AA-2G with high yield and regioselectivity by whole-cell of mutant.

Second coordination sphere effects on the mechanistic pathways for dioxygen activation by a ferritin: involvement of a Tyr radical and the identification of a cation binding site

Ferritins are ubiquitous diiron enzymes involved in iron(II) detoxification and oxidative stress responses and can act as metabolic iron stores. The overall reaction mechanisms of ferritin enzymes are still unclear, particularly concerning the role of the conserved, near catalytic center Tyr residue. Thus, we carried out a computational study of a ferritin using a large cluster model of well over 300 atoms including its first- and second-coordination sphere. The calculations reveal important insight into the structure and reactivity of ferritins. Specifically, the active site Tyr residue delivers a proton and electron in the catalytic cycle prior to iron(II) oxidation. In addition, the calculations highlight a likely cation binding site at Asp65, which through long-range electrostatic interactions, influences the electronic configuration and charge distributions of the metal center. The results are consistent with experimental observations but reveal novel detail of early mechanistic steps that lead to an unusual mixed-valent iron(III)-iron(II) center.

Which is the real oxidant in the competitive ligand self-hydroxylation and substrate oxidation, a biomimetic iron(II)-hydroperoxo species or an oxo-iron(IV)-hydroxy one?

Nonheme iron(II)-hydroperoxo species (FeII-(η2-OOH)) 1 and the concomitant oxo-iron(IV)-hydroxyl one 2 are proposed as the key intermediates of a large class of 2-oxoglutarate dependent dioxygenases (e.g., isopenicillin N synthase). Extensive...

Can the isonitrile biosynthesis enzyme ScoE assist with the biosynthesis of isonitrile groups in drug molecules? A computational study

Computational studies show that the isonitrile synthesizing enzyme ScoE can catalyse the conversion of γ-Gly substituents in substrates to isonitrile. This enables efficient isonitrile substitution into target molecules such as axisonitrile-1.