Hydrogen Abstraction of Camphor Catalyzed by Cytochrome P450cam: A QM/MM Study

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.

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.

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.

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

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.

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.

Stereo-preference of camphor for H-bonding with phenol, methanol and chloroform: A combined matrix isolation IR spectroscopic and quantum chemical investigation

Camphor is known to be held in the substrate pocket of cytochrome P450cam enzyme via H-bond with a tyrosine residue of the enzyme in a unique orientation. This structural exclusivity results in regio- and stereo-specific hydroxylation of camphor by the enzyme. We have carried out a combined IR spectroscopic and quantum chemical investigation to shed light on the factors influencing the conformational exclusivity of 1R-(+)-camphor in the substrate pocket of Cytochrome P450cam, and to determine whether the selectivity is an inherent property of the substrate itself, or is imposed by the enzyme. For this purpose, complexes of camphor have been studied with three H-bond donors namely phenol, methanol and chloroform. Each of the three donors was found to form stable complexes with two distinct conformers; the one mimicking the conformation in enzyme substrate pocket was found to be more stable of the two, for all three donors. Experimentally, both conformers of the H-bonded complexes were identified separately for phenol and methanol in an argon matrix at 8 K, but not for chloroform due to very small energy barrier for interconversion of the two conformers. In room temperature solution phase spectra of camphor with all three donors, the differences in spectral attributes between the two isomeric H-bonded complexes were lost due to thermal motions.

Oxidative Rearrangement Mechanism of Pentalenolactone F Catalyzed by Cytochrome P450 CYP161C2 (PntM)

The CYP161C2 (PntM) from Streptomyces arenae is a member of the cytochrome P450 enzymes, which catalyzes the unusual oxidative rearrangement of pentalenolactone F (1) to the sesquiterpenoid antibiotic pentalenolactone (3). On the basis of the crystal structure of PntM bound with substrate, quantum mechanical/molecular mechanics (QM/MM) calculations have been performed to explore the detailed mechanism of PntM-catalyzed oxidative rearrangement. The conversion from pentalenolactone F (1) to pentalenolactone (3) involves the stereospecific removal of the H-1 si from 1, the syn-1,2-migration of the 2 si methyl group, and the antarafacial loss of H-3 re. The abstraction of H-1 si by Cpd I is calculated to be rate limiting with an energy barrier of 20.3 kcal/mol, which basically agrees with the estimated free energy barrier from experiments (18.6 kcal/mol). It is the unfavorable geometry of Fe-OH-C1 that blocks the oxygen rebound reaction, and the subsequent intramolecular syn-1,2-methyl migration is accompanied by an electron transfer from the substrate to the porphyrin ring via an Fe-OH group, generating the carbocation intermediate. Owing to the positive charge, the intermediate can easily lose a proton to form the final products. Our calculation results indicate that both the carboxyl group of porphyrin and Fe-OH can act as bases to accept the proton of the substrate. The target product pentalenolactone and the three isomeric byproducts correspond to four different modes of deprotonation.