Benzyne Formation in the Mechanism-Based Inactivation of Cytochrome P450 by 1-Aminobenzotriazole and N-Benzyl-1-Aminobenzotriazole: Computational Insights

Mechanistic insights into the metabolic pathways of vanillin: unraveling cytochrome P450 interaction mechanisms and implications for food safety

Vanillin, a key flavor compound found in vanilla beans, is widely used in the food and pharmaceutical industries for its aromatic properties and potential therapeutic benefits. This study presents a comprehensive quantum chemical analysis to elucidate the interaction mechanisms of vanillin with CYP450 enzymes, with a focus on mechanism-based inactivation. Three potential inactivation pathways were evaluated: aldehyde deformylation, methoxy dealkylation, and acetal formation. Aldehyde deformylation was identified as the most energy-efficient, involving the removal of the aldehyde group from vanillin and leading to the formation of benzyne intermediates that could react with the iron porphyrin moiety of CYP450, potentially resulting in enzyme inactivation. Further investigation into the interactions of vanillin with CYP2E1 and CYP1A2 was conducted using molecular docking and molecular dynamics (MD) simulation. The docking analyses supported the findings from DFT studies, wherein vanillin revealed high binding affinities with the studied isozymes. Moreover, vanillin occupied the main binding site in both isozymes, as evidenced by the inclusion of the heme moiety in their binding mechanisms. Employing a 100 ns molecular dynamics simulation, we scrutinized the interaction dynamics between vanillin and the two isozymes of CYP450. The assessment of various MD parameters along with interaction energies revealed that vanillin exhibited stable trajectories and substantial energy stabilization during its interaction with both CYP450 isozymes. These insights can guide future research and ensure the safe application of vanillin, especially in scenarios where it may interact with CYP450 enzymes.

Substrate Conformational Switch Enables the Stereoselective Dimerization in P450 NascB: Insights from Molecular Dynamics Simulations and Quantum Mechanical/Molecular Mechanical Calculations

P450 NascB catalyzes the coupling of cyclo-(l-tryptophan-l-proline) (1) to generate (-)-naseseazine C (2) through intramolecular C-N bond formation and intermolecular C-C coupling. A thorough understanding of its catalytic mechanism is crucial for the engineering or design of P450-catalyzed C-N dimerization reactions. By employing MD simulations, QM/MM calculations, and enhanced sampling, we assessed various mechanisms from recent works. Our study demonstrates that the most favorable pathway entails the transfer of a hydrogen atom from N7-H to Cpd I. Subsequently, there is a conformational change in the substrate radical, shifting it from the Re-face to the Si-face of N7 in Substrate 1. The Si-face conformation of Substrate 1 is stabilized by the protein environment and the π-π stacking interaction between the indole ring and heme porphyrin. The subsequent intermolecular C3-C6' bond formation between Substrate 1 radical and Substrate 2 occurs via a radical attack mechanism. The conformational switch of the Substrate 1 radical not only lowers the barrier of the intermolecular C3-C6' bond formation but also yields the correct stereoselectivity observed in experiments. In addition, we evaluated the reactivity of the ferric-superoxide species, showing it is not reactive enough to initiate the hydrogen atom abstraction from the indole NH group of the substrate. Our simulation provides a comprehensive mechanistic insight into how the P450 enzyme precisely controls both the intramolecular C-N cyclization and intermolecular C-C coupling. The current findings align with the available experimental data, emphasizing the pivotal role of substrate dynamics in governing P450 catalysis.

Deciphering the Molecular Mechanisms of Reactive Metabolite Formation in the Mechanism-Based Inactivation of Cytochrome p450 1B1 by 8-Methoxypsoralen and Assessing the Driving Effect of phe268

This study provides a comprehensive computational exploration of the inhibitory activity and metabolic pathways of 8-methoxypsoralen (8-MP), a furocoumarin derivative used for treating various skin disorders, on cytochrome P450 (P450). Employing quantum chemical DFT calculations, molecular docking, and molecular dynamics (MD) simulations analyses, the biotransformation mechanisms and the active site binding profile of 8-MP in CYP1B1 were investigated. Three plausible inactivation mechanisms were minutely scrutinized. Further analysis explored the formation of reactive metabolites in subsequent P450 metabolic processes, including covalent adduct formation through nucleophilic addition to the epoxide, 8-MP epoxide hydrolysis, and non-CYP-catalyzed epoxide ring opening. Special attention was paid to the catalytic effect of residue Phe268 on the mechanism-based inactivation (MBI) of P450 by 8-MP. Energetic profiles and facilitating conditions revealed a slight preference for the C4'=C5' epoxidation pathway, while recognizing a potential kinetic competition with the 8-OMe demethylation pathway due to comparable energy demands. The formation of covalent adducts via nucleophilic addition, particularly by phenylalanine, and the generation of potentially harmful reactive metabolites through autocatalyzed ring cleavage are likely to contribute significantly to P450 metabolism of 8-MP. Our findings highlight the key role of Phe268 in retaining 8-MP within the active site of CYP1B1, thereby facilitating initial oxygen addition transition states. This research offers crucial molecular-level insights that may guide the early stages of drug discovery and risk assessment related to the use of 8-MP.

Improved methods for the detection of heme and protoporphyrin IX adducts and quantification of heme B from cytochrome P450 containing systems

Heme B is a critical prosthetic group for the function of numerous proteins including the cytochrome P450 (CYP) family of enzymes. CYP enzymes are involved in the metabolism of endogenous and xenobiotic molecules that are of central interest in drug development. Formation of reactive metabolites by CYPs can lead to heme modification and destruction of the enzyme. The structure of the adducted heme can provide key information on the mechanism of inactivation, which is of great interest during preclinical drug discovery. Historically, techniques to extract the modified heme or protoporphyrin IX species involved harsh extraction conditions and esterification of propionate groups to aid chromatography. We have developed a simplified extraction method and LC/MS chromatography system that does not require derivatization to quantify heme B and identify modified heme B species from multiple CYP-containing matrices. The method uses mass defect filter triggered data dependent MS2 scans to rapidly identify heme and protoporphyrin IX adducts. These methods may also be useful for the analysis of other heme variants and hemoproteins.

Mechanistic aspects of reactive metabolite formation in clomethiazole catalyzed biotransformation by cytochrome P450 enzymes

Clomethiazole (CLM), a sedative and anticonvulsant drug, is commonly employed for the treatment of alcohol withdrawal syndrome because it suppresses cytochrome P450 (P450) activity associated with the generation of free radicals and liver damage. The catalyzed biotransformation of thiazole-containing drugs by P450 is known to afford reactive metabolites. These metabolites can alter the biological functions of macromolecules and result in toxicity and adverse drug interactions. Multitargeted molecular modeling and quantum chemical DFT calculations were performed to explore the binding modes and molecular mechanisms underlying the mechanism-based inactivation (MBI) of P450 by CLM. The mechanistic details associated with reactive metabolite formation from further metabolic processes were extensively assessed. Seven possible routes were proposed for CLM-P450 biotransformation including CLM hydroxylation, sulfoxidation, N-oxidation, CN epoxidation (oxaziridine formation), and CC epoxidation. The results revealed a degree of preference for the C-N epoxidation pathway because of the low energy requirements of its rate-determining step (8.74 and 10.07 kcal mol-1 for LS and HS states, respectively). A kinetic competition for the CLM-methyl hydroxylation pathway was detected because the H-abstraction energy barrier was relatively comparable to the thermodynamically prevailing oxaziridine formation rate-determining step (12.58 and 14.52 kcal mol-1 for quartet and doublet states, respectively). Our studies assessed the mechanisms of covalent nucleophilic epoxide adduct formation through nucleophilic addition, hydrolysis of epoxidation products, and nonenzymatic degradation. CLM was shown to display P450-inhibitory activity by forming covalent adducts rather than further metabolization to reactive metabolites. The outcomes of molecular docking allowed assessing the binding profile of CLM with three human P450 isozymes, namely, CYP2E1, CYP3A4, and CYP2D6.

Deciphering the Selectivity of the Electrochemical CO2 Reduction to CO by a Cobalt Porphyrin Catalyst in Neutral Aqueous Solution: Insights from DFT Calculations

Density functional theory (DFT) calculations were conducted to investigate the cobalt porphyrin-catalyzed electro-reduction of CO2 to CO in an aqueous solution. The results suggest that CoII -porphyrin (CoII -L) undertakes a ligand-based reduction to generate the active species CoII -L⋅- , where the CoII center antiferromagnetically interacts with the ligand radical anion. CoII -L⋅- then performs a nucleophilic attack on CO2 , followed by protonation and a reduction to give CoII -L-COOH. An intermolecular proton transfer leads to the heterolytic cleavage of the C-O bond, producing intermediate CoII -L-CO. Subsequently, CO is released from CoII -L-CO, and CoII -L is regenerated to catalyze the next cycle. The rate-determining step of this CO2 RR is the nucleophilic attack on CO2 by CoII -L⋅- , with a total barrier of 20.7 kcal mol-1 . The competing hydrogen evolution reaction is associated with a higher total barrier. A computational investigation regarding the substituent effects of the catalyst indicates that the CoPor-R3 complex is likely to display the highest activity and selectivity as a molecular catalyst.

Energy Decomposition Analysis of the Nature of Coordination Bonding at the Heme Iron Center in Cytochrome P450 Inhibition

Drug compounds or their metabolic intermediates (MIs) sometimes inhibit the function of cytochrome P450 enzymes (P450s) by forming a coordination bond with the Fe(III) heme or Fe(II) heme of P450s. Such inhibition is one of the major causes of drug-drug interactions (DDIs), a subject of longstanding academic and practical interest. However, such coordination bonding is not fully understood at the quantum mechanical level, thus hampering rational improvement of the accuracy of DDI-related predictions. In this work, we employed density functional theory (DFT) and the generalized Kohn-Sham energy decomposition analysis (GKS-EDA) scheme to investigate the nature of the coordination bonding formed in the reversible and quasi-irreversible inhibition of P450s. The GKS-EDA results highlighted a previously unrecognized role of the electron correlation effect in P450 inhibition. The correlation effect tends to be larger in Fe(II) complexes of MI-type inhibitors and is particularly prominent for the nitrosoalkane ligand. An additional natural bond orbital (NBO) analysis provided insight into the relative significance of the σ donation and π backdonation effects in various heme-inhibitor complexes.

The quasi-irreversible inactivation of cytochrome P450 enzymes by paroxetine: a computational approach

The mechanism-based inactivation (MBI) of P450 by paroxetine was investigated by computational analysis. The drug-enzyme interactions were figured out through studying energy profiles of three competing mechanisms. The potency of paroxetine as P450's inhibitor was estimated based on the availability of two active sites for the MBI in the paroxetine structure. The inactivation by the amino site of paroxetine mainly proceeds via the hydrogen atom transfer pathway because of the lower energy demand of its rate determining step. In addition, the low-spin state is the predominant route in the MBI at the methylenedioxo active site as a result of being rebound barrier-free mechanism. Our comparative investigation showed that inactivation at the secondary amine is thermodynamically more favorable because of the lower energy barrier of the dehydration mechanism of the hydroxylated paroxetine complex than its methylenedioxo counterpart. The results of docking analysis coincided with the outputs of DFT calculations since the docking pose with the lowest binding affinity is that for conformation with polar interaction between the amino group of paroxetine and the oxo moiety of P450's active site. Assessment of the molecular dynamics simulations trajectories revealed the favorable interaction of paroxetine with P450.

Reactive Metabolites from Thiazole-Containing Drugs: Quantum Chemical Insights into Biotransformation and Toxicity

Drugs containing thiazole and aminothiazole groups are known to generate reactive metabolites (RMs) catalyzed by cytochrome P450s (CYPs). These RMs can covalently modify essential cellular macromolecules and lead to toxicity and induce idiosyncratic adverse drug reactions. Molecular docking and quantum chemical hybrid DFT study were carried out to explore the molecular mechanisms involved in the biotransformation of thiazole (TZ) and aminothiazole (ATZ) groups leading to RM epoxide, S-oxide, N-oxide, and oxaziridine. The energy barrier required for the epoxidation is 13.63 kcal/mol, that is lower than that of S-oxidation, N-oxidation, and oxaziridine formation (14.56, 17.90, and 20.20, kcal/mol respectively). The presence of the amino group in ATZ further facilitates all the metabolic pathways, for example, the barrier for the epoxidation reaction is reduced by ∼2.5 kcal/mol. Some of the RMs/their isomers are highly electrophilic and tend to form covalent bonds with nucleophilic amino acids, finally leading to the formation of metabolic intermediate complexes (MICs). The energy profiles of these competitive pathways have also been explored.

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.

Computational Insight into the Activation Mechanism of Carcinogenic N'-Nitrosonornicotine (NNN) Catalyzed by Cytochrome P450

Tobacco-specific N'-nitrosonornicotine (NNN), a genotoxic nitrosamine classified as Group 1 carcinogen, is also present in atmospheric particulate matter and has even been detected as a new disinfection byproduct in wastewaters. NNN generally requires metabolic activation by cytochrome P450 enzymes to exert its genotoxicity, but the respective biotransformation pathways have not been described in detail. In this work, we performed density functional theory (DFT) calculations to unravel possible NNN activation pathways including α-hydroxylation, β-hydroxylation, pyridine N-oxidation, and norcotinine formation. The results reveal an initial rate-determining H_α-atom abstraction step for α-hydroxylation, followed by an unexpected kinetic competition between denitrosation and OH rebound, leading to ( iso-)myosmine as a detoxified product and α-hydroxyNNNs as the precursor of carcinogenic diazohydroxides, respectively. Further detoxification routes are given by β-hydroxylation with relative high reaction barrier and N-oxidation with comparable barrier to the toxifying α-hydroxylation. Moreover, we show for the first time how norcotinine can be generated as a minor NNN metabolite that is formed from iso-myosmine through a unique porphyrin-assisted H atom 1,2-transfer mechanism. These results demonstrate that the carcinogenic potential of NNN is subject to a kinetic competition between activating and deactivating metabolic routes, and identify respective biomarkers to inform about the individual risk associated with NNN exposure.

Theoretical Study of the Mechanism of the Nonheme Iron Enzyme EgtB

EgtB is a nonheme iron enzyme catalyzing the C-S bond formation between γ-glutamyl cysteine (γGC) and N-α-trimethyl histidine (TMH) in the ergothioneine biosynthesis. Density functional calculations were performed to elucidate and delineate the reaction mechanism of this enzyme. Two different mechanisms were considered, depending on whether the sulfoxidation or the S-C bond formation takes place first. The calculations suggest that the S-O bond formation occurs first between the thiolate and the ferric superoxide, followed by homolytic O-O bond cleavage, very similar to the case of cysteine dioxygenase. Subsequently, proton transfer from a second-shell residue Tyr377 to the newly generated iron-oxo moiety takes place, which is followed by proton transfer from the TMH imidazole to Tyr377, facilitated by two crystallographically observed water molecules. Next, the S-C bond is formed between γGC and TMH, followed by proton transfer from the imidazole CH moiety to Tyr377, which was calculated to be the rate-limiting step for the whole reaction, with a barrier of 17.9 kcal/mol in the quintet state. The calculated barrier for the rate-limiting step agrees quite well with experimental kinetic data. Finally, this proton is transferred back to the imidazole nitrogen to form the product. The alternative thiyl radical attack mechanism has a very high barrier, being 25.8 kcal/mol, ruling out this possibility.

Mechanistic study of allopurinol oxidation using aldehyde oxidase, xanthine oxidase and cytochrome P450 enzymes

The oxidation reaction of allopurinol to its active metabolite (oxypurinol) is investigated using the AO and P450 enzymes. To the contrary of AO (and XO), the P450 enzyme can metabolize the allopurinol with a not self-inhibitory mechanism.