Recent density functional theory model calculations of drug metabolism by cytochrome P450

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.

Insights on cyclophosphamide metabolism and anticancer mechanism of action: A computational study

The oxazaphosphorine cyclophosphamide (CP) is a DNA-alkylating agent commonly used in cancer chemotherapy. This anticancer agent is administered as a prodrug activated by a liver cytochrome P450-catalyzed 4-hydroxylation reaction that yields the active, cytotoxic metabolite. The primary metabolite, 4-hydroxycyclophosphamide, equilibrates with the ring-open aldophosphamide that undergoes β-elimination to yield the therapeutically active DNA cross-linking phosphoramide mustard and the byproduct acrolein. The present paper presents a DFT investigation of the different metabolic phases and an insight into the mechanism by which CP exerts its cytotoxic action. A detailed computational analysis of the energy profiles describing all the involved transformations and the mechanism of DNA alkylation is given with the aim to contribute to an increase of knowledge that, after more than 60 years of unsuccessful attempts, can lead to the design and development of a new generation of oxazaphosphorines.

Predictive Modelling in pharmacokinetics: from in-silico simulations to personalized medicine

INTRODUCTION

Pharmacokinetic parameters assessment is a critical aspect of drug discovery and development, yet challenges persist due to limited training data. Despite advancements in machine learning and in-silico predictions, scarcity of data hampers accurate prediction of drug candidates' pharmacokinetic properties.

AREAS COVERED

The study highlights current developments in human pharmacokinetic prediction, talks about attempts to apply synthetic approaches for molecular design, and searches several databases, including Scopus, PubMed, Web of Science, and Google Scholar. The article stresses importance of rigorous analysis of machine learning model performance in assessing progress and explores molecular modeling (MM) techniques, descriptors, and mathematical approaches. Transitioning to clinical drug development, article highlights AI (Artificial Intelligence) based computer models optimizing trial design, patient selection, dosing strategies, and biomarker identification. In-silico models, including molecular interactomes and virtual patients, predict drug performance across diverse profiles, underlining the need to align model results with clinical studies for reliability. Specialized training for human specialists in navigating predictive models is deemed critical. Pharmacogenomics, integral to personalized medicine, utilizes predictive modeling to anticipate patient responses, contributing to more efficient healthcare system. Challenges in realizing potential of predictive modeling, including ethical considerations and data privacy concerns, are acknowledged.

EXPERT OPINION

AI models are crucial in drug development, optimizing trials, patient selection, dosing, and biomarker identification and hold promise for streamlining clinical investigations.

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.

Phylogenomic Analysis of Cytochrome P450 Gene Superfamily and Their Association with Flavonoids Biosynthesis in Peanut (Arachis hypogaea L.)

Cytochrome P450s (CYPs) constitute extensive enzyme superfamilies in the plants, playing pivotal roles in a multitude of biosynthetic and detoxification pathways essential for growth and development, such as the flavonoid biosynthesis pathway. However, CYPs have not yet been systematically studied in the cultivated peanuts (Arachis hypogaea L.), a globally significant cash crop. This study addresses this knowledge deficit through a comprehensive genome-wide analysis, leading to the identification of 589 AhCYP genes in peanuts. Through phylogenetic analysis, all AhCYPs were systematically classified into 9 clans, 43 gene families. The variability in the number of gene family members suggests specialization in biological functions. Intriguingly, both tandem duplication and fragment duplication events have emerged as pivotal drivers in the evolutionary expansion of the AhCYP superfamily. Ka/Ks analysis underscored the substantial influence of strong purifying selection on the evolution of AhCYPs. Furthermore, we selected 21 genes encoding 8 enzymes associated with the flavonoid pathway. The results of quantitative real-time PCR (qRT-PCR) experiments unveiled stage-specific expression patterns during the development of peanut testa, with discernible variations between pink and red testa. Importantly, we identified a direct correlation between gene expression levels and the accumulation of metabolites. These findings offer valuable insights into elucidating the comprehensive functions of AhCYPs and the underlying mechanisms governing the divergent accumulation of flavonoids in testa of different colors.

BSR and Full-Length Transcriptome Approaches Identified Candidate Genes for High Seed Ratio in Camellia vietnamensis

(1) Background: C. vietnamensis is very suitable for growth in the low hilly areas of southern subtropical regions. Under appropriate conditions, the oil yield of C. vietnamensis can reach 1125 kg/ha (the existing varieties can reach 750 kg/ha). Moreover, the fruit of C. vietnamensis is large and the pericarp is thick (>5 cm). Therefore, a high seed ratio has become the main target economic trait for the breeding of C. vietnamensis. (2) Methods: A half-sibling population of C. vietnamensis plants with a combination of high and low seed ratios was constructed by crossing a C. vietnamensis female parent. Bulked segregant RNA analysis and full-length transcriptome sequencing were performed to determine the molecular mechanisms underlying a high seed ratio. (3) Results: Seed ratio is a complex quantitative trait with a normal distribution, which is significantly associated with four other traits of fruit (seed weight, seed number, fruit diameter, and pericarp thickness). Two candidate regions related to high seed ratio (HSR) were predicted. One spanned 140.8−148.4 Mb of chromosome 2 and was associated with 97 seed-yield-related candidate genes ranging in length from 278 to 16,628 bp. The other spanned 35.3−37.3 Mb on chromosome 15 and was associated with 38 genes ranging in length from 221 to 16,928 bp. Using the full-length transcript as a template, a total of 115 candidate transcripts were obtained, and 78 transcripts were predicted to be functionally annotated. The DEGs from two set pairs of cDNA sequencing bulks were enriched to cytochrome p450 CYP76F14 (KOG0156; GO:0055114, HSR4, HSR7), the gibberellin phytohormone pathway (GO:0016787, HSR5), the calcium signaling pathway (GO:0005509, HSR6), the polyubiquitin-PPAR signaling pathway (GO:0005515, HSR2, HSR3), and several main transcription factors (bZIP transcription factor, HSR1) in C. vietnamensis.

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.

Efficient with low-cost removal and adsorption mechanisms of norfloxacin, ciprofloxacin and ofloxacin on modified thermal kaolin: experimental and theoretical studies

Quinolone antibiotics (QNs) can be potential hazard to environment and human. Combination of experimental and theoretical studies was used to analyze the adsorption properties of norfloxacin, ciprofloxacin and ofloxacin on modified thermally activated kaolin (KL). Main factors (calcination temperature, dose, pH, cations and regeneration) affecting the adsorption were discussed. Adsorption processes fit the pseudo-second order kinetic and Langmuir model well. The adsorption removal of norfloxacin, ciprofloxacin and ofloxacin can reach 88.53%, 89.43% and 91.46%, respectively. Cations inhibited adsorption, and AlS-KLB can maintain 80% efficiency in five cycles under optimal conditions. Simulations showed that the materials had good adsorption capacity for QNs, and the "①" of KL had the best capacity. Simulations explain the adsorption mechanism: F, H, O atoms of QNs are covalently bonded to O atoms from KL, Al₂O₃ and Al (OH)₃, C atoms from amorphous carbon and H atoms from C-H and Al (OH)₃. The Al atoms of Al₂O₃ and Al, Si atoms of KL are ionically bonded to F, H, O atoms of QNs. This study shed new light on the removal of QNs by providing low-cost and efficient modified KL and elucidating the adsorption mechanism in conjunction with DFT simulations.

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.

Temporal expression of defence and susceptibility genes and tospovirus accumulation in capsicum chlorosis virus-infected capsicum

Yolo Wonder (YW) and Warlock (W), two capsicum cultivars that are susceptible to capsicum chlorosis virus (CaCV), were compared in terms of symptom development, tospovirus accumulation, and host gene expression during the first 12 days post infection (dpi). Temporal expression of selected early CaCV-response genes was used to gain insights into plant-virus interactions and to identify potential targets for CaCV control. Symptoms developed faster in YW during the first seven days of infection, while systemic symptoms were similar in both cultivars at 10 and 12 dpi. CaCV accumulation was higher in YW at 7 dpi despite a lower titre at 3 dpi. At 12 dpi, virus accumulation was similar for both cultivars. Symptom development appears to be correlated to virus accumulation over time for both cultivars. Chalcone synthase (CHS), cytochrome P450 (CYP), and tetraspanin 8-like (TSP8) genes followed a similar expression pattern over time in both cultivars. The thionin gene showed increased expression in CaCV-infected plants at 12 dpi. The WRKY40 gene showed significant differential expression at all time points in YW, but only at 12 dpi in W. The strongest correlation of temporal gene expression and virus titre was seen for CYP, TSP8, thionin, and WRKY40. CHS and CYP may be involved in symptom development, and TSP8 may be involved in virus movement. CHS, CYP, and TSP8 may be good targets for future overexpression or silencing studies to clarify their functions during virus infection and, potentially, for control of CaCV in capsicum.

Hydrogen/Deuterium Transfer from Anisole to Methoxy Radicals: A Theoretical Study of a Deuterium-Labeled Drug Model

Recently, deuterium-labeled drugs, such as deutetrabenazine, have attracted considerable attention. Consequently, understanding the reaction mechanisms of deuterium-labeled drugs is crucial, both fundamentally and for real applications. To understand the mechanisms of H- and D-transfer reactions, in this study, we used deuterated anisole as a deutetrabenazine model and computationally considered the nuclear quantum effects of protons, deuterons, and electrons. We demonstrated that geometrical differences exist in the partially and fully deuterated methoxy groups and hydrogen-bonded structures of intermediates and transition states due to the H/D isotope effect. The observed geometrical features and electronic structures are ascribable to the different nuclear quantum effects of protons and deuterons. Primary and secondary kinetic isotope effects (KIEs) were calculated for H- and D-transfer reactions from deuterated and undeuterated anisole, with the calculated primary KIEs in good agreement with the corresponding experimental data. These results reveal that the nuclear quantum effects of protons and deuterons need to be considered when analyzing the reaction mechanisms of H- and D-transfer reactions and that a theoretical approach that directly includes nuclear quantum effects is a powerful tool for the analysis of H/D isotope effects in H- and D-transfer reactions.

De novo Transcriptome Assembly of Senna occidentalis Sheds Light on the Anthraquinone Biosynthesis Pathway

Senna occidentalis is an annual leguminous herb that is rich in anthraquinones, which have various pharmacological activities. However, little is known about the genetics of S. occidentalis, particularly its anthraquinone biosynthesis pathway. To broaden our understanding of the key genes and regulatory mechanisms involved in the anthraquinone biosynthesis pathway, we used short RNA sequencing (RNA-Seq) and long-read isoform sequencing (Iso-Seq) to perform a spatial and temporal transcriptomic analysis of S. occidentalis. This generated 121,592 RNA-Seq unigenes and 38,440 Iso-Seq unigenes. Comprehensive functional annotation and classification of these datasets using public databases identified unigene sequences related to major secondary metabolite biosynthesis pathways and critical transcription factor families (bHLH, WRKY, MYB, and bZIP). A tissue-specific differential expression analysis of S. occidentalis and measurement of the amount of anthraquinones revealed that anthraquinone accumulation was related to the gene expression levels in the different tissues. In addition, the amounts and types of anthraquinones produced differ between S. occidentalis and S. tora. In conclusion, these results provide a broader understanding of the anthraquinone metabolic pathway in S. occidentalis.

Water biocatalytic effect attenuates cytochrome P450-mediated carcinogenicity of diethylnitrosamine: A computational insight

The mechanism-based mutagenicity and carcinogenicity of diethylnitrosamine (DEN) are believed to act through interactions with cytochrome P450 (P450) enzymes. DFT calculations to explore the conceivable mechanisms underlying the reaction of P450 with DEN with and without water as a biocatalyst were performed. The results shed light on the biocatalytic role of water in lowering the H-abstraction energy barriers because of the electrostatic effect driven by hydrogen bonding. Our DFT analysis revealed how metabolites are formed in the dealkylation (toxification) and denitrosation (detoxification) pathways. Also, our findings uncovered the active position of DEN vulnerable to P450 interactions. Two factors control the toxification and detoxification rates: the stability of denitrosation products and the HS rebound barrier of the α-pathway. Thus, water biocatalytic attenuation of DEN carcinogenicity was attained by stabilizing denitrosation products and slowing the α-HS rebound process. Docking and MD simulations were performed to assess the binding modes of DEN to P450's active site and to inspect the denitrosation and dealkylation processes, respectively.

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 H₂O₂ 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 OleT_JE, 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 OleT_JE 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 OleT_JE guides a reaction to a chemoselectivity.

Interactive Regulation between Aliphatic Hydroxylation and Aromatic Hydroxylation of Thaxtomin D in TxtC: A Theoretical Investigation

TxtC is an unusual bifunctional cytochrome P450 that is able to perform sequential aliphatic and aromatic hydroxylation of the diketopiperazine substrate thaxtomin D in two distinct sites to produce thaxtomin A. Though the X-ray structure of TxtC complexed with thaxtomin D revealed a binding mode for its aromatic hydroxylation, the preferential hydroxylation site is aliphatic C₁₄. It is thus intriguing to unravel how TxtC accomplishes such two-step catalytic hydroxylation on distinct aliphatic and aromatic carbons and why the aliphatic site is preferred in the hydroxylation step. In this work, by employing molecular docking and molecular dynamics (MD) simulation, we revealed that thaxtomin D could adopt two different conformations in the TxtC active site, which were equal in energy with either the aromatic C₂₀-H or aliphatic C₁₄-H pointing toward the active Cpd I oxyferryl moiety. Further ONIOM calculations indicated that the energy barrier for the rate-limiting hydroxylation step on the aliphatic C₁₄ site was 9.6 kcal/mol more favorable than that on the aromatic C₂₀ site. The hydroxyl group on the monohydroxylated intermediate thaxtomin B C₁₄ site formed hydrogen bonds with Ser280 and Thr385, which induced the l-Phe moiety to rotate around the C_β-C_γ bond of the 4-nitrotryptophan moiety. Thus, it adopted an energetically favorable conformation with aromatic C₂₀ adjacent to the oxyferryl moiety. In addition, the hydroxyl group induced solvent water molecules to enter the active site, which propelled thaxtomin B toward the heme plane and resulted in heme distortion. Based on this geometrical layout, the rate-limiting aromatic hydroxylation energy barrier decreased to 15.4 kcal/mol, which was comparable to that of the thaxtomin D aliphatic hydroxylation process. Our calculations indicated that heme distortion lowered the energy level of the lowest Cpd I α-vacant orbital, which promoted electron transfer in the rate-limiting thaxtomin B aromatic hydroxylation step in TxtC.

Mechanism-Based Inactivation of Cytochrome P450 Enzymes: Computational Insights

Mechanism-based inactivation (MBI) refers to the metabolic bioactivation of a xenobiotic by cytochrome P450s to a highly reactive intermediate which subsequently binds to the enzyme and leads to the quasi-irreversible or irreversible inhibition. Xenobiotics, mainly drugs with specific functional units, are the major sources of MBI. Two possible consequences of MBI by medicinal compounds are drug-drug interaction and severe toxicity that are observed and highlighted by clinical experiments. Today almost all of these latent functional groups (e.g., thiophene, furan, alkylamines, etc.) are known, and their features and mechanisms of action, owing to the vast experimental and theoretical studies, are determined. In the past decade, molecular modeling techniques, mostly density functional theory, have revealed the most feasible mechanism that a drug undergoes by P450 enzymes to generate a highly reactive intermediate. In this review, we provide a comprehensive and detailed picture of computational advances toward the elucidation of the activation mechanisms of various known groups with MBI activity. To this aim, we briefly describe the computational concepts to carry out and analyze the mechanistic investigations, and then, we summarize the studies on compounds with known inhibition activity including thiophene, furan, alkylamines, terminal acetylene, etc. This study can be reference literature for both theoretical and experimental (bio)chemists in several different fields including rational drug design, the process of toxicity prevention, and the discovery of novel inhibitors and catalysts.

How external perturbations affect the chemoselectivity of substrate activation by cytochrome P450 OleTJE

Cytochrome P450 enzymes are versatile biocatalysts found in most forms of life. Generally, the cytochrome P450s react with dioxygen and hence are haem-based mono-oxygenases; however, in specific isozymes, H2O2 rather than O2 is used and these P450s act as peroxygenases. The P450 OleTJE is a peroxygenase that binds long to medium chain fatty acids and converts them to a range of products originating from Cα-hydroxylation, Cβ-hydroxylation, Cα-Cβ desaturation and decarboxylation of the substrate. There is still controversy regarding the details of the reaction mechanism of P450 OleTJE; how the products are formed and whether the product distributions can be influenced by external perturbations. To gain further insights into the structure and reactivity of P450 OleTJE, we set up a range of large active site model complexes as well as full enzymatic structures and did a combination of density functional theory studies and quantum mechanics/molecular mechanics calculations. In particular, the work focused on the mechanisms leading to these products under various reaction conditions. Thus, for a small cluster model, we find a highly selective Cα-hydroxylation pathway that is preferred over Cβ-H hydrogen atom abstraction by at least 10 kcal mol-1. Introduction of polar residues to the model, such as an active site protonated histidine residue or through external electric field effects, lowers the Cβ-H hydrogen atom abstraction barriers are lowered, while a full QM/MM model brings the Cα-H and Cβ-H hydrogen atom abstraction barriers within 1 kcal mol-1. Our studies; therefore, implicate that environmental effects in the second-coordination sphere can direct and guide selectivities in enzymatic reaction mechanisms.

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.

How Does Replacement of the Axial Histidine Ligand in Cytochrome c Peroxidase by Nδ-Methyl Histidine Affect Its Properties and Functions? A Computational Study

Heme peroxidases have important functions in nature related to the detoxification of H₂O₂. They generally undergo a catalytic cycle where, in the first stage, the iron(III)-heme-H₂O₂ complex is converted into an iron(IV)-oxo-heme cation radical species called Compound I. Cytochrome c peroxidase Compound I has a unique electronic configuration among heme enzymes where a metal-based biradical is coupled to a protein radical on a nearby Trp residue. Recent work using the engineered N_δ-methyl histidine-ligated cytochrome c peroxidase highlighted changes in spectroscopic and catalytic properties upon axial ligand substitution. To understand the axial ligand effect on structure and reactivity of peroxidases and their axially N_δ-methyl histidine engineered forms, we did a computational study. We created active site cluster models of various sizes as mimics of horseradish peroxidase and cytochrome c peroxidase Compound I. Subsequently, we performed density functional theory studies on the structure and reactivity of these complexes with a model substrate (styrene). Thus, the work shows that the N_δ-methyl histidine group has little effect on the electronic configuration and structure of Compound I and little changes in bond lengths and the same orbital occupation is obtained. However, the N_δ-methyl histidine modification impacts electron transfer processes due to a change in the reduction potential and thereby influences reactivity patterns for oxygen atom transfer. As such, the substitution of the axial histidine by N_δ-methyl histidine in peroxidases slows down oxygen atom transfer to substrates and makes Compound I a weaker oxidant. These studies are in line with experimental work on N_δ-methyl histidine-ligated cytochrome c peroxidases and highlight how the hydrogen bonding network in the second coordination sphere has a major impact on the function and properties of the enzyme.

Computational Approaches in Preclinical Studies on Drug Discovery and Development

Because undesirable pharmacokinetics and toxicity are significant reasons for the failure of drug development in the costly late stage, it has been widely recognized that drug ADMET properties should be considered as early as possible to reduce failure rates in the clinical phase of drug discovery. Concurrently, drug recalls have become increasingly common in recent years, prompting pharmaceutical companies to increase attention toward the safety evaluation of preclinical drugs. In vitro and in vivo drug evaluation techniques are currently more mature in preclinical applications, but these technologies are costly. In recent years, with the rapid development of computer science, in silico technology has been widely used to evaluate the relevant properties of drugs in the preclinical stage and has produced many software programs and in silico models, further promoting the study of ADMET in vitro. In this review, we first introduce the two ADMET prediction categories (molecular modeling and data modeling). Then, we perform a systematic classification and description of the databases and software commonly used for ADMET prediction. We focus on some widely studied ADMT properties as well as PBPK simulation, and we list some applications that are related to the prediction categories and web tools. Finally, we discuss challenges and limitations in the preclinical area and propose some suggestions and prospects for the future.

Cross-linking of aromatic phenolate groups by cytochrome P450 enzymes: a model for the biosynthesis of vancomycin by OxyB

Aromatic cross-linking by cytochrome P450 enzymes was studied computationally. P450 Compound I rapidly abstracts two weak phenolic H-atoms that link up via a rate-determining C–O bond formation.

Reactivity patterns of vanadium(iv/v)-oxo complexes with olefins in the presence of peroxides: a computational study

Vanadium porphyrin complexes are naturally occurring substances found in crude oil and have been shown to have medicinal properties as well. Little is known on their activities with substrates; therefore, we decided to perform a detailed density functional theory study on the properties and reactivities of vanadium(iv)- and vanadium(v)-oxo complexes with a TPPCl8 or 2,3,7,8,12,13,17,18-octachloro-meso-tetraphenylporphyrinato ligand system. In particular, we investigated the reactivity of [VV(O)(TPPCl8)]+ and [VIV(O)(TPPCl8)] with cyclohexene in the presence of H2O2 or HCO4-. The work shows that vanadium(iv)-oxo and vanadium(v)-oxo are sluggish oxidants by themselves and react with olefins slowly. However, in the presence of hydrogen peroxide, these metal-oxo species can be transformed into a side-on vanadium-peroxo complex, which reacts with substrates more efficiently. Particularly with anionic axial ligands, the side-on vanadium-peroxo and vanadium-oxo complexes produced epoxides from cyclohexene via small barrier heights. In addition to olefin epoxidation, we investigated aliphatic hydroxylation mechanisms by the same oxidants and some oxidants show efficient and viable cyclohexene hydroxylation mechanisms. The work implies that vanadium-oxo and vanadium-peroxo complexes can react with double bonds through epoxidation, and under certain conditions also undergo hydroxylation, but the overall reactivity is highly dependent on the equatorial ligand, the local environment and the presence or absence of anionic axial ligands.

Production of metabolites of the anti-cancer drug noscapine using a P450BM3 mutant library

•

Mutants of P450_BM3 can metabolise noscapine.

•

Noscapine is N-demethylated with high selectivity.

•

The metabolites produced are of interest for drug development.

•

The profile of metabolites generated resembles that of mammalian CYP3A4.

Cytochrome P450 enzymes are a promising tool for the late-stage diversification of lead drug candidates and can provide an alternative route to structural modifications that are difficult to achieve with synthetic chemistry. In this study, a library of P450_BM3 mutants was produced using site-directed mutagenesis and the enzymes screened for metabolism of the opium poppy alkaloid noscapine, a drug with anticancer activity. Of the 18 enzyme mutants screened, 12 showed an ability to metabolise noscapine that was not present in the wild-type enzyme. Five noscapine metabolites were detected by LC-MS/MS, with the major metabolite for all mutants being N-demethylated noscapine. The highest observed regioselectivity for N-demethylation was 88%. Two hydroxylated metabolites, a catechol and two C-C cleavage products were also detected. P450-mediated production of hydroxylated and N-demethylated noscapine structures may be useful for the development of noscapine analogues with improved biological activity. The variation in substrate turnover, coupling efficiency and product distribution between the active mutants was considered alongside in silico docking experiments to gain insight into structural and functional effects of the introduced mutations. Selected mutants were identified as targets for further mutagenesis to improve activity and when coupled with an optimised process may provide a route for the preparative-scale production of noscapine metabolites.

Advances in Sustainable Catalysis: A Computational Perspective

The enormous challenge of moving our societies to a more sustainable future offers several exciting opportunities for computational chemists. The first principles approach to "catalysis by design" will enable new and much greener chemical routes to produce vital fuels and fine chemicals. This prospective outlines a wide variety of case studies to underscore how the use of theoretical techniques, from QM/MM to unrestricted DFT and periodic boundary conditions, can be applied to biocatalysis and to both homogeneous and heterogenous catalysts of all sizes and morphologies to provide invaluable insights into the reaction mechanisms they catalyze.

Interplay Between Steric and Electronic Effects: A Joint Spectroscopy and Computational Study of Nonheme Iron(IV)-Oxo Complexes

Iron is an essential element in nonheme enzymes that plays a crucial role in many vital oxidative transformations and metabolic reactions in the human body. Many of those reactions are regio- and stereospecific and it is believed that the selectivity is guided by second-coordination sphere effects in the protein. Here, results are shown of a few engineered biomimetic ligand frameworks based on the N4Py (N,N-bis(2-pyridylmethyl)-N-bis(2-pyridyl)methylamine) scaffold and the second-coordination sphere effects are studied. For the first time, selective substitutions in the ligand framework have been shown to tune the catalytic properties of the iron(IV)-oxo complexes by regulating the steric and electronic factors. In particular, a better positioning of the oxidant and substrate in the rate-determining transition state lowers the reaction barriers. Therefore, an optimum balance between steric and electronic factors mediates the ideal positioning of oxidant and substrate in the rate-determining transition state that affects the reactivity of high-valent reaction intermediates.

How Does the Oxidation State of Palladium Surfaces Affect the Reactivity and Selectivity of Direct Synthesis of Hydrogen Peroxide from Hydrogen and Oxygen Gases? A Density Functional Study

Direct synthesis of H₂O₂ from H₂ and O₂ is an environmentally benign and atom economic process and as such is the ideal pathway in catalysis. However, currently no low-cost pathway of this kind of catalysis exists, although it would be an attractive alternative strategy to the common industrial anthraquinone method for H₂O₂ production. Metal-based catalysts are widely employed in such a direct synthesis process but often need to be oxidized, alloyed, or supplied with additives to make them selective. To understand the metal-oxidation state in heterogeneous catalysis, we studied the selective oxidation of hydrogen by molecular oxygen on Pd(111) and PdO(101) surfaces, leading to either H₂O₂ or H₂O products. Our results demonstrate, for the first time, that the oxidized PdO(101) surface clearly shows better performance and selectivity, as compared to the reduced Pd(111) one. The activation barrier on the oxidized Pd surface is ca. 0.2 eV lower than the one on the reduced Pd surface. On the oxidized surface, the H₂O₂ synthesis route is preferred, while, on the reduced surface, the H₂O route is predominant. The decomposition of H₂O₂ is also greatly inhibited on the oxidized surface. We analyzed the different pathways in detail through thermochemical cycles, which establishes that the oxidized surface shows weaker adsorption ability toward the reagents O₂ and H₂, the key intermediate OOH, and also the product H₂O₂ in comparison with the Pd(111) surface, which we believe affect the selectivity. The work presented here clearly shows that the oxidation state of metal surfaces is one of the most important factors that tunes the catalysis of a chemical reaction and can affect the selectivity and reaction patterns dramatically.

Mechanistic Studies of Fatty Acid Activation by CYP152 Peroxygenases Reveal Unexpected Desaturase Activity

The majority of cytochrome P450 enzymes (CYPs) predominantly operate as monooxygenases, but recently a class of P450 enzymes was discovered, that can act as peroxygenases (CYP152). These enzymes convert fatty acids through oxidative decarboxylation, yielding terminal alkenes, and through α- and β-hydroxylation to yield hydroxy-fatty acids. Bioderived olefins may serve as biofuels, and hence understanding the mechanism and substrate scope of this class of enzymes is important. In this work, we report on the substrate scope and catalytic promiscuity of CYP OleT_JE and two of its orthologues from the CYP152 family, utilizing α-monosubstituted branched carboxylic acids. We identify α,β-desaturation as an unexpected dominant pathway for CYP OleT_JE with 2-methylbutyric acid. To rationalize product distributions arising from α/β-hydroxylation, oxidative decarboxylation, and desaturation depending on the substrate’s structure and binding pattern, a computational study was performed based on an active site complex of CYP OleT_JE containing the heme cofactor in the substrate binding pocket and 2-methylbutyric acid as substrate. It is shown that substrate positioning determines the accessibility of the oxidizing species (Compound I) to the substrate and hence the regio- and chemoselectivity of the reaction. Furthermore, the results show that, for 2-methylbutyric acid, α,β-desaturation is favorable because of a rate-determining α-hydrogen atom abstraction, which cannot proceed to decarboxylation. Moreover, substrate hydroxylation is energetically impeded due to the tight shape and size of the substrate binding pocket.

Predicting the Initial Thermal Decomposition Path of Nitrobenzene Caused by Mode Vibration at Moderate-Low Temperatures: Temperature-Dependent Anti-Stokes Raman Spectra Experiments and First-Principals Calculations

The lack of understanding of the initial decomposition micromechanism of energetic materials subjected to external stimulation has hindered its safe storage, usage, and development. The initial thermal decomposition path of nitrobenzene triggered by molecular thermal motion is investigated using temperature-dependent anti-Stokes Raman spectra experiments and first-principles calculations to clarify the initial thermal decomposition micromechanism. The experiment shows that the symmetric nitro stretching, antisymmetric nitro stretching, and phenyl ring stretching vibration modes are active as increasing temperature below 500 K. The DFT method is used to examine the effects of the three mode vibrations on the initial decomposition of nitrobenzene by relaxed scan for each relevant change in bond lengths and bond angles to obtain the optimal reaction channel leading to initial thermal decomposition of nitrobenzene. The results demonstrate that the initial thermal decomposition is the isomerization of nitrobenzene to phenyl nitrite. The optimal reaction channel leading to the initial isomerization is the increase or decrease of angle O-N-C from the antisymmetric nitro stretching vibration, which causes the torsion of nitro group and the subsequent oxygen atom attacking carbon atom. The scanning energy barrier related to angle O-N-C is about 62.1 kcal/mol, which is very consistent with the calculated activation barrier of isomerization of nitrobenzene. This proves the reliability of our conclusions.

Combined Multistate and Kohn-Sham Density Functional Theory Studies of the Elusive Mechanism of N-Dealkylation of N,N-Dimethylanilines Mediated by the Biomimetic Nonheme Oxidant FeIV(O)(N4Py)(ClO4)2

The oxidative C-H bond activation mediated by heme and nonheme enzymes and related biomimetics is one of the most interesting processes in bioinorganic and oxidative chemistry. However, the mechanisms of these reactions are still elusive and controversy due to the involvement of highly reactive metal-oxo intermediates with multiple spin states, despite extensive experimental efforts, especially for the N-dealkylation of N,N-dialkyalinines. In this work, we employed multistate density functional theory (MSDFT) and the Kohn-Sham DFT to investigate the mechanism of N-demethylation of N,N-dimethyalinines oxidized by the reaction intermediate Fe^IV(O)(N4Py)(ClO₄)₂. The Kohn-Sham DFT study demonstrated that the reaction proceeds via a rate-limiting hydrogen atom transfer (HAT) step and a subsequent barrier-free oxygen rebound step to form the carbinol product. The MSDFT investigation on the first C-H activation further showed that this step is an initial hydrogen atom abstraction that is highly correlated between CEPT and HAT, i.e., both CEPT and HAT processes make significant contributions to the mechanism before reaching the diabatic crossing point, then the valence bond character of the adiabatic ground state is switched to the CEPT product configuration. The findings from this work may be applicable to other hydrogen abstraction process.

Quantum Mechanics/Molecular Mechanics Studies on the Relative Reactivities of Compound I and II in Cytochrome P450 Enzymes

The cytochromes P450 are drug metabolizing enzymes in the body that typically react with substrates through a monoxygenation reaction. During the catalytic cycle two reduction and protonation steps generate a high-valent iron (IV)-oxo heme cation radical species called Compound I. However, with sufficient reduction equivalents present, the catalytic cycle should be able to continue to the reduced species of Compound I, called Compound II, rather than a reaction of Compound I with substrate. In particular, since electron transfer is usually on faster timescales than atom transfer, we considered this process feasible and decided to investigate the reaction computationally. In this work we present a computational study using density functional theory methods on active site model complexes alongside quantum mechanics/molecular mechanics calculations on full enzyme structures of cytochrome P450 enzymes. Specifically, we focus on the relative reactivity of Compound I and II with a model substrate for O⁻H bond activation. We show that generally the barrier heights for hydrogen atom abstraction are higher in energy for Compound II than Compound I for O⁻H bond activation. Nevertheless, for the activation of such bonds, Compound II should still be an active oxidant under enzymatic conditions. As such, our computational modelling predicts that under high-reduction environments the cytochromes P450 can react with substrates via Compound II but the rates will be much slower.

The role of spin states in the catalytic mechanism of the intra- and extradiol cleavage of catechols by O2

Iron-dependent enzymes and biomimetic iron complexes can catalyze the ring cleavage of very inert, aromatic compounds. The mechanisms of these transformations and the factors that lead either to extradiol cleavage or intradiol cleavage have not been fully understood. By using density functional theory we have elucidated the mechanism of the catalytic cycle for two biomimetic complexes, and explained the difference in the experimentally obtained products.

Insights into the NHC-catalyzed cascade Michael/aldol/lactamization reaction: mechanism and origin of stereoselectivity

A DFT study of the NHC-catalyzed cascade Michael/aldol/lactamization reaction has been performed for the first time.