Drug metabolism by cytochrome p450 enzymes: what distinguishes the pathways leading to substrate hydroxylation over desaturation?

Study of phenanthrenes from their unique mass spectrometric behavior through quantum chemical calculations to liquid chromatographic quantitation

Phenanthrenes and their derivatives have biological relevance owing to their antimicrobial, antioxidant, and cytotoxic effects on cancer cells. They can be efficiently analyzed through ultrahigh-performance liquid chromatography coupled to tandem high-resolution mass spectrometry (UHPLC-MS/HRMS). Herein, we first studied the unique fragmentation behavior of phenanthrenes based on direct infusion MS/HRMS analysis. As a newly described phenomenon, "organ pipe distribution", we found a structural connection linking their unique fragmentation pattern to serial H radical losses. The bonds responsible for this behavior were identified through quantum chemical calculations using a stepwise approach. Furthermore, the chromatographic aspect of this study was enhanced by developing, validating, and applying a new unscheduled targeted UHPLC-MS/HRMS method for quantifying phenanthrenes in Juncus compressus herb. Targeted compounds were efficiently separated within 4 min upon utilizing the Accucore C30 column, and the unscheduled targeted analytical approach afforded five new isomers. Compounds 1 (effususol), 3 (dehydroeffusol), and 6 (7-hydroxy-1-methyl-2-methoxy-5-vinyl-9,10-dihydrophenanthrene) had their linearity limits determined within 10-5000 nM, and Compounds 2 (effusol), 4 (juncusol), 5 (effususin A), and 7 (compressin A) within 25-5000 nM. The coefficients of variation for precision ranged from 1.4 % to 15.2 %. The obtained matrix effects and accuracy values were also within acceptable ranges. Compounds 2 (effusol) and 3 (dehydroeffusol) were present in both methanolic and dichloromethanolic extracts of Plants 1 and 3 at the highest concentrations. Furthermore, the relationship between phenanthrene fingerprints, obtained through ANOVA statistical analysis of quantitative data, and the geographical location of herbs was also established.

Evidence for the formation of 6PPD-quinone from antioxidant 6PPD by cytochrome P450

N-(1,3-dimethylbutyl)-N'-phenyl-p-phenylenediamine (6PPD) as a rubber antioxidant has attracted global concern, since its ozone-oxidation product 6PPD-quinone (6PPDQ) was found to be the primary toxicant responsible for urban runoff mortality syndrome in coho salmon. However, the biotransformation fate and associated toxicological mechanism of 6PPD have not received much study yet. In this work, the in vitro assays showed 6PPD can be transformed into 6PPDQ by cytochromes P450 (CYP450) in human liver microsomes (HLMs) with 0.98 % production rate, and the adducts of 6PPDQ with calf thymus DNA and the N-N coupling product between 6PPD and 6PPDQ were further identified after 6PPD incubation in HLMs. Further evidence for the 6PPDQ formation can be obtained from the in vivo assays that the 6PPDQ-DNA adducts and 6PPD-N-N-6PPDQ dimer were detected in mice by oral gavage with 6PPD, and the latter dimer species was detected as well in 6PPD exposure to zebrafish larvae. Especially, the bioaccumulation property and high reactivity of 6PPDQ result in the continuous formation of the significant DNA adducts and 6PPD-N-N-6PPDQ dimer even in case of low production rate of biotransformation of 6PPD to 6PPDQ, which may provide potentially effective biomarkers for such process. DFT computations revealed the formation mechanism of 6PPDQ is the (N)H-abstraction of 6PPD by CYP450, followed by amino radical rebound at the nearby ortho-carbon, yielding a quinol intermediate due to spin delocalization, that might readily undergo further oxidation by CYP450 into 6PPDQ.

Insight into chemically reactive metabolites of aliphatic amine pollutants: A de novo prediction strategy and case study of sertraline

The uncommon metabolic pathways of organic pollutants are easily overlooked, potentially leading to idiosyncratic toxicity. Prediction of their biotransformation associated with the toxic effects is the very purpose that this work focuses, to develop a de novo method to mechanistically predict the reactive toxicity pathways of uncommon metabolites from start aliphatic amine molecules, which employed sertraline triggered by CYP450 enzymes as a model system, as there are growing concerns about the effects on human health posed by antidepressants in the aquatic environment. This de novo prediction strategy combines computational and experimental methods, involving DFT calculations upon sequential growth, in vitro and in vivo assays, dissecting chemically reactive mechanism relevant to toxicity, and rationalizing the fundamental factors. Significantly, desaturation and debenzylation-aromatization as the emerging metabolic pathways of sertraline have been elucidated, with the detection of DNA adducts of oxaziridine metabolite in mice, highlighting the potential reactive toxicity. Molecular orbital analysis supports the reactivity preference for toxicological-relevant C-N desaturation over N-hydroxylation of sertraline, possibly extended to several other aliphatic amines based on the Bell-Evans-Polanyi principle. It was further validated toward some other wide-concerned aliphatic amine pollutants involving atrazine, ε-caprolactam, 6PPD via in silico and in vitro assays, thereby constituting a complete path for de novo prediction from case study to general applications.

Bioconcentration, biotransformation, and transcriptomic impact of the UV-filter 4-MBC in the manila clam Ruditapes philippinarum

Ultraviolet filters (UV-filters) are compounds extensively used in personal care products. These compounds are produced at increasing rates and discharged into marine ecosystems in unknown quantities and with no regulation, making them emerging contaminants. Among those, the UV-filter 4-Methylbenzylidene camphor (4-MBC) is used in a variety of personal care products such as sunscreens, soaps, or lipsticks. This high consumption has resulted in its presence in various environmental matrices at in concentrations ranging from ng to μg L-1. Very little is known, however, about the possible adverse effects in exposed non-target organisms. Our study presents novel data on the bioconcentration, toxicokinetics, and molecular effects of 4-MBC in a marine bivalve species of commercial interest, Ruditapes philippinarum (Manila clam). Organisms were exposed at two different concentrations (1.34 and 10.79 μg L-1) of 4-MBC for 7 days, followed by a 3-day depuration period (clean sea waters). Bioconcentration factors (BCF) were 3562 and 2229 L kg-1 for the low and high exposure concentrations, respectively, making this pollutant bioaccumulative according to REACH criteria. Up to six 4-MBC biotransformation products (BTPs)were identified, 2 of them for the first time. Transcriptomic analysis revealed between 658 and 1310 differently expressed genes (DEGs) after 4-MBC exposure. Functional and enrichment analysis of the DEGs showed the activation of the detoxification pathway to metabolize and excrete the bioconcentrated 4-MBC, which also involved energy depletion and caused an impact on the metabolism of carbohydrates and lipids and in the oxidative phosphorylation pathways. Oxidative stress and immune response were also evidenced through the activation of cathepsins and the complement system. Such elucidation of the mode of action of a ubiquitous pollutant such as 4-MBC at the molecular level is valuable both from an environmental point of view and for the sustainable production of Manila clam, one of the most cultivated mollusk species worldwide.

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.

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.

Emerging Metabolic Profiles of Sulfonamide Antibiotics by Cytochromes P450: A Computational-Experimental Synergy Study on Emerging Pollutants

Metabolism, especially by CYP450 enzymes, is the main reason for mediating the toxification and detoxification of xenobiotics in humans, while some uncommon metabolic pathways, especially for emerging pollutants, probably causing idiosyncratic toxicity are easily overlooked. The pollution of sulfonamide antibiotics in aqueous system has attracted increasing public attention. Hydroxylation of the central amine group can trigger a series of metabolic processes of sulfonamide antibiotics in humans; however, this work parallelly reported the coupling and fragmenting initiated by amino H-abstraction of sulfamethoxazole (SMX) catalyzed by human CYP450 enzymes. Elucidation of the emerging metabolic profiles was mapped via a multistep synergy between computations and experiments, involving preliminary DFT computations and in vitro and in vivo assays, profiling adverse effects, and rationalizing the fundamental factors via targeted computations. Especially, the confirmed SMX dimer was shown to potentially act as a metabolism disruptor in humans, while spin aromatic delocalization resulting in the low electron donor ability of amino radicals was revealed as the fundamental factor to enable coupling of sulfonamide antibiotics by CYP450 through the nonconventional nonrebound pathway. This work may further strengthen the synergistic use of computations prior to experiments to avoid wasteful experimental screening efforts in environmental chemistry and toxicology.

Theoretical Kinetic Isotope Effects in Establishing the Precise Biodegradation Mechanisms of Organic Pollutants

Compound-specific isotope analysis (CSIA) for natural isotope ratios has been recognized as a promising tool to elucidate biodegradation pathways of organic pollutants by microbial enzymes by relating reported kinetic isotope effects (KIEs) to apparent KIEs (AKIEs) derived from bulk isotope fractionations (ε_bulk). However, for many environmental reactions, neither are the reference KIE ranges sufficiently narrow nor are the mechanisms elucidated to the point that rate-determining steps have been identified unequivocally. In this work, besides providing reference KIEs and rationalizing AKIEs, good relationships have been explained by DFT computations for diverse biodegradation pathways with known enzymatic models between the theoretical isotope fractionations (ε_bulk^') from intrinsic KIEs on the rate-determining steps and the observed ε_bulk. (1) To confirm the mechanistic details of previously reported pathway-dependent CSIA, it includes isotope changes in MTBE biodegradation between hydroxylation by CYP450 and S_N2 reaction by cobalamin-dependent methyltransferase, the regioselectivity of toluene biodegradation by CYP450, and the rate-determining step in toluene biodegradation by benzylsuccinate synthase. (2) To yield new fundamental insights into some unclear biodegradation pathways, it consists of the oxidative function of toluene dioxygenase in biodegradation of TCE, the epoxidation mode in biodegradation of TCE by toluene 4-monooxygenase, and the weighted average mechanism in biodegradation of cDCE by CYP450.

Exploring the Factors which Result in Cytochrome P450 Catalyzed Desaturation Versus Hydroxylation

The cytochrome P450 family of monooxygenase enzymes have essential biological roles involving the selective oxidation of carbon-hydrogen bonds. They can also catalyze other important metabolic reactions including desaturation to form alkenes. Currently the factors that control the partition between P450 hydroxylation and desaturation pathways are poorly defined. The CYP199A4 enzyme from the bacterium Rhodopseudomonas palustris HaA2 catalyzes the oxidation of 4-ethyl- and 4-isopropyl- benzoic acids with hydroxylation and desaturation occurring in significant quantities. Here we demonstrate that 4-cyclopropylbenzoic acid is regioselectively hydroxylated by CYP199A4 at the benzylic carbon. In contrast, the oxidation of 4-n-propylbenzoic acid by CYP199A4 results in three major metabolites: an alkene from desaturation and two hydroxylation products at the benzylic (Cα) and Cβ carbons in similar quantities. Extending the length of the alkyl substituent resulted in 4-n-butylbenzoic acid being oxidized at the benzylic position (45%) and desaturated (55%). In contrast, 4-isobutylbenzoic generated very little alkene (5%) but was hydroxylated at the benzylic position (54%) and at the tertiary Cβ position (41%). The oxidation of 4-n-propylbenzoic acid by the F298 V mutant of CYP199A4 occurred with no hydroxylation at Cβ and a significant increase in metabolites arising from desaturation (73%). The X-ray crystal structures of CYP199A4 with each substrate revealed that they bind in the active site with the alkyl substituent positioned over the heme. However, the longer alkylbenzoic acids were bound in a different conformation as was 4-n-propylbenzoic acid in the F298 V mutant. Overall, the changes in metabolite distribution could be ascribed to bond strength differences and the position of the alkyl group relative to the heme.

Mechanism of Melatonin Metabolism by CYP1A1: What Determines the Bifurcation Pathways of Hydroxylation versus Deformylation?

Melatonin, a widely applied cosmetic active ingredient, has a variety of uses as a skin protector through antioxidant and anti-inflammatory functions as well as giving the body UV-induced defenses and immune system support. In the body, melatonin is synthesized from a tryptophan amino acid in a cascade of reactions, but as melatonin is toxic at high concentrations, it is metabolized in the human skin by the cytochrome P450 enzymes. The P450s are diverse heme-based mono-oxygenases that catalyze oxygen atom-transfer processes that trigger metabolism and detoxification reactions in the body. In the catalytic cycle of the P450s, a short-lived high-valent iron(IV)-oxo heme cation radical is formed that has been proposed to be the active oxidant. How and why it activates melatonin in the human body and what the origin of the product distributions is, are unknown. This encouraged us to do a detailed computational study on a typical human P450 isozyme, namely CYP1A1. We initially did a series of molecular dynamics simulations with substrate docked into several orientations. These simulations reveal a number of stable substrate-bound positions in the active site, which may lead to differences in substrate activation channels. Using tunneling analysis on the full protein structures, we show that two of the four binding conformations lead to open substrate-binding pockets. As a result, in these open pockets, the substrate is not tightly bound and can escape back into the solution. In the closed conformations, in contrast, the substrate is mainly oriented with the methoxy group pointing toward the heme, although under a different angle. We then created large quantum cluster models of the enzyme and focused on the chemical reaction mechanisms for melatonin activation, leading to competitive O-demethylation and C⁶-aromatic hydroxylation pathways. The calculations show that active site positioning determines the product distributions, but the bond that is activated is not necessarily closest to the heme in the enzyme-substrate complex. As such, the docking and molecular dynamics positioning of the substrate versus oxidant can give misleading predictions on product distributions. In particular, in quantum mechanics cluster model I, we observe that through a tight hydrogen bonding network, a preferential 6-hydroxylation of melatonin is obtained. However, O-demethylation becomes possible in alternative substrate-binding orientations that have the C⁶-aromatic ring position shielded. Finally, we investigated enzymatic and non-enzymatic O-demethylation processes and show that the hydrogen bonding network in the substrate-binding pocket can assist and perform this step prior to product release from the enzyme.

Origin of metabolites diversity and selectivity of P450 catalyzed benzo[a]pyrene metabolic activation

Polycyclic Aromatic Hydrocarbon (PAHs) presents one of the most abundant class of environmental pollutants. Recent study shows a lab-synthesized PAHs derivative, helicenium, can selectively kill cancer cells rather than normal cells, calling for the in-depth understanding of the metabolic process. However, the origin of metabolites diversity and selectivity of P450 catalyzed PAHs metabolic activation is still unclear to a great extent. Here we systematically investigated P450 enzymes catalyzed activation mechanism of a representative PAHs, benzo[a]pyrene (BaP), and found the corresponding activation process mainly involves two elementary steps: electrophilic addition and epoxidation. Electrophilic addition step is evidenced to be rate determining step. Two representative binding modes of BaP with P450 were found, which enables the electrophilic addition of Heme (FeO) to almost all the carbons of BaP. This electrophilic addition was proposed to be accelerated by the P450 enzyme environment when compared with the gas phase and water solvent. To dig deeper on the origin of metabolites diversity, we built several linear regression models to explore the structural-energy relationships. The selectivity was eventually attributed to the integrated effects of structural (e.g. O-C distance and O-C-Fe angle) and electrostatic parameters (e.g. charge of C and O) from both BaP and P450.

Using Physical Organic Chemistry Knowledge to Predict Unusual Metabolites of Synthetic Phenolic Antioxidants by Cytochrome P450

Biotransformation, especially by human CYP450 enzymes, plays a crucial role in regulating the toxicity of organic compounds in organisms, but is poorly understood for most emerging pollutants, as their numerous "unusual" biotransformation reactions cannot retrieve examples from the textbooks. Therefore, in order to predict the unknown metabolites with altering toxicological profiles, there is a realistic need to develop efficient methods to reveal the "unusual" metabolic mechanism of emerging pollutants. Combining experimental work with computational predictions has been widely accepted as an effective approach in studying complex metabolic reactions; however, the full quantum chemical computations may not be easily accessible for most environmentalists. Alternatively, this work practiced using the concepts from physical organic chemistry for studying the interrelationships between structure and reactivity of organic molecules, to reveal the "unusual" metabolic mechanism of synthetic phenolic antioxidants catalyzed by CYP450, for which the simple pencil-and-paper and property-computation methods based on physical organic chemistry were performed. The phenol-coupling product of butylated hydroxyanisole (BHA) (based on spin aromatic delocalization) and ipso-addition quinol metabolite of butylated hydroxytoluene (BHT) (based on hyperconjugative effect) were predicted as two "unusual" metabolites, which were further confirmed by our in vitro analysis. We hope this easily handled approach will promote environmentalists to attach importance to physical organic chemistry, with an eye to being able to use the knowledge gained to efficiently predict the fates of substantial unknown synthesized organic compounds in the future.

Determining the inherent selectivity for carbon radical hydroxylation versus halogenation with high-spin oxoiron(iv)-halide complexes: a concerted rebound step

A synthetic iron model can process both halogenation and hydroxylation with vague selectivity, which is different from halogenase even though these structures are used for the simulation of halogenase. The key factor of the synthetic oxo-iron model mediated hydroxylation or the halogenation is still under debate. Herein density functional theory calculation is used to investigate the hydroxylation versus halogenation of propylene by the complex [Fe^IV(O)(TQA)(X)]⁺ (X = F, Cl, Br). Our results suggest that a concerted rebound mechanism (between the -X and the hydroxyl ligands after the hydrogen abstraction) leads to the formation of two different kinds of products.

DFT calculation for the hydroxylation versus halogenation of propylene by [Fe^IV(O)(TQA)X]⁺ (X = F, Cl and Br) reveals that after hydrogen abstraction, halogen and oxygen rebound reactions are a synergistic process.

Development of Next-Generation Probiotics by Investigating the Interrelationships between Gastrointestinal Microbiota and Diarrhea in Preruminant Holstein Calves

Simple Summary

The present study investigated the relationship between gastrointestinal microbiota and diarrhea in preruminant calves by using immune-related markers and further isolating specific bacterial strains, enriched in clinically healthy individuals, for potential next-generation probiotics. The gathering of microbiomic data strongly indicated the possible beneficial effects of Bifidobacterium longum subsp. longum. With further screening and isolating with immunomodulatory and antagonistic effects, two Bifidobacterium longum subsp. longum strains might be expected to emerge as next-generation probiotics. The finding here might provide a solution for preventing gastrointestinal disorders for preruminant calves without sustained periods of administration through inhibiting the infectious bacteria, immunomodulatory effect and possible modulating microbiota.

Abstract

(1) Background: We aimed to isolate and identify potential next-generation probiotics (NGP) by investigating the interrelationships between gastrointestinal microbiota and diarrhea in preruminant Holstein calves. (2) Material and methods: Twenty preruminant Holstein calves were divided into healthy and diarrheic groups after the combination outcomes of veterinary diagnosis and fecal scores. The fecal microbiome, plasma cytokines, plasma immunoglobulin (Ig) G and haptoglobin were analyzed. The potential probiotic bacteria were identified by comparing the microbiota difference between healthy and diarrheic calves and correlation analysis with fecal scores and inflammatory markers. The identified bacteria were also isolated for further evaluation for antimicrobial activities and immunoregulatory effects. (3) Results: Microbiota analysis suggested that Ruminococcaceae_UCG_014, Bifidobacterium and Pseudoflavonifractor positively correlated with bovine IgG and negatively correlated with fecal score; inflammatory factors, bovine HP, and IL-8 were classified as beneficial bacteria contributing to the health of the calves. The alternation of gut microbial composition also induced changes in the functional gene enrichment of gut microbiota in calves. The gathering of microbiomic data strongly indicated the possible beneficial effects of Bifidobacterium longum subsp. longum, expected to develop as NGP. After isolation and evaluation of the potential functionality in vitro, two specific bifidobacterial strains demonstrated antimicrobial activities and immunoregulatory effects. (4) Conclusions: The results provide a new probiotic searching approach for preventing gastrointestinal disorders in preruminant calves. Further animal study is necessary to verify the results.

Computational Investigation of the Bisphenolic Drug Metabolism by Cytochrome P450: What Factors Favor Intramolecular Phenol Coupling

Intramolecular phenol coupling reactions of alkaloids can lead to active metabolites catalyzed by the mammalian cytochrome P450 enzyme (P450); however, the mechanistic knowledge of such an "unusual" process is lacking. This work performs density functional theory computations to reveal the P450-mediated metabolic pathway leading from R-reticuline to the morphine precursor salutaridine by exploring possible intramolecular phenol coupling mechanisms involving diradical coupling, radical addition, and electron transfer. The computed results show that the outer-sphere electron transfer with a high barrier (>20.0 kcal/mol) is unlikely to happen. However, for inter-sphere intramolecular phenol coupling, it reveals that intramolecular phenol coupling of R-reticuline proceeds via the diradical mechanism consecutively by compound I and protonated compound II of P450 rather than the radical addition mechanism. The existence of a much higher radical rebound barrier than that of H-abstraction in the quartet high-spin state can endow the R-reticuline phenoxy radical with a sufficient lifetime to enable intramolecular phenol coupling, while the H-abstraction/radical rebound mode with a negligible rebound barrier leading to phenol hydroxylation can only happen in the doublet low-spin state. Therefore, the ratio [coupling]/[hydroxylation] can be approximately reflected by the relative yield of the high-spin and low-spin H-abstraction by P450, which thus can provide a theoretical ratio of 16:1 for R-reticuline, which is in accordance with previous experimental results. Especially, the high rebound barrier of the phenoxy radical derived from the weak electron-donating ability of the phenoxy radical is revealed as an intrinsic nature. Therefore, the revealed intramolecular phenol coupling mechanism can be potentially extended to several other bisphenolic drugs to infer groups of unexpected metabolites in organisms.

Precision Biotransformation of Emerging Pollutants by Human Cytochrome P450 Using Computational-Experimental Synergy: A Case Study of Tris(1,3-dichloro-2-propyl) Phosphate

Precision biotransformation is an envisioned strategy offering detailed insights into biotransformation pathways in real environmental settings using experimentally guided high-accuracy quantum chemistry. Emerging pollutants, whose metabolites are easily overlooked but may cause idiosyncratic toxicity, are important targets of such a strategy. We demonstrate here that complex metabolic reactions of tris(1,3-dichloro-2-propyl) phosphate (TDCIPP) catalyzed by human CYP450 enzymes can be mapped via a three-step synergy strategy: (i) screening the possible metabolites via high-throughout (moderate-accuracy) computations; (ii) analyzing the proposed metabolites in vitro by human liver microsomes and recombinant human CYP450 enzymes; and (iii) rationalizing the experimental data via precise mechanisms using high-level targeted computations. Through the bilateral dialogues from qualitative to semi-quantitative to quantitative levels, we show how TDCIPP metabolism especially by CYP3A4 generates bis(1,3-dichloro-2-propyl) phosphate (BDCIPP) as an O-dealkylation metabolite and bis(1,3-dichloro-2-propyl) 3-chloro-1-hydroxy-2-propyl phosphate (alcohol_β-dehalogen) as a dehalogenation/reduction metabolite via the initial rate-determining H-abstraction from αC- and βC-positions. The relative yield ratio [dehalogenation/reduction]/[O-dealkylation] is derived from the relative barriers of H-abstraction at the βC- and αC-positions by CYP3A4, estimated as 0.002 to 0.23, viz., an in vitro measured ratio of 0.04. Importantly, alcohol_β-dehalogen formation points to a new mechanism involving successive oxidation and reduction functions of CYP450, with its precursor aldehyde_β-dehalogen being a key intermediate detected by trapping assays and rationalized by computations. We conclude that the proposed three-step synergy strategy may meet the increasing challenge of elucidating biotransformation mechanisms of substantial synthesized organic compounds in the future.

The stereoselective oxidation of para-substituted benzenes by a cytochrome P450 biocatalyst

The serine 244 to aspartate (S244D) variant of the cytochrome P450 enzyme CYP199A4 was used to expand its substrate range beyond benzoic acids. Substrates, in which the carboxylate group of the benzoic acid moiety is replaced, were oxidised with high activity by the S244D mutant (product formation rates > 60 nmol.(nmol-CYP) -1 .min -1 ) and with total turnover numbers of up to 20,000. Ethyl α-hydroxylation was more rapid than methyl oxidation, styrene epoxidation and S -oxidation. The S244D mutant catalysed the ethyl hydroxylation, epoxidation and sulfoxidation reactions with an excess of one stereoisomer (in some instances up to >98%). The crystal structure of 4-methoxybenzoic acid-bound CYP199A4 S244D showed that the active site architecture and the substrate orientation were similar to that of the WT enzyme. Overall, this work demonstrates that CYP199A4 can catalyse the stereospecific hydroxylation, epoxidation or sulfoxidation of substituted benzene substrates under mild conditions resulting in more sustainable transformations using this heme monooxygenase enzyme.

Theoretical Investigation on the Elusive Reaction Mechanism of Spirooxindole Formation Mediated by Cytochrome P450s: A Nascent Feasible Charge-Shift C-O Bond Makes a Difference

Spirooxindoles are pivotal biofunctional groups widely distributed in natural products and clinic drugs. However, construction of such subtle chiral skeletons is a long-standing challenge to both organic and bioengineering scientists. The knowledge of enzymatic spirooxindole formation in nature may inspire rational design of new catalysts. To this end, we presented a theoretical investigation on the elusive mechanism of the spiro-ring formation at the 3-position of oxindole mediated by cytochrome P450 enzymes (P450). Our calculated results demonstrated that the electrophilic attack of CpdI, the active species of P450, to the substrate, shows regioselectivity, i.e., the attack at the C9 position forms a tetrahedral intermediate involving an unusual feasible charge-shift C9^δ+-O^δ- bond, while the attack at the C1 position forms an epoxide intermediate. The predominant route is the first route with the charge-shift bonding intermediate due to holding a relatively lower barrier by >5 kcal mol^-1 than the epoxide route, which fits the experimental observations. Such a delocalized charge-shift bond facilitates the formation of a spiro-ring mainly through elongation of the C1-C9 bond to eliminate the aromatization of the tricyclic beta-carboline. Our theoretical results shed profound mechanistic insights for the first time into the elusive spirooxindole formation mediated by P450s.

Product Distributions of Cytochrome P450 OleTJE with Phenyl-Substituted Fatty Acids: A Computational Study

There are two types of cytochrome P450 enzymes in nature, namely, the monooxygenases and the peroxygenases. Both enzyme classes participate in substrate biodegradation or biosynthesis reactions in nature, but the P450 monooxygenases use dioxygen, while the peroxygenases take H2O2 in their catalytic cycle instead. By contrast to the P450 monooxygenases, the P450 peroxygenases do not require an external redox partner to deliver electrons during the catalytic cycle, and also no external proton source is needed. Therefore, they are fully self-sufficient, which affords them opportunities in biotechnological applications. One specific P450 peroxygenase, namely, P450 OleTJE, reacts with long-chain linear fatty acids through oxidative decarboxylation to form hydrocarbons and, as such, has been implicated as a suitable source for the biosynthesis of biofuels. Unfortunately, the reactions were shown to produce a considerable amount of side products originating from Cα and Cβ hydroxylation and desaturation. These product distributions were found to be strongly dependent on whether the substrate had substituents on the Cα and/or Cβ atoms. To understand the bifurcation pathways of substrate activation by P450 OleTJE leading to decarboxylation, Cα hydroxylation, Cβ hydroxylation and Cα-Cβ desaturation, we performed a computational study using 3-phenylpropionate and 2-phenylbutyrate as substrates. We set up large cluster models containing the heme, the substrate and the key features of the substrate binding pocket and calculated (using density functional theory) the pathways leading to the four possible products. This work predicts that the two substrates will react with different reaction rates due to accessibility differences of the substrates to the active oxidant, and, as a consequence, these two substrates will also generate different products. This work explains how the substrate binding pocket of P450 OleTJE guides a reaction to a chemoselectivity.

Density Functional Theory Study into the Reaction Mechanism of Isonitrile Biosynthesis by the Nonheme Iron Enzyme ScoE

The nonheme iron enzyme ScoE catalyzes the biosynthesis of an isonitrile substituent in a peptide chain. To understand details of the reaction mechanism we created a large active site cluster model of 212 atoms that contains substrate, the active oxidant and the first- and second-coordination sphere of the protein and solvent. Several possible reaction mechanisms were tested and it is shown that isonitrile can only be formed through two consecutive catalytic cycles that both use one molecule of dioxygen and α-ketoglutarate. In both cycles the active species is an iron(IV)-oxo species that in the first reaction cycle reacts through two consecutive hydrogen atom abstraction steps: first from the N–H group and thereafter from the C–H group to desaturate the NH-CH2 bond. The alternative ordering of hydrogen atom abstraction steps was also tested but found to be higher in energy. Moreover, the electronic configurations along that pathway implicate an initial hydride transfer followed by proton transfer. We highlight an active site Lys residue that is shown to donate charge in the transition states and influences the relative barrier heights and bifurcation pathways. A second catalytic cycle of the reaction of iron(IV)-oxo with desaturated substrate starts with hydrogen atom abstraction followed by decarboxylation to give isonitrile directly. The catalytic cycle is completed with a proton transfer to iron(II)-hydroxo to generate the iron(II)-water resting state. The work is compared with experimental observation and previous computational studies on this system and put in a larger perspective of nonheme iron chemistry.

A comprehensive insight into aldehyde deformylation: mechanistic implications from biology and chemistry

Aldehyde deformylation is an important reaction in biology, organic chemistry and inorganic chemistry and the process has been widely applied and utilized. For instance, in biology, the aldehyde deformylation reaction has wide differences in biological function, whereby cyanobacteria convert aldehydes into alkanes or alkenes, which are used as natural products for, e.g., defense mechanisms. By contrast, the cytochromes P450 catalyse the biosynthesis of hormones, such as estrogen, through an aldehyde deformylation reaction step. In organic chemistry, the aldehyde deformylation reaction is a common process for replacing functional groups on a molecule, and as such, many different synthetic methods and procedures have been reported that involve an aldehyde deformylation step. In bioinorganic chemistry, a variety of metal(iii)-peroxo complexes have been synthesized as biomimetic models and shown to react efficiently with aldehydes through deformylation reactions. This review paper provides an overview of the various aldehyde deformylation reactions in organic chemistry, biology and biomimetic model systems, and shows a broad range of different chemical reaction mechanisms for this process. Although a nucleophilic attack at the carbonyl centre is the consensus reaction mechanism, several examples of an alternative electrophilic reaction mechanism starting with hydrogen atom abstraction have been reported as well. There is still much to learn and to discover on aldehyde deformylation reactions, as deciphered in this review paper.

Catalytic approach to in vivo metabolism of atractylenolide III using biomimetic iron–porphyrin complexes

In vivo oxygenation of atractylenolide III during the metabolism and in situ oxidation mechanism by an iron–porphyrin iron(iv)-oxo porphyrin π-cation–radical complex.

Regioselectivity of hyoscyamine 6β-hydroxylase-catalysed hydroxylation as revealed by high-resolution structural information and QM/MM calculations

Hyoscyamine 6β-hydroxylase (H6H) is a bifunctional non-heme 2-oxoglutarate/Fe2+-dependent dioxygenase that catalyzes the two final steps in the biosynthesis of scopolamine. Based on high resolution crystal structures of H6H from Datura metel, detailed information on substrate binding was obtained that provided insights into the onset of the enzymatic process. In particular, the role of two prominent residues was revealed - Glu-116 that interacts with the tertiary amine located on the hyoscyamine tropane moiety and Tyr-326 that forms CH-π hydrogen bonds with the hyoscyamine phenyl ring. The structures were used as the basis for QM/MM calculations that provided an explanation for the regioselectivity of the hydroxylation reaction on the hyoscyamine tropane moiety (C6 vs. C7) and quantified contributions of active site residues to respective barrier heights.

Computational studies of DNA base repair mechanisms by nonheme iron dioxygenases: selective epoxidation and hydroxylation pathways

DNA base repair mechanisms of alkylated DNA bases is an important reaction in chemical biology and particularly in the human body. It is typically catalyzed by an α-ketoglutarate-dependent nonheme iron dioxygenase named the AlkB repair enzyme. In this work we report a detailed computational study into the structure and reactivity of AlkB repair enzymes with alkylated DNA bases. In particular, we investigate the aliphatic hydroxylation and C[double bond, length as m-dash]C epoxidation mechanisms of alkylated DNA bases by a high-valent iron(iv)-oxo intermediate. Our computational studies use quantum mechanics/molecular mechanics methods on full enzymatic structures as well as cluster models on active site systems. The work shows that the iron(iv)-oxo species is rapidly formed after dioxygen binding to an iron(ii) center and passes a bicyclic ring structure as intermediate. Subsequent cluster models explore the mechanism of substrate hydroxylation and epoxidation of alkylated DNA bases. The work shows low energy barriers for substrate activation and consequently energetically feasible pathways are predicted. Overall, the work shows that a high-valent iron(iv)-oxo species can efficiently dealkylate alkylated DNA bases and return them into their original form.

Deciphering the origin of million-fold reactivity observed for the open core diiron [HO-FeIII-O-FeIV[double bond, length as m-dash]O]2+ species towards C-H bond activation: role of spin-states, spin-coupling, and spin-cooperation

High-valent metal-oxo species have been characterised as key intermediates in both heme and non-heme enzymes that are found to perform efficient aliphatic hydroxylation, epoxidation, halogenation, and dehydrogenation reactions. Several biomimetic model complexes have been synthesised over the years to mimic both the structure and function of metalloenzymes. The diamond-core [Fe2(μ-O)2] is one of the celebrated models in this context as this has been proposed as the catalytically active species in soluble methane monooxygenase enzymes (sMMO), which perform the challenging chemical conversion of methane to methanol at ease. In this context, a report of open core [HO(L)FeIII-O-FeIV(O)(L)]2+ (1) gains attention as this activates C-H bonds a million-fold faster compared to the diamond-core structure and has the dual catalytic ability to perform hydroxylation as well as desaturation with organic substrates. In this study, we have employed density functional methods to probe the origin of the very high reactivity observed for this complex and also to shed light on how this complex performs efficient hydroxylation and desaturation of alkanes. By modelling fifteen possible spin-states for 1 that could potentially participate in the reaction mechanism, our calculations reveal a doublet ground state for 1 arising from antiferromagnetic coupling between the quartet FeIV centre and the sextet FeIII centre, which regulates the reactivity of this species. The unusual stabilisation of the high-spin ground state for FeIV[double bond, length as m-dash]O is due to the strong overlap of with the orbital, reducing the antibonding interactions via spin-cooperation. The electronic structure features computed for 1 are consistent with experiments offering confidence in the methodology chosen. Further, we have probed various mechanistic pathways for the C-H bond activation as well as -OH rebound/desaturation of alkanes. An extremely small barrier height computed for the first hydrogen atom abstraction by the terminal FeIV[double bond, length as m-dash]O unit was found to be responsible for the million-fold activation observed in the experiments. The barrier height computed for -OH rebound by the FeIII-OH unit is also smaller suggesting a facile hydroxylation of organic substrates by 1. A strong spin-cooperation between the two iron centres also reduces the barrier for second hydrogen atom abstraction, thus making the desaturation pathway competitive. Both the spin-state as well as spin-coupling between the two metal centres play a crucial role in dictating the reactivity for species 1. By exploring various mechanistic pathways, our study unveils the fact that the bridged μ-oxo group is a poor electrophile for both C-H activation as well for -OH rebound. As more and more evidence is gathered in recent years for the open core geometry of sMMO enzymes, the idea of enhancing the reactivity via an open-core motif has far-reaching consequences.

Hydroxyl Transfer to Carbon Radicals by Mn(OH) vs Fe(OH) Corrole Complexes

The transfer of •OH from metal-hydroxo species to carbon radicals (R•) to give hydroxylated products (ROH) is a fundamental process in metal-mediated heme and nonheme C-H bond oxidations. This step, often referred to as the hydroxyl "rebound" step, is typically very fast, making direct study of this process challenging if not impossible. In this report, we describe the reactions of the synthetic models M(OH)(ttppc) (M = Fe (1), Mn (3); ttppc = 5,10,15-tris(2,4,6-triphenyl)phenyl corrolato^3-) with a series of triphenylmethyl carbon radical (R•) derivatives ((4-X-C₆H₄)₃C•; X = OMe, tBu, Ph, Cl, CN) to give the one-electron reduced M^III(ttppc) complexes and ROH products. Rate constants for 3 for the different radicals ranged from 11.4(1) to 58.4(2) M^-1 s^-1, as compared to those for 1, which fall between 0.74(2) and 357(4) M^-1 s^-1. Linear correlations for Hammett and Marcus plots for both Mn and Fe were observed, and the small magnitudes of the slopes for both correlations imply a concerted •OH transfer reaction for both metals. Eyring analyses of reactions for 1 and 3 with (4-X-C₆H₄)₃C• (X = tBu, CN) also give good linear correlations, and a comparison of the resulting activation parameters highlight the importance of entropy in these •OH transfer reactions. Density functional theory calculations of the reaction profiles show a concerted process with one transition state for all radicals investigated and help to explain the electronic features of the OH rebound process.

Lignin Biodegradation by a Cytochrome P450 Enzyme: A Computational Study into Syringol Activation by GcoA

A recently characterized cytochrome P450 isozyme GcoA activates lignin components through a selective O-demethylation or alternatively an acetal formation reaction. These are important reactions in biotechnology and, because lignin is readily available; it being the main component in plant cell walls. In this work we present a density functional theory study on a large active site model of GcoA to investigate syringol activation by an iron(IV)-oxo heme cation radical oxidant (Compound I) leading to hemiacetal and acetal products. Several substrate-binding positions were tested and full energy landscapes calculated. The study shows that substrate positioning determines the product distributions. Thus, with the phenol group pointing away from the heme, an O-demethylation is predicted, whereas an initial hydrogen-atom abstraction of the weak phenolic O-H group would trigger a pathway leading to ring-closure to form acetal products. Predictions on how to engineer P450 GcoA to get more selective product distributions are given.

Catalytic Mechanism of Aromatic Nitration by Cytochrome P450 TxtE: Involvement of a Ferric-Peroxynitrite Intermediate

The cytochromes P450 are heme-dependent enzymes that catalyze many vital reaction processes in the human body related to biodegradation and biosynthesis. They typically act as mono-oxygenases; however, the recently discovered P450 subfamily TxtE utilizes O₂ and NO to nitrate aromatic substrates such as L-tryptophan. A direct and selective aromatic nitration reaction may be useful in biotechnology for the synthesis of drugs or small molecules. Details of the catalytic mechanism are unknown, and it has been suggested that the reaction should proceed through either an iron(III)-superoxo or an iron(II)-nitrosyl intermediate. To resolve this controversy, we used stopped-flow kinetics to provide evidence for a catalytic cycle where dioxygen binds prior to NO to generate an active iron(III)-peroxynitrite species that is able to nitrate l-Trp efficiently. We show that the rate of binding of O₂ is faster than that of NO and also leads to l-Trp nitration, while little evidence of product formation is observed from the iron(II)-nitrosyl complex. To support the experimental studies, we performed density functional theory studies on large active site cluster models. The studies suggest a mechanism involving an iron(III)-peroxynitrite that splits homolytically to form an iron(IV)-oxo heme (Compound II) and a free NO₂ radical via a small free energy of activation. The latter activates the substrate on the aromatic ring, while compound II picks up the ipso-hydrogen to form the product. The calculations give small reaction barriers for most steps in the catalytic cycle and, therefore, predict fast product formation from the iron(III)-peroxynitrite complex. These findings provide the first detailed insight into the mechanism of nitration by a member of the TxtE subfamily and highlight how the enzyme facilitates this novel reaction chemistry.

Elucidating the effect of the ionic liquid type and alkyl chain length on the stability of ionic liquid-iron porphyrin complexes

The present study is motivated by the long-term objective of understanding how ionic liquids are biodegraded by cytochrome P450, which contains iron porphyrin (FeP) serving as the catalytic center. To this end, the current study is designed to elucidate the impact of types and conformations of ionic liquids on the binding energy with FeP, the key interactions that stabilize the ionic liquid-FeP complex, and how the electron uptake ability of FeP is altered in the presence of ionic liquids. Four classes of ionic liquids are considered: 1-alkyl-3-methylimidazolium, 1-alkyl-pyridinium, 1-alkylsulfonium, and N-methyl-N-alkylpyrrolidinium. The influence of linear alkyl chains of ethyl, butyl, hexyl, octyl, and decyl is examined on the favorable binding modes with FeP, considering two widely different conformations: tail up and tail down with respect to FeP. Electronic structure calculations are performed at the M06 level of theory with the 6-31G(d,p) basis set for C, H, and N atoms, while the Lanl2DZ basis set is employed for Fe. Donor-acceptor interactions contributing to the binding of ionic liquids to FeP are unraveled through the natural bond orbital analysis. The results from this study indicate that the binding energies are dependent not only on the class of ionic liquids but also on the conformations presented to FeP. The propensity of FeP to acquire an electron is significantly enhanced in the presence of ionic liquid cations, irrespective of the type and the alkyl chain length.

Pharmacokinetics and Metabolism of Liposome-Encapsulated 2,4,6-Trihydroxygeranylacetophenone in Rats Using High-Resolution Orbitrap Liquid Chromatography Mass Spectrometry

2,4,6-trihydroxy-3-geranylacetophenone (tHGA) is a bioactive compound that shows excellent anti-inflammatory properties. However, its pharmacokinetics and metabolism have yet to be evaluated. In this study, a sensitive LC-HRMS method was developed and validated to quantify tHGA in rat plasma. The method showed good linearity (0.5–80 ng/mL). The accuracy and precision were within 10%. Pharmacokinetic investigations were performed on three groups of six rats. The first two groups were given oral administrations of unformulated and liposome-encapsulated tHGA, respectively, while the third group received intraperitoneal administration of liposome-encapsulated tHGA. The maximum concentration (C_max), the time required to reach C_max (t_max), elimination half-life (t_1/2) and area under curve (AUC_0–24) values for intraperitoneal administration were 54.6 ng/mL, 1.5 h, 6.7 h, and 193.9 ng/mL·h, respectively. For the oral administration of unformulated and formulated tHGA, C_max values were 5.4 and 14.5 ng/mL, t_max values were 0.25 h for both, t_1/2 values were 6.9 and 6.6 h, and AUC_0–24 values were 17.6 and 40.7 ng/mL·h, respectively. The liposomal formulation improved the relative oral bioavailability of tHGA from 9.1% to 21.0% which was a 2.3-fold increment. Further, a total of 12 metabolites were detected and structurally characterized. The metabolites were mainly products of oxidation and glucuronide conjugation.

Bioengineering of Cytochrome P450 OleTJE: How Does Substrate Positioning Affect the Product Distributions?

The cytochromes P450 are versatile enzymes found in all forms of life. Most P450s use dioxygen on a heme center to activate substrates, but one class of P450s utilizes hydrogen peroxide instead. Within the class of P450 peroxygenases, the P450 OleT_JE isozyme binds fatty acid substrates and converts them into a range of products through the α-hydroxylation, β-hydroxylation and decarboxylation of the substrate. The latter produces hydrocarbon products and hence can be used as biofuels. The origin of these product distributions is unclear, and, as such, we decided to investigate substrate positioning in the active site and find out what the effect is on the chemoselectivity of the reaction. In this work we present a detailed computational study on the wild-type and engineered structures of P450 OleT_JE using a combination of density functional theory and quantum mechanics/molecular mechanics methods. We initially explore the wild-type structure with a variety of methods and models and show that various substrate activation transition states are close in energy and hence small perturbations as through the protein may affect product distributions. We then engineered the protein by generating an in silico model of the double mutant Asn242Arg/Arg245Asn that moves the position of an active site Arg residue in the substrate-binding pocket that is known to form a salt-bridge with the substrate. The substrate activation by the iron(IV)-oxo heme cation radical species (Compound I) was again studied using quantum mechanics/molecular mechanics (QM/MM) methods. Dramatic differences in reactivity patterns, barrier heights and structure are seen, which shows the importance of correct substrate positioning in the protein and the effect of the second-coordination sphere on the selectivity and activity of enzymes.

Computational Biotransformation Profile of Emerging Phenolic Pollutants by Cytochromes P450: Phenol-Coupling Mechanism

Phenols are ubiquitous environmental pollutants, whose biotransformation involving phenol coupling catalyzed by cytochromes P450 may produce more lipophilic and toxic metabolites. Density functional theory (DFT) computations were performed to explore the debated phenol-coupling mechanisms, taking triclosan as a model substrate. We find that a diradical pathway facilitated by compound I and protonated compound II of P450 is favored vs alternative radical addition or electron-transfer mechanisms. The identified diradical coupling resembles a "two-state reactivity" from compound I characterized by significantly high rebound barriers of the phenoxy radicals, which can be formulated into three equations for calculating the ratio [coupling]/[hydroxylation]. A higher barrier for rebound than for H-abstraction in high-spin triclosan can facilitate the phenoxy radical dissociation and thus enable phenol coupling, while H-abstraction/radical rebound causing phenol hydroxylation via minor rebound barriers mostly occurs via the low-spin state. Therefore, oxidation of triclosan by P450 fits the first equation with a ratio [coupling]/[hydroxylation] of 1:4, consistent with experimental data indicating different extents of triclosan coupling (6-40%). The high rebound barrier of phenoxy radicals, as a key for the mechanistic identification of phenol coupling vs hydroxylation, originates from their weak electron donor ability due to spin aromatic delocalization. We envision that the revealed mechanism can be extended to the cross-coupling reactions between different phenolic pollutants, and the coupling reactions of several other aromatic pollutants, to infer unknown metabolites.

Sluggish reactivity by a nonheme iron(iv)-tosylimido complex as compared to its oxo analogue

A comparative spectroscopic and computational study of reactivity between ferryl-tosylimido and ferryl-oxo complexes of two biomimetic model systems. The Fe(iv)-tosylimido complex was found to be sluggish in comparison to its fellow oxo counterpart.

Selective Hydrogen Atom Abstraction from Dihydroflavonol by a Nonheme Iron Center Is the Key Step in the Enzymatic Flavonol Synthesis and Avoids Byproducts

The plant non-heme iron dioxygenase flavonol synthase performs a regioselective desaturation reaction as part of the biosynthesis of the signaling molecule flavonol that triggers the growing of leaves and flowers. These compounds also have health benefits for humans. Desaturation of aliphatic compounds generally proceeds through two consecutive hydrogen atom abstraction steps from two adjacent carbon atoms and in nature often is performed by a high-valent iron(IV)-oxo species. We show that the order of the hydrogen atom abstraction steps, however, is opposite of those expected from the C-H bond strengths in the substrate and determines the product distributions. As such, flavonol synthase follows a negative catalysis mechanism. Using density functional theory methods on large active-site model complexes, we investigated pathways for desaturation and hydroxylation by an iron(IV)-oxo active-site model. Contrary to thermochemical predictions, we find that the oxidant abstracts the hydrogen atom from the strong C²-H bond rather than the weaker C³-H bond of the substrate first. We analyze the origin of this unexpected selective hydrogen atom abstraction pathway and find that the alternative C³-H hydrogen atom abstraction would be followed by a low-energy and competitive substrate hydroxylation mechanism hence, should give considerable amount of byproducts. Our computational modeling studies show that substrate positioning in flavonol synthase is essential, as it guides the reactivity to a chemo- and regioselective substrate desaturation from the C²-H group, leading to desaturation products efficiently.

Small molecule IVQ, as a prodrug of gluconeogenesis inhibitor QVO, efficiently ameliorates glucose homeostasis in type 2 diabetic mice

Gluconeogenesis is a major source of hyperglycemia in patients with type 2 diabetes mellitus (T2DM), thus targeting gluconeogenesis to suppress glucose production is a promising strategy for anti-T2DM drug discovery. In our preliminary in vitro studies, we found that a small-molecule (E)-3-(2-(quinoline-4-yl)vinyl)-1H-indol-6-ol (QVO) inhibited the hepatic glucose production (HGP) in primary hepatocytes. We further revealed that QVO suppressed hepatic gluconeogenesis involving calmodulin-dependent protein kinase kinase β- and liver kinase B1-adenosine monophosphate-activated protein kinase (AMPK) pathways as well as AMPK-independent mitochondrial function-related signaling pathway. To evaluate QVO's anti-T2DM activity in vivo, which was impeded by the complicated synthesis route of QVO with a low yield, we designed and synthesized 4-[2-(1H-indol-3-yl)vinyl]quinoline (IVQ) as a prodrug with easier synthesis route and higher yield. IVQ did not inhibit the HGP in primary hepatocytes in vitro. Pharmacokinetic studies demonstrated that IVQ was quickly converted to QVO in mice and rats following administration. In both db/db and ob/ob mice, oral administration of IVQ hydrochloride (IVQ-HCl) (23 and 46 mg/kg every day, for 5 weeks) ameliorated hyperglycemia, and suppressed hepatic gluconeogenesis and activated AMPK signaling pathway in the liver tissues. Furthermore, IVQ caused neither cardiovascular system dysfunction nor genotoxicity. The good druggability of IVQ has highlighted its potential in the treatment of T2DM and the prodrug design for anti-T2DM drug development.

DFT study on the effect of proximal residues on the Mycobacterium tuberculosis catalase-peroxidase (katG) heme compound I intermediate and its bonding interaction with isoniazid

Isoniazid (INH) is converted into isonicotinyl radical through its interaction with the catalase-peroxidase (katG) enzyme present in the cells of Mycobacterium tuberculosis (M. tb.), the bacteria that causes the tuberculosis disease. This process is important because resistance of M. tb. cells to INH treatment has been associated with the failure of completion of this process. However, this process is poorly understood and there are a variety of conflicting theories about the details of the mechanism of this process. One theory suggests that INH binds to katG and transfers a single electron to the heme while the heme is in its two electron oxidized state, compound I [Fe(iv)Por˙+] (CpdI). In this study, DFT calculations at the UB3LYP/6-31g(d)-LANL2DZ level of theory are used to study the M. tb. katG CpdI molecule. Different models of the M. tb. CpdI molecule were prepared and the calculations revealed the impact of Trp321 on the electronic properties of the heme. Without Trp321 the heme assumed a triradical state with single electrons on the πxy and πyz orbitals of Fe and another on the a2u orbital of the porphyrin ring that can either be coupled with the first two, to assume a quartet state, or decoupled to form a doublet state. With Trp321, however, a transfer of an electron from the πTrp orbital to a2u porphyrin orbital leads to loss of radical character of the porphyrin and leaves the Trp321 group with a radical character. INH was observed to have strong interaction with CpdI and the absence of Trp321 significantly decreased the binding energy by 2 kcal mol-1 explaining the importance of Trp321 in the binding of INH. The results in this study show the importance of Trp321 in the binding of INH and its effect on its electronic properties, which is consistent with previous experimental findings that mutation of Trp321 results in an increase in drug resistance.

Bridging the Homogeneous-Heterogeneous Divide: Modeling Spin for Reactivity in Single Atom Catalysis

Single atom catalysts (SACs) are emergent catalytic materials that have the promise of merging the scalability of heterogeneous catalysts with the high activity and atom economy of homogeneous catalysts. Computational, first-principles modeling can provide essential insight into SAC mechanism and active site configuration, where the sub-nm-scale environment can challenge even the highest-resolution experimental spectroscopic techniques. Nevertheless, the very properties that make SACs attractive in catalysis, such as localized d electrons of the isolated transition metal center, make them challenging to study with conventional computational modeling using density functional theory (DFT). For example, Fe/N-doped graphitic SACs have exhibited spin-state dependent reactivity that remains poorly understood. However, spin-state ordering in DFT is very sensitive to the nature of the functional approximation chosen. In this work, we develop accurate benchmarks from correlated wavefunction theory (WFT) for relevant octahedral complexes. We use those benchmarks to evaluate optimal DFT functional choice for predicting spin state ordering in small octahedral complexes as well as models of pyridinic and pyrrolic nitrogen environments expected in larger SACs. Using these guidelines, we determine Fe/N-doped graphene SAC model properties and reactivity as well as their sensitivities to DFT functional choice. Finally, we conclude with broad recommendations for computational modeling of open-shell transition metal single-atom catalysts.