Can a single oxidant with two spin states masquerade as two different oxidants? A study of the sulfoxidation mechanism by cytochrome p450

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.

Substrate Conformation Regulates Aromatic C-H Vs C-F Bond Activation in Heme-Dependent Tyrosine Hydroxylase

A recently discovered heme-dependent enzyme tyrosine hydroxylase (TyrH) offers a green approach for functionalizing the high-strength C-H and C-F bonds in aromatic compounds. However, there is ambiguity regarding the nature of the oxidant (compound 0 or compound I) involved in activating these bonds. Herein, using comprehensive molecular dynamics (MD) simulations and hybrid quantum mechanical/molecular mechanical calculations, we reveal that it is compound I (Cpd I) that acts as the primary oxidant involved in the functionalization of both C-F and C-H bonds. The energy barrier for C-H and C-F activation using compound 0 (Cpd 0) as an oxidant was very high, indicating that Cpd 0 cannot be an oxidant. Consistent with the previous experimental finding, our simulation shows two different conformations of the substrate, where one orientation favors the C-H activation, while the other conformation prefers the C-F activation. As such, our mechanistic study shows that nature utilizes just one oxidant, that is, Cpd I, but it is the active site conformation that decides whether it selects C-F or C-H functionalization which may resemble involvement of two different oxidants.

Cytochrome P450-catalyzed oxidation of halogen-containing substrates

Cytochrome P450 (CYP) enzymes are heme-thiolate monooxygenases which catalyze the oxidation of aliphatic and aromatic C-H bonds and other reactions. The oxidation of halogens by cytochrome P450 enzymes has also been reported. Here we use CYP199A4, from the bacterium Rhodopseudomonas palustris strain HaA2, with a range of para-substituted benzoic acid ligands, which contain halogens, to assess if this enzyme can oxidize these species or if the presence of these electronegative atoms can alter the outcome of P450-catalyzed reactions. Despite binding to the enzyme, there was no detectable oxidation of any of the 4-halobenzoic acids. CYP199A4 was, however, able to efficiently catalyze the oxidation of both 4-chloromethyl- and 4-bromomethyl-benzoic acid to 4-formylbenzoic acid via hydroxylation of the α‑carbon. The 4-chloromethyl substrate bound in the enzyme active site in a similar manner to 4-ethylbenzoic acid. This places the benzylic α‑carbon hydrogens in an unfavorable position for abstraction indicating a degree of substrate mobility must be possible within the active site. CYP199A4 catalyzed oxidations of 4-(2'-haloethyl)benzoic acids yielding α-hydroxylation and desaturation metabolites. The α-hydroxylation product was the major metabolite. The desaturation pathway is significantly disfavored compared to 4-ethylbenzoic acid. This may be due to the electron-withdrawing halogen atom or a different positioning of the substrate within the active site. The latter was demonstrated by the X-ray crystal structures of CYP199A4 with these substrates. Overall, the presence of a halogen atom positioned close to the heme iron can alter the binding orientation and outcomes of enzyme-catalyzed oxidation.

Investigating the Active Oxidants Involved in Cytochrome P450 Catalyzed Sulfoxidation Reactions

Cytochrome P450 (CYP) heme-thiolate monooxygenases catalyze the hydroxylation of the C-H bonds of organic molecules. This reaction is initiated by a ferryl-oxo heme radical cation (Cpd I). These enzymes can also catalyze sulfoxidation reactions and the ferric-hydroperoxy complex (Cpd 0) and the Fe(III)-H₂ O₂ complex have been proposed as alternative oxidants for this transformation. To investigate this, the oxidation of 4-alkylthiobenzoic acids and 4-methoxybenzoic acid by the CYP199A4 enzyme from Rhodopseudomonas palustris HaA2 was compared using both monooxygenase and peroxygenase pathways. By examining mutants at the mechanistically important, conserved acid alcohol-pair (D251N, T252A and T252E) the relative amounts of the reactive intermediates that would form in these reactions were disturbed. Substrate binding and X-ray crystal structures helped to understand changes in the activity and enabled an attempt to evaluate whether multiple oxidants can participate in these reactions. In peroxygenase reactions the T252E mutant had higher activity towards sulfoxidation than O-demethylation but in the monooxygenase reactions with the WT enzyme the activity of both reactions was similar. The peroxygenase activity of the T252A mutant was greater for sulfoxidation reactions than the WT enzyme, which is the reverse of the activity changes observed for O-demethylation. The monooxygenase activity and coupling efficiency of sulfoxidation and oxidative demethylation were reduced by similar degrees with the T252A mutant. These observations infer that while Cpd I is required for O-dealkylation, another oxidant may contribute to sulfoxidation. Based on the activity of the CYP199A4 mutants it is proposed that this is the Fe(III)-H₂ O₂ complex which would be more abundant in the peroxide-driven reactions.

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

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

Tyrosyl Radical-Mediated Sequential Oxidative Decarboxylation of Coproporphyrinogen III through PCET: Theoretical Insights into the Mechanism of Coproheme Decarboxylase ChdC

The peroxide-dependent coproheme decarboxylase ChdC from Geobacillus stearothermophilus catalyzes two key steps in the synthesis of heme b, i.e., two sequential oxidative decarboxylations of coproporphyrinogen III (coproheme III) at propionate groups P2 and P4. In the binding site of coproheme III, P2 and P4 are anchored by different residues (Tyr144, Arg217, and Ser222 for P2 and Tyr113, Lys148, and Trp156 for P4); however, strong experimental evidence supports that the generated Tyr144 radical acts as an unique intermediary for hydrogen atom transfer (HAT) from both reactive propionates. So far, the reaction details are still unclear. Herein, we carried out quantum mechanics/molecular mechanics calculations to explore the decarboxylation mechanism of coproheme III. In our calculations, the coproheme Cpd I, Fe(IV) = O coupled to a porphyrin radical cation (por^•+) with four propionate groups, was used as a reactant model. Our calculations reveal that Tyr144 is directly involved in the decarboxylation of propionate group P2. First, the proton-coupled electron transfer (PCET) occurs from Tyr144 to P2, generating a Tyr144 radical, which then abstracts a hydrogen atom from the Cβ of P2. The β-H extraction was calculated to be the rate-limiting step of decarboxylation. It is the porphyrin radical cation (por^•+) that makes the PCET from Tyr144 to P2 to be quite easy to initiate the decarboxylation. Finally, the electron transfers from the Cβ^• through the porphyrin to the iron center, leading to the decarboxylation of P2. Importantly, the decarboxylation of P4 mediated by Lys148 was calculated to be very difficult, which suggests that after the P2 decarboxylation, the generated harderoheme III intermediate should rebind or rotate in the active site so that the propionate P4 occupies the binding site of P2, and Tyr144 again mediates the decarboxylation of P4. Thus, our calculations support the fact that Tyr144 is responsible for the decarboxylation of both P2 and P4.

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.

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

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

Structural insights into the role of the acid-alcohol pair of residues required for dioxygen activation in cytochrome P450 enzymes

The cytochrome P450 heme monooxygenases commonly use an acid-alcohol pair of residues, within the I-helix, to activate iron-bound dioxygen. This work aims to clarify conflicting reports on the importance of the alcohol functionality in this process. Mutants of the P450, CYP199A4 (CYP199A4_D251N and CYP199A4_T252A), were prepared, characterised and their crystal structures were solved. The acid residue of CYP199A4 is not part of a salt bridge network, a key feature of paradigmatic model system P450cam. Instead, there is a direct proton delivery network, via a chain of water molecules, extending to the surface. Nevertheless, CYP199A4_D251N dramatically reduced the activity of the enzyme consistent with a role in proton delivery. CYP199A4_T252A decreased the coupling efficiency of the enzyme with a concomitant increase in the hydrogen peroxide uncoupling pathway. However, the effect of this mutation was much less pronounced than reported with P450cam. Its crystal structures revealed fewer changes at the I-helix, compared to the P450cam system. The structural changes observed within the I-helix of P450cam during oxygen activation do not seem to be required in this P450. These differences are due to the presence of a second threonine residue at position 253, which is absent in P450cam. This threonine forms part of the hydrogen bonding network, resulting in subtle structural changes and is also present across the majority of the P450 superfamily. Overall, the results suggest that while the acid-alcohol pair is important for dioxygen activation this process and the method of proton delivery can differ across P450s.Graphic abstract.

Biophysical Techniques for Distinguishing Ligand Binding Modes in Cytochrome P450 Monooxygenases

The cytochrome P450 superfamily of heme monooxygenases catalyzes important chemical reactions across nature. The changes in the optical spectra of these enzymes, induced by the addition of substrates or inhibitors, are critical for assessing how these molecules bind to the P450, enhancing or inhibiting the catalytic cycle. Here we use the bacterial CYP199A4 enzyme (Uniprot entry Q2IUO2), from Rhodopseudomonas palustris HaA2, and a range of substituted benzoic acids to investigate different binding modes. 4-Methoxybenzoic acid elicits an archetypal type I spectral response due to a ≥95% switch from the low- to high-spin state with concomitant dissociation of the sixth aqua ligand. 4-(Pyridin-3-yl)- and 4-(pyridin-2-yl)benzoic acid induced different type II ultraviolet-visible (UV-vis) spectral responses in CYP199A4. The former induced a greater red shift in the Soret wavelength (424 nm vs 422 nm) along with a larger overall absorbance change and other differences in the α-, β-, and δ-bands. There were also variations in the ferrous UV-vis spectra of these two substrate-bound forms with a spectrum indicative of Fe-N bond formation with 4-(pyridin-3-yl)benzoic acid. The crystal structures of CYP199A4, with the pyridinyl compounds bound, revealed that while the nitrogen of 4-(pyridin-3-yl)benzoic acid is coordinated to the heme, with 4-(pyridin-2-yl)benzoic acid an aqua ligand remains. Continuous wave and pulse electron paramagnetic resonance data in frozen solution revealed that the substrates are bound in the active site in a form consistent with the crystal structures. The redox potential of each CYP199A4-substrate combination was measured, allowing correlation among binding modes, spectroscopic properties, and the observed biochemical activity.

P450-BM3-Catalyzed Sulfoxidation versus Hydroxylation: A Common or Two Different Catalytically Active Species?

While the mechanism of the P450-catalyzed oxidative hydroxylation of organic compounds has been studied in detail for many years, less is known about sulfoxidation. Depending upon the structure of the respective substrate, heme-Fe═O (Cpd I), heme-Fe(III)-OOH (Cpd 0), and heme-Fe(III)-H2O2 (protonated Cpd 0) have been proposed as reactive intermediates. In the present study, we consider the transformation of isosteric substrates via sulfoxidation and oxidative hydroxylation, respectively, catalyzed by regio- and enantioselective mutants of P450-BM3 which were constructed by directed evolution. 1-Thiochromanone and 1-tetralone were used as the isosteric substrates because, unlike previous studies involving fully flexible compounds such as thia-fatty acids and fatty acids, respectively, these compounds are rigid and cannot occur in a multitude of different conformations and binding modes in the large P450-BM3 binding pocket. The experimental results comprising activity and regio- and enantioselectivity, flanked by molecular dynamics computations within a time scale of 300 ns and QM/MM calculations of transition-state energies, unequivocally show that heme-Fe═O (Cpd I) is the common catalytically active intermediate in both sulfoxidation and oxidative hydroxylation.

Quantum-Mechanical/Molecular-Mechanical Studies of CYP11A1-Catalyzed Biosynthesis of Pregnenolone from Cholesterol Reveal a C-C Bond Cleavage Reaction That Occurs by a Compound I-Mediated Electron Transfer

We explore here a long-standing mechanistic question by using quantum-mechanical/molecular-mechanical (QM/MM) methodology. The question concerns the mechanism of steroid hormone biosynthesis, whereby the P450 enzyme, CYP11A1, catalyzes the C20-C22 bond-cleavage in the 20,22-hydroxylated cholesterol, 20R,22R-DiOHCH, leading to pregnenolone, which is critical for the subsequent production of all steroid hormones. This is an unusual feat whereby the P450 enzyme breaks two O-H bonds and one C-C bond, while making two C═O bonds. How does the enzyme perform such a complex and highly energy-demanding reaction? Our computational results rule out the previously proposed Compound I (Cpd I) electrophilic attack mechanism via the formation of a peroxide intermediate as well as the H-abstraction-mediated C-C cleavage mechanism. Notably, oxygen-rebound cannot transpire, in spite of the fact that the classical active species, Cpd I, participates in the catalytic process. Our findings reveal a mechanism whereby C-C bond cleavage is mediated by an electron transfer from the C22-O^--deprotonated substrate to Cpd I. As such, our QM/MM calculations demonstrate that Cpd I acts as an electron sink that facilitates the C-C bond cleavage.

Investigation of the requirements for efficient and selective cytochrome P450 monooxygenase catalysis across different reactions

The cytochrome P450 metalloenzyme (CYP) CYP199A4 from Rhodopseudomonas palustris HaA2 catalyzes the highly efficient oxidation of para-substituted benzoic acids. Here we determined crystal structures of CYP199A4, and the binding and turnover parameters, with different meta-substituted benzoic acids in order to establish which criteria are important for efficient catalysis. When compared to the para isomers, the meta-substituted benzoic acids were less efficiently oxidized. For example, 3-formylbenzoic acid was oxidized with lower activity than the equivalent para isomer and 3-methoxybenzoic acid did not undergo O-demethylation by CYP199A4. The structural data highlighted that the meta-substituted benzoic acids bound in the enzyme active site in a modified position with incomplete loss of the distal water ligand of the heme moiety. However, for both sets of isomers the meta- or para-substituent pointed towards, and was in close proximity, to the heme iron. The absence of oxidation activity with 3-methoxybenzoic acid was assigned to the observation that the CH bonds of this molecule point away from the heme iron. In contrast, in the para isomer they are in an ideal location for abstraction. These findings were confirmed by using the bulkier 3-ethoxybenzoic acid as a substrate which removed the water ligand and reoriented the meta-substituent so that the methylene hydrogens pointed towards the heme, enabling more efficient oxidation. Overall we show relatively small changes in substrate structure and position in the active site can have a dramatic effect on the activity.

Synthetic Fe/Cu Complexes: Toward Understanding Heme-Copper Oxidase Structure and Function

Heme-copper oxidases (HCOs) are terminal enzymes on the mitochondrial or bacterial respiratory electron transport chain, which utilize a unique heterobinuclear active site to catalyze the 4H+/4e- reduction of dioxygen to water. This process involves a proton-coupled electron transfer (PCET) from a tyrosine (phenolic) residue and additional redox events coupled to transmembrane proton pumping and ATP synthesis. Given that HCOs are large, complex, membrane-bound enzymes, bioinspired synthetic model chemistry is a promising approach to better understand heme-Cu-mediated dioxygen reduction, including the details of proton and electron movements. This review encompasses important aspects of heme-O2 and copper-O2 (bio)chemistries as they relate to the design and interpretation of small molecule model systems and provides perspectives from fundamental coordination chemistry, which can be applied to the understanding of HCO activity. We focus on recent advancements from studies of heme-Cu models, evaluating experimental and computational results, which highlight important fundamental structure-function relationships. Finally, we provide an outlook for future potential contributions from synthetic inorganic chemistry and discuss their implications with relevance to biological O2-reduction.

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.

How Are Substrate Binding and Catalysis Affected by Mutating Glu127 and Arg161 in Prolyl-4-hydroxylase? A QM/MM and MD Study

Prolyl-4-hydroxylase is a vital enzyme for human physiology involved in the biosynthesis of 4-hydroxyproline, an essential component for collagen formation. The enzyme performs a unique stereo- and regioselective hydroxylation at the C4 position of proline despite the fact that the C5 hydrogen atoms should be thermodynamically easier to abstract. To gain insight into the mechanism and find the origin of this regioselectivity, we have done a quantum mechanics/molecular mechanics (QM/MM) study on wildtype and mutant structures. In a previous study (Timmins et al., 2017) we identified several active site residues critical for substrate binding and positioning. In particular, the Glu127 and Arg161 were shown to form multiple hydrogen bonding and ion-dipole interactions with substrate and could thereby affect the regio- and stereoselectivity of the reaction. In this work, we decided to test that hypothesis and report a QM/MM and molecular dynamics (MD) study on prolyl-4-hydroxylase and several active site mutants where Glu127 or Arg161 are mutated for Asp, Gln, or Lys. Thus, the R161D and R161Q mutants give very high barriers for hydrogen atom abstraction from any proline C-H bond and therefore will be inactive. The R161K mutant, by contrast, sees the regio- and stereoselectivity of the reaction change but still is expected to hydroxylate proline at room temperature. By contrast, the Glu127 mutants E127D and E127Q show possible changes in regioselectivity with the former being more probable to react compared to the latter.

MD simulations and QM/MM calculations show that single-site mutations of cytochrome P450BM3 alter the active site's complexity and the chemoselectivity of oxidation without changing the active species

A new water channel appears in the T268A mutant of P450_BM3 and plays a role in the enzyme’s chemoselectivity.

It is a long-standing mechanistic consensus that the mutation of the proton-shuttle mediator Threonine (T) in Cytochrome P450 enzymes severs the water channel and thereby quenches the formation of the active species: the high-valent iron(iv)-oxo porphyrin π-cation radical species, compound I (Cpd I). Using MD simulations and hybrid QM/MM calculations of P450_BM3 we demonstrate that this is not the case. Thus, while the original water channel is disrupted in the T268A mutant of the enzyme, a new channel is formed that generates Cpd I. With this new understanding, we address the puzzling regiochemical and kinetic-isotope effect (KIE) results (Volz et al., J. Am. Chem. Soc., 2002, 124, 9724–9725) on the sulfoxidation and N-dealkylation of dimethyl-(4-methylsulfanyl-phenyl)-amine by wild type (WT) P450_BM3 and its T268A vs. F87A mutants. We show that the observed variable ratio of S/Me oxidation for these enzymes, vis-à-vis the constant KIE, originates from Cpd I being the sole oxidant. Thus, while the conserved KIE probes the conserved nature of the transition state, the variable regiochemical S/Me ratio reflects the active-site reorganization in the mutants: the shifted location of the new water channel in T268A tightens the binding of the S-end by Cpd I and increases the S/Me ratio, whereas the absence of π-interaction with the S-end in F87A creates a looser binding that lowers the S/Me ratio. Our results match the experimental findings. As such, this study sheds light on puzzling experimental results, and may shift a central paradigm in P450 research. The broader implication on enzymatic research is that a single-site mutation is not a localised alteration but one that may lead to a profound change in the active site, sufficiently so as to change the chemoselectivity of catalyzed reactions.

Substrate Sulfoxidation by an Iron(IV)-Oxo Complex: Benchmarking Computationally Calculated Barrier Heights to Experiment

High-valent metal-oxo oxidants are common reactive species in synthetic catalysts as well as heme and nonheme iron enzymes. In general, they efficiently react with substrates through oxygen atom transfer, and for a number of cases, experimental rate constants have been determined. However, because these rate constants are generally measured in a polar solution, it has been found difficult to find computational methodologies to reproduce experimental trends and reactivities. In this work, we present a detailed computational study into para-substituted thioanisole sulfoxidation by a nonheme iron(IV)-oxo complex. A range of density functional theory methods and basis sets has been tested for their suitability to describe the reaction mechanism and compared with experimentally obtained free energies of activation. It is found that the enthalpy of activation is reproduced well, but all methods overestimate the entropy of activation by about 50%, for which we recommend a correction factor. The effect of solvent and dispersion on the barrier heights is explored both at the single-point level and also through inclusion in geometry optimizations, and particularly, solvent is seen as highly beneficial to reproduce experimental free energies of activation. Interestingly, in general, experimental trends and Hammett plots are reproduced well with almost all methods and procedures, and only a systematic error seems to apply for these chemical systems. Very good agreement between experiment and theory is found for a number of different methods, including B3LYP and PBE0, and procedures that are highlighted in the paper.

Predicting Regioselectivity and Lability of Cytochrome P450 Metabolism Using Quantum Mechanical Simulations

We describe methods for predicting cytochrome P450 (CYP) metabolism incorporating both pathway-specific reactivity and isoform-specific accessibility considerations. Semiempirical quantum mechanical (QM) simulations, parametrized using experimental data and ab initio calculations, estimate the reactivity of each potential site of metabolism (SOM) in the context of the whole molecule. Ligand-based models, trained using high-quality regioselectivity data, correct for orientation and steric effects of the different CYP isoform binding pockets. The resulting models identify a SOM in the top 2 predictions for between 82% and 91% of compounds in independent test sets across seven CYP isoforms. In addition to predicting the relative proportion of metabolite formation at each site, these methods estimate the activation energy at each site, from which additional information can be derived regarding their lability in absolute terms. We illustrate how this can guide the design of compounds to overcome issues with rapid CYP metabolism.

Oxidation of Aryl Diphenylmethyl Sulfides Promoted by a Nonheme Iron(IV)-Oxo Complex: Evidence for an Electron Transfer-Oxygen Transfer Mechanism

The oxidation of a series of aryl diphenylmethyl sulfides (4-X-C6H4SCH(C6H5)2, where X = OCH3 (1), X = CH3 (2), X = H (3), and X = CF3 (4)) promoted by the nonheme iron(IV)-oxo complex [(N4Py)Fe(IV)═O](2+) occurs by an electron transfer-oxygen transfer (ET-OT) mechanism as supported by the observation of products (diphenylmethanol, benzophenone, and diaryl disulfides) deriving from α-C-S and α-C-H fragmentation of radical cations 1(+•)-4(+•), formed besides the S-oxidation products (aryl diphenylmethyl sulfoxides). The fragmentation/S-oxidation product ratios regularly increase through a decrease in the electron-donating power of the aryl substituents, that is, by increasing the fragmentation rate constants of the radical cations as indicated by a laser flash photolysis (LFP) study of the photochemical oxidation of 1-4 carried out in the presence of N-methoxyphenanthridinium hexafluorophosphate (MeOP(+)PF6(-)).

Toxicity Originating from Thiophene Containing Drugs: Exploring the Mechanism using Quantum Chemical Methods

Drug metabolism of thiophene containing substrates by cytochrome P450s (CYP450) leads to toxic side effects, for example, nephrotoxicity (suprofen, ticlopidine), hepatotoxicity (tienilic acid), thrombotic thrombocytopenic purpura (clopidogrel), and aplastic anemia (ticlopidine). The origin of toxicity in these cases has been attributed to two different CYP450 mediated metabolic reactions: S-oxidation and epoxidation. In this work, the molecular level details of the bioinorganic chemistry associated with the generation of these competitive reactions are reported. Density functional theory was utilized (i) to explore the molecular mechanism for S-oxidation and epoxidation using the radical cationic center Cpd I [(iron(IV)-oxo-heme porphine system with SH(-) as the axial ligand, to mimic CYP450s] as the model oxidant, (ii) to establish the 3D structures of the reactants, transition states, and products on both the metabolic pathways, and (iii) to examine the potential energy (PE) profile for both the pathways to determine the energetically preferred toxic metabolite formation. The energy barrier required for S-oxidation was observed to be 14.75 kcal/mol as compared to that of the epoxidation reaction (13.23 kcal/mol) on the doublet PE surface of Cpd I. The formation of the epoxide metabolite was found to be highly exothermic (-23.24 kcal/mol), as compared to S-oxidation (-8.08 kcal/mol). Hence, on a relative scale the epoxidation process was observed to be thermodynamically and kinetically more favorable. The energy profiles associated with the reactions of the S-oxide and epoxide toxic metabolites were also explored. This study helps in understanding the CYP450-catalyzed toxic reactions of drugs containing the thiophene ring at the atomic level.

Quantum mechanical/molecular mechanical calculated reactivity networks reveal how cytochrome P450cam and Its T252A mutant select their oxidation pathways

Quantum mechanical/molecular mechanical calculations address the longstanding-question of a "second oxidant" in P450 enzymes wherein the proton-shuttle, which leads to formation of the "primary-oxidant" Compound I (Cpd I), was severed by mutating the crucial residue (in P450cam: Threonine-252-to-Alanine, hence T252A). Investigating the oxidant candidates Cpd I, ferric hydroperoxide, and ferric hydrogen peroxide (Fe(III)(O2H2)), and their reactions, generates reactivity networks which enable us to rule out a "second oxidant" and at the same time identify an additional coupling pathway that is responsible for the epoxidation of 5-methylenylcamphor by the T252A mutant. In this "second-coupling pathway", the reaction starts with the Fe(III)(O2H2) intermediate, which transforms to Cpd I via a O-O homolysis/H-abstraction mechanism. The persistence of Fe(III)(O2H2) and its oxidative reactivity are shown to be determined by interplay of substrate and protein. The substrate 5-methylenylcamphor prevents H2O2 release, while the protein controls the Fe(III)(O2H2) conversion to Cpd I by nailing-through hydrogen-bonding interactions-the conformation of the HO(•) radical produced during O-O homolysis. This conformation prevents HO(•) attack on the porphyrin's meso position, as in heme oxygenase, and prefers H-abstraction from Fe(IV)OH thereby generating H2O + Cpd I. Cpd I then performs substrate oxidations. Camphor cannot prevent H2O2 release and hence the T252A mutant does not oxidize camphor. This "second pathway" transpires also during H2O2 shunting of the cycle of wild-type P450cam, where the additional hydrogen-bonding with Thr252 prevents H2O2 release, and contributes to a successful Cpd I formation. The present results lead to a revised catalytic cycle of Cytochrome P450cam.

ortho-Hydroxylation of aromatic acids by a non-heme Fe(V)=O species: how important is the ligand design?

There is a growing interest in probing the mechanism of catalytic transformations effected by non-heme iron-oxo complexes as these reactions set a platform for understanding the relevant enzymatic reactions. The ortho-hydroxylation of aromatic compounds is one such reaction catalysed by iron-oxo complexes. Experimentally [Fe(II)(BPMEN)(CH3CN)2](2+) (1) and [Fe(II)(TPA)(CH3CN)2](2+) (2) (where TPA = tris(2-pyridylmethyl)amine and BPMEN = N,N′-dimethyl-N,N′-bis(2-pyridylmethyl)ethane-1,2-diamine) complexes containing amino pyridine ligands along with H2O2 are employed to carry out these transformations where complex 1 is found to be more reactive than complex 2. Herein, using density functional methods employing B3LYP and dispersion corrected B3LYP (B3LYP-D) functionals, we have explored the mechanism of this reaction to reason out the importance of ligand design in fine-tuning the reactivity of such catalytic transformations. Dispersion corrected B3LYP is found to be superior to B3LYP in predicting the correct ground state of these species and also yields lower barrier heights than the B3LYP functional. Starting the reaction from the Fe(III)–OOH species, both homolytic and heterolytic cleavage of the O···O bond is explored leading to the formation of the transient Fe(IV)=O and Fe(V)=O species. For both the ligand systems, heterolytic cleavage was energetically preferable and our calculations suggest that both the reactions are catalyzed by an elusive high-valent Fe(V)=O species. The Fe(V)=O species undergoes the reaction via an electrophilic attack of the benzene ring to effect the ortho-hydroxylation reaction. The reactivity pattern observed for 1 and 2 are reflected in the computed barrier heights for the ortho-hydroxylation reaction. Electronic structure analysis reveals that the difference in reactivity between the ligand architectures described in complex 1 and 2 arise due to orientation of the pyridine ring(s) parallel or perpendicular to the Fe(V)=O bond. The parallel orientation of the pyridine ring is found to mix with the (πFe(dyz)–O(py))* orbital of the Fe-oxo bond leading to a reduction in the electrophilicity of the ferryl oxygen atom. Our calculations highlight the importance of ligand design in this chemistry and suggest that this concept can be used to (i) stabilize high-valent intermediates which can be trapped and thoroughly characterized (ii) enhance the reactivity and efficiency of the oxidants by increasing the electrophilicity of the ferryl oxygen containing FeVO species. Our computed results are in general agreement with the experimental results.

Monooxygenase, peroxidase and peroxygenase properties and reaction mechanisms of cytochrome P450 enzymes

This review examines the monooxygenase, peroxidase and peroxygenase properties and reaction mechanisms of cytochrome P450 (CYP) enzymes in bacterial, archaeal and mammalian systems. CYP enzymes catalyze monooxygenation reactions by inserting one oxygen atom from O2 into an enormous number and variety of substrates. The catalytic versatility of CYP stems from its ability to functionalize unactivated carbon-hydrogen (C-H) bonds of substrates through monooxygenation. The oxidative prowess of CYP in catalyzing monooxygenation reactions is attributed primarily to a porphyrin π radical ferryl intermediate known as Compound I (CpdI) (Por•+FeIV=O), or its ferryl radical resonance form (FeIV-O•). CYP-mediated hydroxylations occur via a consensus H atom abstraction/oxygen rebound mechanism involving an initial abstraction by CpdI of a H atom from the substrate, generating a highly-reactive protonated Compound II (CpdII) intermediate (FeIV-OH) and a carbon-centered alkyl radical that rebounds onto the ferryl hydroxyl moiety to yield the hydroxylated substrate. CYP enzymes utilize hydroperoxides, peracids, perborate, percarbonate, periodate, chlorite, iodosobenzene and N-oxides as surrogate oxygen atom donors to oxygenate substrates via the shunt pathway in the absence of NAD(P)H/O2 and reduction-oxidation (redox) auxiliary proteins. It has been difficult to isolate the historically elusive CpdI intermediate in the native NAD(P)H/O2-supported monooxygenase pathway and to determine its precise electronic structure and kinetic and physicochemical properties because of its high reactivity, unstable nature (t½~2 ms) and short life cycle, prompting suggestions for participation in monooxygenation reactions of alternative CYP iron-oxygen intermediates such as the ferric-peroxo anion species (FeIII-OO-), ferric-hydroperoxo species (FeIII-OOH) and FeIII-(H2O2) complex.

Mechanistic insight into peroxo-shunt formation of biomimetic models for compound II, their reactivity toward organic substrates, and the influence of N-methylimidazole axial ligation

High-valent iron-oxo species have been invoked as reactive intermediates in catalytic cycles of heme and nonheme enzymes. The studies presented herein are devoted to the formation of compound II model complexes, with the application of a water soluble (TMPS)Fe(III)(OH) porphyrin ([meso-tetrakis(2,4,6-trimethyl-3-sulfonatophenyl)porphinato]iron(III) hydroxide) and hydrogen peroxide as oxidant, and their reactivity toward selected organic substrates. The kinetics of the reaction of H2O2 with (TMPS)Fe(III)(OH) was studied as a function of temperature and pressure. The negative values of the activation entropy and activation volume for the formation of (TMPS)Fe(IV)=O(OH) point to the overall associative nature of the process. A pH-dependence study on the formation of (TMPS)Fe(IV)=O(OH) revealed a very high reactivity of OOH(-) toward (TMPS)Fe(III)(OH) in comparison to H2O2. The influence of N-methylimidazole (N-MeIm) ligation on both the formation of iron(IV)-oxo species and their oxidising properties in the reactions with 4-methoxybenzyl alcohol or 4-methoxybenzaldehyde, was investigated in detail. Combined experimental and theoretical studies revealed that among the studied complexes, (TMPS)Fe(III)(H2O)(N-MeIm) is highly reactive toward H2O2 to form the iron(IV)-oxo species, (TMPS)Fe(IV)=O(N-MeIm). The latter species can also be formed in the reaction of (TMPS)Fe(III)(N-MeIm)2 with H2O2 or in the direct reaction of (TMPS)Fe(IV)=O(OH) with N-MeIm. Interestingly, the kinetic studies involving substrate oxidation by (TMPS)Fe(IV)=O(OH) and (TMPS)Fe(IV)=O(N-MeIm) do not display a pronounced effect of the N-MeIm axial ligand on the reactivity of the compound II mimic in comparison to the OH(-) substituted analogue. Similarly, DFT computations revealed that the presence of an axial ligand (OH(-) or N-MeIm) in the trans position to the oxo group in the iron(IV)-oxo species does not significantly affect the activation barriers calculated for C-H dehydrogenation of the selected organic substrates.

Use of density functional theory in drug metabolism studies

INTRODUCTION

The cytochrome P450 enzymes (CYPs) metabolize many drug compounds. They catalyze a wide variety of reactions, and potentially, a large number of different metabolites can be generated. Density functional theory (DFT) has, over the past decade, been shown to be a powerful tool to rationalize and predict the possible metabolites generated by the CYPs as well as other drug-metabolizing enzymes.

AREAS COVERED

We review applications of DFT on reactions performed by the CYPs and other drug-metabolizing enzymes able to perform oxidation reactions, with an emphasis on predicting which metabolites are produced. We also cover calculations of binding energies for complexes in which the ligands interact directly with the heme iron atom.

EXPERT OPINION

DFT is a useful tool for prediction of the site of metabolism. The use of small models of the enzymes work surprisingly well for most CYP isoforms. This is probably due to the fact that the binding of the substrates is not the major determinant. When binding of the substrate plays a significant role, the well-known issue of determining the free energy of binding is the challenge. How approaches taking the protein environment into account, like docking, MD and QM/MM, can be used are discussed.

Prediction of cytochrome P450 xenobiotic metabolism: tethered docking and reactivity derived from ligand molecular orbital analysis

Metabolism of xenobiotic and endogenous compounds is frequently complex, not completely elucidated, and therefore often ambiguous. The prediction of sites of metabolism (SoM) can be particularly helpful as a first step toward the identification of metabolites, a process especially relevant to drug discovery. This paper describes a reactivity approach for predicting SoM whereby reactivity is derived directly from the ground state ligand molecular orbital analysis, calculated using Density Functional Theory, using a novel implementation of the average local ionization energy. Thus each potential SoM is sampled in the context of the whole ligand, in contrast to other popular approaches where activation energies are calculated for a predefined database of molecular fragments and assigned to matching moieties in a query ligand. In addition, one of the first descriptions of molecular dynamics of cytochrome P450 (CYP) isoforms 3A4, 2D6, and 2C9 in their Compound I state is reported, and, from the representative protein structures obtained, an analysis and evaluation of various docking approaches using GOLD is performed. In particular, a covalent docking approach is described coupled with the modeling of important electrostatic interactions between CYP and ligand using spherical constraints. Combining the docking and reactivity results, obtained using standard functionality from common docking and quantum chemical applications, enables a SoM to be identified in the top 2 predictions for 75%, 80%, and 78% of the data sets for 3A4, 2D6, and 2C9, respectively, results that are accessible and competitive with other recently published prediction tools.

Density functional study on the cytochrome-mediated S-oxidation: identification of crucial reactive intermediate on the metabolic path of thiazolidinediones

S-Oxidation is an important cytochrome P450 (CYP450)-catalyzed reaction, and the structural and energetic details of this process can only be studied by using quantum chemical methods. Thiazolidinedione (TZD) ring metabolism involving initial S-oxidation leads to the generation of reactive metabolites (RMs) and subsequent toxicity forcing the withdrawal of the glitazone class of drugs, thus, the study of the biochemical pathway of TZD ring metabolism is a subject of interest. The S-oxidation of the TZD ring and the formation of the isocyanate intermediate (ISC) was implicated as a possible pathway; however, there are several questions still unanswered in this biochemical pathway. The current study focuses on the CYP450-mediated S-oxidation, fate of the sulfoxide product (TZDSO), ring cleavage to ISC, and formation of nucleophilic adducts. The process of S-oxidation was explored by using Cpd I (iron(IV)-oxo porphyrin, to mimic CYP450) at TZVP/6-311+G(d) basis set. The barriers were calculated after incorporating dispersion and solvent corrections. The metabolic conversion from TZDSO to ISC (studied at B3LYP/6-311++G(2df,3pd)//B3LYP/6-31+G(d)) required a novel protonated intermediate, TZDSOH(+). The effect of higher basis sets (6-311+G(d,p), aug-cc-pvqz) on this conversion was studied. TZDSOH(+) was observed to be more reactive and thermodynamically accessible than ISC, indicating that TZDSOH(+) is the actual reactive intermediate leading to toxicity of the TZD class of compounds.

Metabolic-intermediate complex formation with cytochrome P450: theoretical studies in elucidating the reaction pathway for the generation of reactive nitroso intermediate

Mechanism-based inhibition (MBI) of cytochrome P450 (CYP) can lead to drug-drug interactions and often to toxicity. Some aliphatic and aromatic amines can undergo biotransformation reactions to form reactive metabolites such as nitrosoalkanes, leading to MBI of CYPs. It has been proposed that the nitrosoalkanes coordinate with the heme iron, forming metabolic-intermediate complex (MIC), resulting in the quasi-irreversible inhibition of CYPs. Limited mechanistic details regarding the formation of reactive nitroso intermediate and its coordination with heme-iron have been reported. A quantum chemical analysis was performed to elucidate potential reaction pathways for the generation of nitroso intermediate and the formation of MIC. Elucidation of the energy profile along the reaction path, identification of three-dimensional structures of reactive intermediates and transition states, as well as charge and spin density analyses, were performed using the density functional B3LYP method. The study was performed using Cpd I [iron (IV-oxo] heme porphine with SH(-) as the axial ligand) to represent the catalytic domain of CYP, simulating the biotransformation process. Three pathways: (i) N-oxidation followed by proton shuttle, (ii) N-oxidation followed by 1,2-H shift, and (iii) H-abstraction followed by rebound mechanism, were studied. It was observed that the proton shuttle pathway was more favorable over the whole reaction leading to reactive nitroso intermediate. This study revealed that the MIC formation from a primary amine is a favorable exothermic process, involving eight different steps and preferably takes place on the doublet spin surface of Cpd I. The rate-determining step was identified to be the first N-oxidation of primary amine.