Electrostatic Perturbations in the Substrate-Binding Pocket of Taurine/α-Ketoglutarate Dioxygenase Determine its Selectivity

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.

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.

Rational Design of the Spatial Effect in a Fe(II)/α-Ketoglutarate-Dependent Dioxygenase Reverses the Regioselectivity of C(sp3)-H Bond Hydroxylation in Aliphatic Amino Acids

The hydroxylation of remote C(sp3)-H bonds in aliphatic amino acids yields crucial precursors for the synthesis of high-value compounds. However, accurate regulation of the regioselectivity of remote C(sp3)-H bonds hydroxylation in aliphatic amino acids continues to be a common challenge in chemosynthesis and biosynthesis. In this study, the Fe(II)/α-ketoglutarate-dependent dioxygenase from Bacillus subtilis (BlAH) was mined and found to catalyze hydroxylation at the γ and δ sites of aliphatic amino acids. Crystal structure analysis, molecular dynamics simulations, and quantum chemical calculations revealed that regioselectivity was regulated by the spatial effect of BlAH. Based on these results, the spatial effect of BlAH was reconstructed to stabilize the transition state at the δ site of aliphatic amino acids, thereby successfully reversing the γ site regioselectivity to the δ site. For example, the regioselectivity of L-Homoleucine (5 a) was reversed from the γ site (1 : 12) to the δ site (>99 : 1). The present study not only expands the toolbox of biocatalysts for the regioselective functionalization of remote C(sp3)-H bonds, but also provides a theoretical guidance for the precision-driven modification of similarly remote C(sp3)-H bonds in complex molecules.

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.

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.

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.

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.

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.

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.

Identification of isobenzofuranone derivatives as promising antidiabetic agents: Synthesis, in vitro and in vivo inhibition of α-glucosidase and α-amylase, computational docking analysis and molecular dynamics simulations

Diabetes mellitus, one of the major health challenges of the 21st century, is associated with numerous biomedical complications including retinopathy, neuropathy, nephropathy, cardiovascular diseases and liver disorders. To control the chronic hyperglycemic condition, the development of potential inhibitors of drug targets such as α-glucosidase and α-amylase remains a promising strategy and focus of continuous efforts. Therefore, in the present work, a concise library of isobenzofuranone derivatives (3a-q) was designed and synthesized using Suzuki-Miyaura cross-coupling approach. The biological potential of these heterocyclic compounds against carbohydrate-hydrolyzing enzymes; α-glucosidase and α-amylase, was examined. In vitro inhibitory results demonstrated that the tested isobenzofuranones were considerably more effective and potent inhibitors than the standard drug, acarbose. Compound 3d having an IC50 value of 6.82 ± 0.02 μM was emerged as the lead candidate against α-glucosidase with ⁓127-folds strong inhibition than acarbose. Similarly, compound 3g demonstrated ⁓11-folds higher inhibition strength against α-amylase when compared with acarbose. Both compounds were tested in vivo and results demonstrate that the treatment of diabetic rats with α-amylase inhibitor show more pronounced histopathological normalization in kidney and liver than with α-glucosidase inhibitor. The Lineweaver-Burk plot revealed an uncompetitive mode of inhibition for 3d against α-glucosidase whereas compound 3g exhibited mixed inhibition against α-amylase. Furthermore, in silico molecular docking and dynamics simulations validated the in vitro data for these compounds whereas pharmacokinetics profile revealed the druglike properties of potent inhibitors.

Discovery of druggable potent inhibitors of serine proteases and farnesoid X receptor by ligand-based virtual screening to obstruct SARS-CoV-2

The coronavirus, a subfamily of the coronavirinae family, is an RNA virus with over 40 variations that can infect humans, non-human mammals and birds. There are seven types of human coronaviruses, including SARS-CoV-2, is responsible for the recent COVID-19 pandemic. The current study is focused on the identification of drug molecules for the treatment of COVID-19 by targeting human proteases like transmembrane serine protease 2 (TMPRSS2), furin, cathepsin B, and a nuclear receptor named farnesoid X receptor (FXR). TMPRSS2 and furin help in cleaving the spike protein of the SARS-CoV-2 virus, while cathepsin B plays a critical role in the entry and pathogenesis. FXR, on the other hand, regulates the expression of ACE2, and its inhibition can reduce SARS-CoV-2 infection. By inhibiting these four protein targets with non-toxic inhibitors, the entry of the infectious agent into host cells and its pathogenesis can be obstructed. We have used the BioSolveIT suite for pharmacophore-based computational drug designing. A total of 1611 ligands from the ligand library were docked with the target proteins to obtain potent inhibitors on the basis of pharmacophore. Following the ADMET analysis and protein ligand interactions, potent and druggable inhibitors of the target proteins were obtained. Additionally, toxic substructures and the less toxic route of administration of the most potent inhibitors in rodents were also determined computationally. Compounds namely N-(diaminomethylene)-2-((3-((1R,3R)-3-(2-(methoxy(methyl)amino)-2-oxoethyl)cyclopentyl)propyl)amino)-2-oxoethan-1-aminium (26), (1R,3R)-3-(((2-ammonioethyl)ammonio)methyl)-1-((4-propyl-1H-imidazol-2-yl)methyl)piperidin-1-ium (29) and (1R,3R)-3-(((2-ammonioethyl)ammonio)methyl)-1-((1-propyl-1H-pyrazol-4-yl)methyl)piperidin-1-ium (30) were found as the potent inhibitors of TMPRSS2, whereas, 1-(1-(1-(1H-tetrazol-1-yl)cyclopropane-1‑carbonyl)piperidin-4-yl)azepan-2-one (6), (2R)-4-methyl-1-oxo-1-((7R,11S)-4-oxo-6,7,8,9,10,11-hexahydro-4H-7,11-methanopyrido[1,2-a]azocin-9-yl)pentan-2-aminium (12), 4-((1-(3-(3,5-dimethylisoxazol-4-yl)propanoyl)piperidin-4-yl)methyl)morpholin-4-ium (13), 1-(4,6-dimethylpyrimidin-2-yl)-N-(3-oxocyclohex-1-en-1-yl)piperidine-4-carboxamide (14), 1-(4-(1,5-dimethyl-1H-1,2,4-triazol-3-yl)piperidin-1-yl)-3-(3,5-dimethylisoxazol-4-yl)propan-1-one (25) and N,N-dimethyl-4-oxo-4-((1S,5R)-8-oxo-5,6-dihydro-1H-1,5-methanopyrido[1,2-a][1,5]diazocin-3(2H,4H,8H)-yl)butanamide (31) inhibited the FXR preferentially. In case of cathepsin B, N-((5-benzoylthiophen-2-yl)methyl)-2-hydrazineyl-2-oxoacetamide (2) and N-([2,2'-bifuran]-5-ylmethyl)-2-hydrazineyl-2-oxoacetamide (7) were identified as the most druggable inhibitors whereas 1-amino-2,7-diethyl-3,8-dioxo-6-(p-tolyl)-2,3,7,8-tetrahydro-2,7-naphthyridine-4‑carbonitrile (5) and (R)-6-amino-2-(2,3-dihydroxypropyl)-1H-benzo[de]isoquinoline-1,3(2H)-dione (20) were active against furin.

Identification of potential inhibitors targeting yellow fever virus helicase through ligand and structure-based computational studies

Yellow fever is a flavivirus having plus-sensed RNA which encodes a single polyprotein. Host proteases cut this polyprotein into seven nonstructural proteins including a vital NS3 protein. The present study aims to identify the most effective inhibitor against the helicase (NS3) using different advanced ligand and structure-based computational studies. A set of 300 ligands was selected against helicase by chemical structural similarity model, which are similar to S-adenosyl-l-cysteine using infiniSee. This tool screens billions of compounds through a similarity search from in-built chemical spaces (CHEMriya, Galaxi, KnowledgeSpace and REALSpace). The pharmacophore was designed from ligands in the library that showed same features. According to the sequence of ligands, six compounds (29, 87, 99, 116, 148, and 208) were taken for pharmacophore designing against helicase protein. Subsequently, compounds from the library which showed the best pharmacophore shared-features were docked using FlexX functionality of SeeSAR and their optibrium properties were analyzed. Afterward, their ADME was improved by replacing the unfavorable fragments, which resulted in the generation of new compounds. The selected best compounds (301, 302, 303 and 304) were docked using SeeSAR and their pharmacokinetics and toxicological properties were evaluated using SwissADME. The optimal inhibitor for yellow fever helicase was 2-amino-N-(4-(dimethylamino)thiazol-2-yl)-4-methyloxazole-5-carboxamide (302), which exhibits promising potential for drug development.Communicated by Ramaswamy H. Sarma.

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.

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.

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.

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.

Melatonin Activation by Cytochrome P450 Isozymes: How Does CYP1A2 Compare to CYP1A1?

Cytochrome P450 enzymes are versatile enzymes found in most biosystems that catalyze mono-oxygenation reactions as a means of biosynthesis and biodegradation steps. In the liver, they metabolize xenobiotics, but there are a range of isozymes with differences in three-dimensional structure and protein chain. Consequently, the various P450 isozymes react with substrates differently and give varying product distributions. To understand how melatonin is activated by the P450s in the liver, we did a thorough molecular dynamics and quantum mechanics study on cytochrome P450 1A2 activation of melatonin forming 6-hydroxymelatonin and N-acetylserotonin products through aromatic hydroxylation and O-demethylation pathways, respectively. We started from crystal structure coordinates and docked substrate into the model, and obtained ten strong binding conformations with the substrate in the active site. Subsequently, for each of the ten substrate orientations, long (up to 1 μs) molecular dynamics simulations were run. We then analyzed the orientations of the substrate with respect to the heme for all snapshots. Interestingly, the shortest distance does not correspond to the group that is expected to be activated. However, the substrate positioning gives insight into the protein residues it interacts with. Thereafter, quantum chemical cluster models were created and the substrate hydroxylation pathways calculated with density functional theory. These relative barrier heights confirm the experimental product distributions and highlight why certain products are obtained. We make a detailed comparison with previous results on CYP1A1 and identify their reactivity differences with melatonin.

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 C⁴-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.

An iron(ii)-based metalloradical system for intramolecular amination of C(sp2)-H and C(sp3)-H bonds: synthetic applications and mechanistic studies

A catalytic system for intramolecular C(sp2)-H and C(sp3)-H amination of substituted tetrazolopyridines has been successfully developed. The amination reactions are developed using an iron-porphyrin based catalytic system. It has been demonstrated that the same iron-porphyrin based catalytic system efficiently activates both the C(sp2)-H and C(sp3)-H bonds of the tetrazole as well as azide-featuring substrates with a high level of regioselectivity. The method exhibited an excellent functional group tolerance. The method affords three different classes of high-value N-heterocyclic scaffolds. A number of important late-stage C-H aminations have been performed to access important classes of molecules. Detailed studies (experimental and computational) showed that both the C(sp2)-H and C(sp3)-H amination reactions involve a metalloradical activation mechanism, which is different from the previously reported electro-cyclization mechanism. Collectively, this study reports the discovery of a new class of metalloradical activation modes using a base metal catalyst that should find wide application in the context of medicinal chemistry, drug discovery and industrial applications.

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.

What Drives Radical Halogenation versus Hydroxylation in Mononuclear Nonheme Iron Complexes? A Combined Experimental and Computational Study

Nonheme iron halogenases are unique enzymes in nature that selectively activate an aliphatic C–H bond of a substrate to convert it into C–X (X = Cl/Br, but not F/I). It is proposed that they generate an Fe^III(OH)(X) intermediate in their catalytic cycle. The analogous Fe^III(OH) intermediate in nonheme iron hydroxylases transfers OH^• to give alcohol product, whereas the halogenases transfer X^• to the carbon radical substrate. There remains significant debate regarding what factors control their remarkable selectivity of the halogenases. The reactivity of the complexes Fe^III(BNPA^Ph2O)(OH)(X) (X = Cl, Br) with a secondary carbon radical (R^•) is described. It is found that X^• transfer occurs with a secondary carbon radical, as opposed to OH^• transfer with tertiary radicals. Comprehensive computational studies involving density functional theory were carried out to examine the possible origins of this selectivity. The calculations reproduce the experimental findings, which indicate that halogen transfer is not observed for the tertiary radicals because of a nonproductive equilibrium that results from the endergonic nature of these reactions, despite a potentially lower reaction barrier for the halogenation pathway. In contrast, halogen transfer is favored for secondary carbon radicals, for which the halogenated product complex is thermodynamically more stable than the reactant complex. These results are rationalized by considering the relative strengths of the C–X bonds that are formed for tertiary versus secondary carbon centers. The computational analysis also shows that the reaction barrier for halogen transfer is significantly dependent on secondary coordination sphere effects, including steric and H-bonding interactions.

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.

Critical Roles of Exchange and Superexchange Interactions in Dictating Electron Transfer and Reactivity in Metalloenzymes

Electron transfer (ET) is a fundamental process in transition-metal-dependent metalloenzymes. In these enzymes, the spin-spin interactions within the same metal center and/or between different metal sites can play a pivotal role in the catalytic cycle and reactivity. This Perspective highlights that the exchange and/or superexchange interactions can intrinsically modulate the inner-sphere and long-range electron transfer, thus controlling the mechanism and activity of metalloenzymes. For mixed-valence diiron oxygenases, the spin-regulated inner-sphere ET can be dictated by exchange interactions, leading to efficient O-O bond activations. Likewise, the spin-regulated inner-sphere ET can be enhanced by both exchange and superexchange interactions in [Fe4S4]-dependent SAM enzymes, which enable the efficient cleavage of the S─C(γ) or S─C5' bond of SAM. In addition to inner-sphere ET, superexchange interactions may modulate the long-range ET between metalloenzymes. We anticipate that the exchange and superexchange enhanced reactivity could be applicable in other important metalloenzymes, such as Photosystem II and nitrogenases.

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.