Fundamental reaction pathways for cytochrome P450-catalyzed 5'-hydroxylation and N-demethylation of nicotine

New insight into biodegradation mechanism of phenylurea herbicides by cytochrome P450 enzymes: Successive N-demethylation mechanism

Phenylurea herbicides (PUHs) present one of the most important herbicides, which have cause serious effects on ecological environment and humans. Nowadays enzyme strategy shows great advantages in degradation of PUHs. Here density functional theory (DFT), quantitative structure - activity relationship (QSAR) and quantum mechanics/molecular mechanics (QM/MM) approaches are used to investigate the degradation mechanism of PUHs catalyzed by P450 enzymes. Two successive N-demethylation pathways are identified and two hydrogen abstraction (H-abstraction) reaction pathways are identified as the rate-determining step through high-throughput DFT calculations. The Boltzmann-weighted average energy barrier of the second H-abstraction pathway (19.95 kcal/mol) is higher than that of the first H-abstraction pathway (16.80 kcal/mol). Two QSAR models are established to predict the energy barriers of the two H-abstraction pathways based on the quantum chemical descriptors and mordred molecular descriptors. The determination coefficient (R2) values of QSAR models are > 0.9, which reveal that the established QSAR models have great predictive capability. QM/MM calculations indicate that human P450 enzymes are more efficient in degradation of PUHs than crop and weed P450 enzymes. Correlations between energy barriers and key structural/charge parameters are revealed and key parameters that have influence on degradation efficiency of PUHs are identified. This study provides lateral insights into the biodegradation strategy and removal method of PUHs and valuable information for designing or engineering of highly efficient degradation enzymes and genetically modified crops.

Biotransformation of Bisphenol by Human Cytochrome P450 2C9 Enzymes: A Density Functional Theory Study

Bisphenol A (BPA, 2,2-bis-(4-hydroxyphenyl)propane) is used as a precursor in the synthesis of polycarbonate and epoxy plastics; however, its availability in the environment is causing toxicity as an endocrine-disrupting chemical. Metabolism of BPA and their analogues (substitutes) is generally performed by liver cytochrome P450 enzymes and often leads to a mixture of products, and some of those are toxic. To understand the product distributions of P450 activation of BPA, we have performed a computational study into the mechanisms and reactivities using large model structures of a human P450 isozyme (P450 2C9) with BPA bound. Density functional theory (DFT) calculations on mechanisms of BPA activation by a P450 compound I model were investigated, leading to a number of possible products. The substrate-binding pocket is tight, and as a consequence, aliphatic hydroxylation is not feasible as the methyl substituents of BPA cannot reach compound I well due to constraints of the substrate-binding pocket. Instead, we find low-energy pathways that are initiated with phenol hydrogen atom abstraction followed by OH rebound to the phenolic ortho- or para-position. The barriers of para-rebound are well lower in energy than those for ortho-rebound, and consequently, our P450 2C9 model predicts dominant hydroxycumyl alcohol products. The reactions proceed through two-state reactivity on competing doublet and quartet spin state surfaces. The calculations show fast and efficient substrate activation on a doublet spin state surface with a rate-determining electrophilic addition step, while the quartet spin state surface has multiple high-energy barriers that can also lead to various side products including C⁴-aromatic hydroxylation. This work shows that product formation is more feasible on the low spin state, while the physicochemical properties of the substrate govern barrier heights of the rate-determining step of the reaction. Finally, the importance of the second-coordination sphere is highlighted that determines the product distributions and guides the bifurcation pathways.

The Antioxidant Properties of Alfalfa (Medicago sativa L.) and Its Biochemical, Antioxidant, Anti-Inflammatory, and Pathological Effects on Nicotine-Induced Oxidative Stress in the Rat Liver

Medicago sativa Linn or alfalfa is a tonic plant rich in proteins, vitamins, and minerals that is used to treat many diseases due to its pharmacological properties such as anti-inflammatory and antioxidant activities. So, the aim of this study was to evaluate the efficacy of alfalfa methanolic extract (AME) on the prevention of liver damage caused by nicotine. The total phenols, flavonoids levels, and the free radical scavenging activity of its extract (IC50) were measured. In this study, 30 Wistar rats were randomly divided into 5 groups as control (untreated), N (nicotine only), T1, T2, and T3 (nicotine + AME 100, 250, and 500 mg/kg/day, respectively). AME (orally) and nicotine (intraperitoneal injection, 0.5 mg/kg/day) were then administered for 21 days. Weight gain, the liver-to-body weight ratio, liver functional enzymes, and the lipid profile were measured. Moreover, we evaluated oxidative stress, proinflammatory parameters, and histopathological changes in the liver. Total phenols, flavonoids, and IC50 were determined as $51.68\pm 0.62$  mg GAE/g, $18.55\pm 1.01$  mg QE/g, and $350.91\pm 16.46$  μg/ml, respectively. Nicotine changed the measured parameters to abnormal. AME increased weight gain, the liver-to-body weight ratio, and enzymatic antioxidant levels and decreased malondialdehyde, liver functional enzymes, and proinflammatory cytokine levels. The lipid profile and histopathological changes have also been improved by AME in a dose-dependent manner. The results showed that AME in a dose-dependent manner by improving the inflammation and oxidative damage could improve the liver damage caused by nicotine.

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.

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.

How do mutations affect the structural characteristics and substrate binding of CYP21A2? An investigation by molecular dynamics simulations

CYP21A2 mutations affect the activity of the protein leading to CAH disease.

Insight into the Interaction Mechanism of Nicotine, NNK, and NNN with Cytochrome P450 2A13 Based on Molecular Dynamics Simulation

Tobacco smoke contains various cancer-causing toxic substances, including nicotine and nitrosamines 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) and N'-nitrosonornicotine (NNN). The cytochrome 2A13 is involved in nicotine metabolism and in the activation of the pro-carcinogenic agents NNK and NNN, by means of α-hydroxylation reactions. Despite the significance of cytochrome 2A13 in the biotransformation of these molecules, its conformational mechanism and the molecular basis involved in the process are not fully understood. In this study, we used molecular dynamics and principal component analysis simulations for an in-depth analysis of the essential protein motions involved in the interaction of cytochrome 2A13 with its substrates. We also evaluated the interaction of these substrates with the amino acid residues in the binding pocket of cytochrome 2A13. Furthermore, we quantified the nature of these chemical interactions from free energy calculations using the Molecular Mechanics/Generalized Born Surface Area method. The ligands remained favorably oriented toward compound I (cytochrome P450 O═Fe^IV state), to undergo α-hydroxylation. The hydrogen bond with asparagine 297 was essential to maintaining the substrates in a favorable catalytic orientation. The plot of first principal motion vs second principal motion revealed that the enzyme's interaction with nicotine and NNK involved different conformational subgroups, whereas the conformational subgroups in the interaction with NNN are more similar. These results provide new mechanistic insights into the mode of interaction of the substrates with the active site of cytochrome 2A13, in the presence of compound I, which is essential for α-hydroxylation.

DFT investigation on the metabolic mechanisms of theophylline by cytochrome P450 monooxygenase

Theophylline, one of the most commonly used bronchodilators and respiratory stimulators for the treatment of acute and chronic asthmatic conditions, can cause permanent neurological damage through chronic or excessive ingestion. In this work, DFT calculation was performed to identify the metabolic mechanisms of theophylline by cytochrome P450 (CYP450) monooxygenase. Two main metabolic pathways were investigated, namely, N₁- (path A) and N₃- (path B) demethylations, which proceeded through N-methyl hydroxylation followed by the decomposition of the generated carbinolamine species. N-methyl hydroxylation involved a hydrogen atom transfer (HAT) mechanism, which can be generalized as the N-demethylation mechanism of xanthine derivatives. The energy gap between the low-spin double state (LS) and the high-spin quartet state (HS) was low (<1 kcal mol^-1), indicating a two-state reactivity (TSR) mechanism. The generated carbinolamine species preferred to decompose through the adjacent heteroatom (O₆ for path A and O₂ for path B) mediated mechanism. Path B was kinetically more feasible than path A attributed to its relatively lower activation energy. 1-Methylxanthine therefore was the energetically favorable metabolite of theophylline. The observations obtained in the work were in agreement with the experimental observation, which can offer important implications for further pharmacological and clinic studies.

Theoretical study on the N-demethylation mechanism of theobromine catalyzed by P450 isoenzyme 1A2

Theobromine, a widely consumed pharmacological active substance, can cause undesirable muscle stiffness, nausea and anorexia in high doses ingestion. The main N-demethylation metabolic mechanism of theobromine catalyzed by P450 isoenzyme 1A2 (CYP1A2) has been explored in this work using the unrestricted hybrid density functional method UB3LYP in conjunction with the LACVP(Fe)/6-31G (H, C, N, O, S, Cl) basis set. Single-point calculations including empirical dispersion corrections were carried out at the higher 6-311++G** basis set. Two N-demethylation pathways were characterized, i.e., 3-N and 7-N demethylations, which involve the initial N-methyl hydroxylation to form carbinolamines and the subsequent carbinolamines decomposition to yield monomethylxanthines and formaldehydes. Our results have shown that the rate-limiting N-methyl hydroxylation occurs via a hydrogen atom transfer (HAT) mechanism, which proceeds in a spin-selective mechanism (SSM) in the gas phase. The carbinolamines generated are prone to decomposition via the contiguous heteroatom-assisted proton-transfer. Strikingly, 3-N demethylation is more favorable than 7-N demethylation due to its lower free energy barrier and 7-methylxanthine therefore is the optimum product reported for the demethylation of theobromine catalyzed by CYP1A2, which are in good agreement with the experimental observation. This work has first revealed the detail N-demethylation mechanisms of theobromine at the theoretical level. It can offer more significant information for the metabolism of purine alkaloid.

Catalytic mechanism of cytochrome P450 for N-methylhydroxylation of nicotine: reaction pathways and regioselectivity of the enzymatic nicotine oxidation

The fundamental reaction mechanism of cytochrome P450 2A6 (CYP2A6)-catalyzed N-methylhydroxylation of (S)-(-)-nicotine and the free energy profile have been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical (QM/MM) reaction-coordinate calculations. In the CYP2A6-(S)-(-)-nicotine binding structures that allow for 5'-hydroxylation, the N-methyl group is also sufficiently close to the oxygen of Cpd I for the N-methylhydroxylation reaction to occur. It has been demonstrated that the CYP2A6-catalyzed N-methylhydroxylation reaction is a concerted process involving a hydrogen-transfer transition state on both the quartet and the doublet states. The N-methylhydroxylation reaction proceeds mainly in the doublet state, since the free energy barriers on the doublet state are lower than the corresponding ones on the quartet state. The calculated free energy barriers indicate that (S)-(-)-nicotine oxidation catalyzed by CYP2A6 proceeds with a high regioselective abstraction of the hydrogen at the 5'-position, rather than the hydrogen at the N-methyl group. The predicted regioselectivity of 93% is in agreement with the most recent experimentally reported regioselectivity of 95%. The binding mode of (S)-(-)-nicotine in the active site of CYP2A6 is an important determinant for the stereoselectivity of nicotine (S)-(-)-oxidation, whereas the regioselectivity of (S)-(-)-nicotine oxidation is determined mainly by the free energy barrier difference between the 5'-hydroxylation and N-methylhydroxylation reactions.

Catalytic mechanism of cytochrome P450 for 5'-hydroxylation of nicotine: fundamental reaction pathways and stereoselectivity

A series of computational methods were used to study how cytochrome P450 2A6 (CYP2A6) interacts with (S)-(-)-nicotine, demonstrating that the dominant molecular species of (S)-(-)-nicotine in CYP2A6 active site exists in the free base state (with two conformations, SR(t) and SR(c)), despite the fact that the protonated state is dominant for the free ligand in solution. The computational results reveal that the dominant pathway of nicotine metabolism in CYP2A6 is through nicotine free base oxidation. Further, first-principles quantum mechanical/molecular mechanical free energy (QM/MM-FE) calculations were carried out to uncover the detailed reaction pathways for the CYP2A6-catalyzed nicotine 5'-hydroxylation reaction. In the determined CYP2A6-(S)-(-)-nicotine binding structures, the oxygen of Compound I (Cpd I) can abstract a hydrogen from either the trans-5'- or the cis-5'-position of (S)-(-)-nicotine. CYP2A6-catalyzed (S)-(-)-nicotine 5'-hydroxylation consists of two reaction steps, that is, the hydrogen transfer from the 5'-position of (S)-(-)-nicotine to the oxygen of Cpd I (the H-transfer step), followed by the recombination of the (S)-(-)-nicotine moiety with the iron-bound hydroxyl group to generate the 5'-hydroxynicotine product (the O-rebound step). The H-transfer step is rate-determining. The 5'-hydroxylation proceeds mainly with the stereoselective loss of the trans-5'-hydrogen, that is, the 5'-hydrogen trans to the pyridine ring. The calculated overall stereoselectivity of ∼97% favoring the trans-5'-hydroxylation is close to the observed stereoselectivity of 89-94%. This is the first time it has been demonstrated that a CYP substrate exists dominantly in one protonation state (cationic species) in solution, but uses its less-favorable protonation state (neutral free base) to perform the enzymatic reaction.