Molecular Dynamics and QM/MM Calculations Predict the Substrate-Induced Gating of Cytochrome P450 BM3 and the Regio- and Stereoselectivity of Fatty Acid Hydroxylation

Cytochrome P450 Enzyme Design by Constraining the Catalytic Pocket in a Diffusion Model

Although cytochrome P450 enzymes are the most versatile biocatalysts in nature, there is insufficient comprehension of the molecular mechanism underlying their functional innovation process. Here, by combining ancestral sequence reconstruction, reverse mutation assay, and progressive forward accumulation, we identified 5 founder residues in the catalytic pocket of flavone 6-hydroxylase (F6H) and proposed a "3-point fixation" model to elucidate the functional innovation mechanisms of P450s in nature. According to this design principle of catalytic pocket, we further developed a de novo diffusion model (P450Diffusion) to generate artificial P450s. Ultimately, among the 17 non-natural P450s we generated, 10 designs exhibited significant F6H activity and 6 exhibited a 1.3- to 3.5-fold increase in catalytic capacity compared to the natural CYP706X1. This work not only explores the design principle of catalytic pockets of P450s, but also provides an insight into the artificial design of P450 enzymes with desired functions.

Mechanism of 7H-Dibenzo[c,g]carbazole metabolism in cytochrome P450 1A1: Insights from computational studies

7H-Dibenzo[c,g]carbazole (DBC) is a prevalent environmental contaminant that induces tumorigenesis in several experimental animals. Recently, it has been utilized to develop high-performance solar cells and organic phosphorescent materials. It is imperative to strengthen investigations of DBC metabolism to understand its potential risks to human health. In this study, human CYP1A1 was employed as the metabolic enzyme to investigate the metabolic mechanism of DBC by molecular docking, molecular dynamics (MD) simulation, and quantum mechanical (QM) calculation. The results indicate that DBC binds to CYP1A1 in two modes (mode 1 and mode 2) mainly through nonpolar solvation energies (ΔGnonpolar). The formation of the two binding modes is attributed to the anchoring effect of the hydrogen bond formed by DBC with Asp320 (mode 1) or Ser116 (mode 2). Mode 1 is a "reactive" conformation, while mode 2 is not considered a "reactive" conformation. C5 is identified as the dominant site, and the pyrrole nitrogen cannot participate in the metabolism. DBC is metabolized mainly by a distinct electrophilic addition-rearrangement mechanism, with an energy barrier of 21.74 kcal/mol. The results provide meaningful insights into the biometabolic process of DBC and contribute to understanding its environmental effects and health risks.

Computational Characterization of the Reactivity of Compound I in Unspecific Peroxygenases

Unspecific peroxygenases (UPOs) are emerging as promising biocatalysts for selective oxyfunctionalization of unactivated C-H bonds. However, their potential in large-scale synthesis is currently constrained by suboptimal chemical selectivity. Improving the selectivity of UPOs requires a deep understanding of the molecular basis of their catalysis. Recent molecular simulations have sought to unravel UPO's selectivity and inform their design principles. However, most of these studies focused on substrate-binding poses. Few researchers have investigated how the reactivity of CpdI, the principal oxidizing intermediate in the catalytic cycle, influences selectivity in a realistic protein environment. Moreover, the influence of protein electrostatics on the reaction kinetics of CpdI has also been largely overlooked. To bridge this gap, we used multiscale simulations to interpret the regio- and enantioselective hydroxylation of the n-heptane substrate catalyzed by Agrocybe aegerita UPO (AaeUPO). We comprehensively characterized the energetics and kinetics of the hydrogen atom-transfer (HAT) step, initiated by CpdI, and the subsequent oxygen rebound step forming the product. Notably, our approach involved both free energy and potential energy evaluations in a quantum mechanics/molecular mechanics (QM/MM) setting, mitigating the dependence of results on the choice of initial conditions. These calculations illuminate the thermodynamics and kinetics of the HAT and oxygen rebound steps. Our findings highlight that both the conformational selection and the distinct chemical reactivity of different substrate hydrogen atoms together dictate the regio- and enantio-selectivity. Building on our previous study of CpdI's formation in AaeUPO, our results indicate that the HAT step is the rate-limiting step in the overall catalytic cycle. The subsequent oxygen rebound step is swift and retains the selectivity determined by the HAT step. We also pinpointed several polar and charged amino acid residues whose electrostatic potentials considerably influence the reaction barrier of the HAT step. Notably, the Glu196 residue is pivotal for both the CpdI's formation and participation in the HAT step. Our research offers in-depth insights into the catalytic cycle of AaeUPO, which will be instrumental in the rational design of UPOs with enhanced properties.

Computationally guided bioengineering of the active site, substrate access pathway, and water channels of thermostable cytochrome P450, CYP175A1, for catalyzing the alkane hydroxylation reaction

Understanding structure-function relationships in proteins is pivotal in their development as industrial biocatalysts. In this regard, rational engineering of protein active site access pathways and various tunnels and channels plays a central role in designing competent enzymes with high stability and enhanced efficiency. Here, we report the rational evolution of a thermostable cytochrome P450, CYP175A1, to catalyze the C-H activation reaction of longer-chain alkanes. A strategy combining computational tools with experiments has shown that the substrate scope and enzymatic activity can be enhanced by rational engineering of certain important channels such as the substrate entry and water channels along with the active site of the enzyme. The evolved enzymes showed an improved catalytic rate for hexadecane hydroxylation with high regioselectivity. The Q67L/Y68F mutation showed binding of the substrate in the active site, water channel mutation L80F/V220T showed improved catalytic activity through the peroxide shunt pathway and substrate entry channel mutation W269F/I270A showed better substrate accessibility to the active pocket. All-atom MD simulations provided the rationale for the inactivity of the wild-type CYP175A1 for hexadecane hydroxylation and predicted the above hot-spot residues to enhance the activity. The reaction mechanism was studied by QM/MM calculations for enzyme-substrate complexes and reaction intermediates. Detailed thermal and thermodynamic stability of all the mutants were analyzed and the results showed that the evolved enzymes were thermally stable. The present strategy showed promising results, and insights gained from this work can be applied to the general enzymatic system to expand substrate scope and improve catalytic activity.

A Simple Access to γ- and ε-Keto Arenes via Enzymatic Divergent C─H Bond Oxyfunctionalization

Performing divergent C─H bond functionalization on molecules with multiple reaction sites is a significant challenge in organic chemistry. Biocatalytic oxyfunctionalization reactions of these compounds to the corresponding ketones/aldehydes are typically hindered by selectivity issues. To address these challenges, the catalytic performance of oxidoreductases is explored. The results show that combining the peroxygenase-catalyzed propargylic C─H bond oxidation with the Old Yellow Enzyme-catalyzed reduction of conjugated C─C triple bonds in one-pot enables the regio- and chemoselective oxyfunctionalization of sp3 C─H bonds that are distant from benzylic sites. This enzymatic approach yielded a variety of γ-keto arenes with diverse structural and electronic properties in yields of up to 99% and regioselectivity of 100%, which are difficult to achieve using other chemocatalysis and enzymes. By adjusting the C─C triple bond, the carbonyl group's position can be further tuned to yield ε-keto arenes. This enzymatic approach can be combined with other biocatalysts to establish new synthetic pathways for accessing various challenging divergent C─H bond functionalization reactions.

Understanding the Multifaceted Mechanism of Compound I Formation in Unspecific Peroxygenases through Multiscale Simulations

Unspecific peroxygenases (UPOs) can selectively oxyfunctionalize unactivated hydrocarbons by using peroxides under mild conditions. They circumvent the oxygen dilemma faced by cytochrome P450s and exhibit greater stability than the latter. As such, they hold great potential for industrial applications. A thorough understanding of their catalysis is needed to improve their catalytic performance. However, it remains elusive how UPOs effectively convert peroxide to Compound I (CpdI), the principal oxidizing intermediate in the catalytic cycle. Previous computational studies of this process primarily focused on heme peroxidases and P450s, which have significant differences in the active site from UPOs. Additionally, the roles of peroxide unbinding in the kinetics of CpdI formation, which is essential for interpreting existing experiments, have been understudied. Moreover, there has been a lack of free energy characterizations with explicit sampling of protein and hydration dynamics, which is critical for understanding the thermodynamics of the proton transport (PT) events involved in CpdI formation. To bridge these gaps, we employed multiscale simulations to comprehensively characterize the CpdI formation in wild-type UPO from Agrocybe aegerita (AaeUPO). Extensive free energy and potential energy calculations were performed in a quantum mechanics/molecular mechanics setting. Our results indicate that substrate-binding dehydrates the active site, impeding the PT from H2O2 to a nearby catalytic base (Glu196). Furthermore, the PT is coupled with considerable hydrogen bond network rearrangements near the active site, facilitating subsequent O-O bond cleavage. Finally, large unbinding free energy barriers kinetically stabilize H2O2 at the active site. These findings reveal a delicate balance among PT, hydration dynamics, hydrogen bond rearrangement, and cosubstrate unbinding, which collectively enable efficient CpdI formation. Our simulation results are consistent with kinetic measurements and offer new insights into the CpdI formation mechanism at atomic-level details, which can potentially aid the design of next-generation biocatalysts for sustainable chemical transformations of feedstocks.

Enantioselective metabolism of novel chiral insecticide Paichongding by human cytochrome P450 3A4: A computational insight

As a novel chiral neonicotinoid insecticide, Paichongding (IPP) has been widely applied in agriculture due to its excellent insecticidal activity. However, the enantioselective metabolism of IPP stereoisomers (5R7R-IPP, 5S7S-IPP, 5R7S-IPP, and 5S7R-IPP) mediated by enzymes in non-target organisms, especially the cytochrome P450s (CYPs), remains unknown. To address this knowledge gap, we developed an integrated computational framework to elucidate the binding interactions and enantioselective metabolism of IPP stereoisomers in human CYP3A4. The results reveal that 5R7R-IPP shows much stronger binding affinity to CYP3A4 than 5S7S-IPP, while enantiomers 5R7S-IPP and 5S7R-IPP have no essential difference in their binding potential, owing to their specific interactions with key CYP3A4 residues. Although enantiomers 5R7R-IPP and 5S7S-IPP feature distinct binding modes resulting from the chiral differences, their transformation activities are slightly different, with C5 and C13 being the primary metabolic sites, respectively. In contrast, CYP3A4 preferably metabolizes 5R7S-IPP over 5S7R-IPP. The metabolism of epimers 5R7R-IPP and 5R7S-IPP share C5-hydroxylation routes due to the conserved 5R-conformaitons, but differ with the transformation routes at C11/C13 and C3 sites. The 7R-chirality of 5S7R-IPP significantly reduces the metabolic potency compared to 5S7S-IPP. CYP3A4-catalyzed hydroxylation and desaturation of IPP stereoisomers generate various chiral metabolites, with C5- and C13-hydroxyIPPs further transforming into depropylated products. Furthermore, the toxicity assessment reveals that IPP, along with the majority of its hydroxylated, desaturated, and depropylated metabolites, can potentially induce adverse effects on human health, specifically hepatotoxicity, respiratory toxicity, and carcinogenicity. This study provides valuable insights into the enantioselective fate of chiral IPP metabolism by CYP3A4, and the identified metabolites can serve as potential biomarkers for monitoring IPP exposure and associated health risk in human body.

Catalysis by KDM6 Histone Demethylases - A Synergy between the Non-Heme Iron(II) Center, Second Coordination Sphere, and Long-Range Interactions

KDM6A (UTX) and KDM6B (JMJD3) are human non-heme Fe(II) and 2-oxoglutarate (2OG) dependent JmjC oxygenases that catalyze the demethylation of trimethylated lysine 27 in the N-terminal tail of histone H3, a post-translational modification that regulates transcription. A Combined Quantum Mechanics/ Molecular Mechanics (QM/MM) and Molecular Dynamics (MD) study on the catalytic mechanism of KDM6A/B reveals that the transition state for the rate-limiting hydrogen atom transfer (HAT) reaction in KDM6A catalysis is stabilized by polar (Asn217) and aromatic (Trp369)/non-polar (Pro274) residues in contrast to KDM4, KDM6B and KDM7 demethylases where charged residues (Glu, Arg, Asp) are involved. KDM6A employs both σ- and π-electron transfer pathways for HAT, whereas KDM6B employs the σ-electron pathway. Differences in hydrogen bonding of the Fe-chelating Glu252(KDM6B) contribute to the lower energy barriers in KDM6B vs. KDM6A. The study reveals a dependence of the activation barrier of the rebound hydroxylation on the Fe-O-C angle in the transition state of KDM6A. Anti-correlation of the Zn-binding domain with the active site residues is a key factor distinguishing KDM6A/B from KDM7/4s. The results reveal the importance of communication between the Fe center, second coordination sphere, and long-range interactions in catalysis by KDMs and, by implication, other 2OG oxygenases.

QM/MM Modeling Aided Enzyme Engineering in Natural Products Biosynthesis

Natural products and their derivatives are widely used across various industries, particularly pharmaceuticals. Modern engineered biosynthesis provides an alternative way of producing and meeting the growing need for diverse natural products. Natural enzymes, on the other hand, often exhibit unsatisfactory catalytic characteristics and necessitate further enzyme engineering modifications. QM/MM, as a powerful and extensively used computational tool in the field of enzyme catalysis, has been increasingly applied in rational enzyme engineering over the past decade. In this review, we summarize recent advances in QM/MM computational investigation on enzyme catalysis and enzyme engineering for natural product biosynthesis. The challenges and perspectives for future QM/MM applications aided enzyme engineering in natural product biosynthesis will also be discussed.

Regioselective oxidation of heterocyclic aromatic hydrocarbons catalyzed by cytochrome P450: A case study of carbazole

Recently there are increasing interests in accurately evaluating the health effects of heterocyclic PAHs. However, the activation mechanism and possible metabolites of heterocyclic PAHs catalyzed by human CYP1A1 is still elusive to a great extent. Here, leveraged to high level QM/MM calculations, the corresponding activation pathways of a representative heterocyclic PAHs, carbazole, were systematically explored. The first stage is electrophilic addition or hydrogen abstraction from N-H group. Electrophilic addition was evidenced to be more feasible and regioselectivity at C3 and C4 sites were identified. Correlations between energy barriers and key structural/electrostatic parameters reveal that O-Cα distance and Fe-O-Cα angle are the main origin for the catalytic regioselectivity. Electrophilic addition was determined as the rate-determining step and the subsequent possible reactions include epoxidation, NIH shift (the hydrogen migration from the site of hydroxylation to the adjacent carbon) and proton shuttle. The corresponding products are epoxides, ketones and hydroxylated carbazoles, respectively. The main metabolites (hydroxylated carbazoles) are estimated to be more toxic than carbazole. The regioselectivity of carbazole activated by CYP1A1 is different from the environmental processes (gas and aqueous phase). Collectively, these results will inform the in-depth understanding the metabolic processes of heterocyclic PAHs and aid the accurate evaluation of their health effects.

Enantioselectivity and origin of enhanced efficiency in polyethylene terephthalate hydrolases catalyzed depolymerization

Biotechnology is one of the most promising strategies to resolve the global crisis of plastic pollution. A clear understanding of the core enzyme mechanisms in the biotransformation process is critical for rational enzyme engineering and for practical, industrial-scale applications. Herein, we systematically examined and evidenced a largely unexplored piece in the depolymerization mechanism catalyzed by polyethylene terephthalate (PET) hydrolases: their enantioselectivity. We found that all the short-lived tetrahedron intermediates (IM3 and IM8) possess S-type chirality in six representative PET hydrolases. For instance, the binding percentage ratio of pro-S:pro-R is 57:21 in FAST-PETase, while pro-S binding leads to a much lower average energy barrier (5.2 kcal/mol) than pro-R binding (33.1 kcal/mol). Key structural features (e.g. the angle for Ser@H1-His@N1-PET@O2 and distance for His@N1-PET@O2) that significantly modulate the enantioselectivity were identified. The origin of the energy landscape variation between wild-type IsPETase and mutant FAST-PETase was also unveiled via analysis of key features, the distortion/interaction energy, and non-covalent bond interactions. This study supplies the missing piece in the mechanism for depolymerization catalyzed by PET hydrolases, and will aid in the rational design of enzymes for industrial recycling of PET plastic waste.

Insight into the Role of Active Site Protonation States and Water Molecules in the Catalytic Inhibition of DPP4 by Vildagliptin

Vildagliptin (VIL) is an antidiabetic drug that inhibits dipeptidyl peptidase-4 (DPP4) through a covalent mechanism. The molecular bases for this inhibitory process have been addressed experimentally and computationally. Nevertheless, relevant issues remain unknown such as the roles of active site protonation states and conserved water molecules nearby the catalytic center. In this work, molecular dynamics simulations were applied to examine the structures of 12 noncovalent VIL-DPP4 complexes encompassing all possible protonation states of three noncatalytic residues (His126, Asp663, Asp709) that were inconclusively predicted by different computational tools. A catalytically competent complex structure was only achieved in the system with His126 in its ε-form and nonconventional neutral states for Asp663/Asp709. This complex suggested the involvement of one water molecule in the catalytic process of His740/Ser630 activation, which was confirmed by QM/MM simulations. Our findings support the suitability of a novel water-mediated mechanism in which His740/Ser630 activation occurs concertedly with the nucleophilic attack on VIL and the imidate protonation by Tyr547. Then, the restoration of His740/ Tyr547 protonation states occurs via a two-water hydrogen bonding network in a low-barrier process, thus describing the final step of the catalytic cycle for the first time. Additionally, two hydrolytic mechanisms were proposed based on the hydrogen bonding networks formed by water molecules and the catalytic residues along the inhibitory mechanism. These findings are valuable to unveil the molecular features of the covalent inhibition of DPP4 by VIL and support the future development of novel derivatives with improved structural or mechanistic profiles.

Human Cytochrome P450 1, 2, 3 Families as Pharmacogenes with Emphases on Their Antimalarial and Antituberculosis Drugs and Prevalent African Alleles

Precision medicine gives individuals tailored medical treatment, with the genotype determining the therapeutic strategy, the appropriate dosage, and the likelihood of benefit or toxicity. Cytochrome P450 (CYP) enzyme families 1, 2, and 3 play a pivotal role in eliminating most drugs. Factors that affect CYP function and expression have a major impact on treatment outcomes. Therefore, polymorphisms of these enzymes result in alleles with diverse enzymatic activity and drug metabolism phenotypes. Africa has the highest CYP genetic diversity and also the highest burden of malaria and tuberculosis, and this review presents current general information on CYP enzymes together with variation data concerning antimalarial and antituberculosis drugs, while focusing on the first three CYP families. Afrocentric alleles such as CYP2A6*17, CYP2A6*23, CYP2A6*25, CYP2A6*28, CYP2B6*6, CYP2B6*18, CYP2C8*2, CYP2C9*5, CYP2C9*8, CYP2C9*9, CYP2C19*9, CYP2C19*13, CYP2C19*15, CYP2D6*2, CYP2D6*17, CYP2D6*29, and CYP3A4*15 are implicated in diverse metabolic phenotypes of different antimalarials such as artesunate, mefloquine, quinine, primaquine, and chloroquine. Moreover, CYP3A4, CYP1A1, CYP2C8, CYP2C18, CYP2C19, CYP2J2, and CYP1B1 are implicated in the metabolism of some second-line antituberculosis drugs such as bedaquiline and linezolid. Drug-drug interactions, induction/inhibition, and enzyme polymorphisms that influence the metabolism of antituberculosis, antimalarial, and other drugs, are explored. Moreover, a mapping of Afrocentric missense mutations to CYP structures and a documentation of their known effects provided structural insights, as understanding the mechanism of action of these enzymes and how the different alleles influence enzyme function is invaluable to the advancement of precision medicine.

The Performance of Different Water Models on the Structure and Function of Cytochrome P450 Enzymes

Modeling approaches and modern simulations to investigate the biomolecular structure and function rely on various methods. Since water molecules play a crucial role in all sorts of chemistry, the accurate modeling of water molecules is vital for such simulations. In cytochrome P450 (CYP450), in particular, water molecules play a key role in forming active oxidant that ultimately performs oxidation and metabolism. In the present study, we have highlighted the behavior of the three most widely used water models─TIP3P, SPC/E, and OPC─for three different CYP450 enzymes─CYP450_BM3, CYP450_OleT, and CYP450_BSβ─during MD simulations and QM/MM calculations. We studied the various properties, such as RMSD, RMSF, H-bond, water occupancy, and hydrogen atom transfer (HAT), using QM/MM calculations and compared them for all three water models. Our study shows that the stabilities of the enzyme complexes are well maintained in all three water models. However, the OPC water model performs well for the polar active sites, that is, in CYP450_OleT and CYP450_BSβ, while the TIP3P water model is superior for the hydrophobic site, such as CYP450_BM3.

Mechanisms of Electron Transfer Rate Modulations in Cytochrome P450 BM3

Bacterial cytochromes P450 BM3 (CYP450 BM3) catalyze reactions of industrial importance. Despite many successful biotransformations, robust (re)design for novel applications remains challenging. Rational design and evolutionary approaches are not always successful, highlighting a lack of complete understanding of the mechanisms of electron transfer (ET) modulations. Thus, the full potential of CYP450 reactions remains under-exploited. In this work, we report the first molecular dynamics (MD)-based explicit prediction of BM3 ET parameters (reorganization energies; λ and ET free energies; ΔG°), and log ET rates (log k_ET) using the Marcus theory. Overall, the calculated ET rates for the BM3 wild-type (WT), mutants (F393 and L86), ligand-bound state, and ion concentrations agree well with experimental data. In ligand-free (LF) BM3, mutations modulate k_ET via ET ΔG°. Simulations show that the experimental ET rate enhancement is due to increased driving force (more negative ΔG°) upon ligation. This increase is related to the protein reorganization required to accommodate the ligand in the binding pocket rather than binding interactions with the ligand. Our methodology (CYPWare 1.0) automates all the stages of the MD simulation step-up, energy calculations, and estimation of ET parameters. CYPWare 1.0 and this work thus represent an important advancement in the CYP450 ET rate predictions, which has the potential to guide the redesign of ET enzymes. This program and a Web tool are available on GitHub for academic research.

Predicting Regioselectivity of AO, CYP, FMO, and UGT Metabolism Using Quantum Mechanical Simulations and Machine Learning

Unexpected metabolism in modification and conjugation phases can lead to the failure of many late-stage drug candidates or even withdrawal of approved drugs. Thus, it is critical to predict the sites of metabolism (SoM) for enzymes, which interact with drug-like molecules, in the early stages of the research. This study presents methods for predicting the isoform-specific metabolism for human AOs, FMOs, and UGTs and general CYP metabolism for preclinical species. The models use semi-empirical quantum mechanical simulations, validated using experimentally obtained data and DFT calculations, to estimate the reactivity of each SoM in the context of the whole molecule. Ligand-based models, trained and tested using high-quality regioselectivity data, combine the reactivity of the potential SoM with the orientation and steric effects of the binding pockets of the different enzyme isoforms. The resulting models achieve κ values of up to 0.94 and AUC of up to 0.92.

The involvement of CYP1A2 in biodegradation of dioxins in pigs

2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is one of the most harmful chemicals showing resistance to biodegradation. The majority of TCDD effects is mediated by the aryl hydrocarbon receptor (AhR) pathway. TCDD binding to AhR results in the activation of cytochrome P450 enzymes (CYP1A1, CYP1A2, CYP1B1) involved in dioxin biodegradation. The goal of the study was to explore the potentialrole of CYP1A2 in the metabolism of TCDD. We investigated a molecular structure of CYP1A2 and the binding selectivity and affinity between the pig CYP1A2 and: 1/ DiCDD or TCDD (dioxins differing in toxicity and biodegradability) or 2/ their selected metabolites. pCYP1A2 demonstrated higher affinity towards DiCDD and TCDD than other pCYP1 enzymes. All dioxin-pCYP1A2 complexes were found to be stabilized by hydrophobic interactions. The calculated distances between the heme oxygen and the dioxin carbon nearest to the oxygen, reflecting the hydroxylating potential of CYP1A2, were higher than in other pCYP1 enzymes. The distances between the heme iron and the nearest dioxin carbon exceeded 5 Å, a distance sufficient to allow the metabolites to leave the active site. However, the molecular dynamics simulations revealed that two access channels of CYP1A2 were closed upon binding the majority of the examined dioxins. Moreover, the binding of dioxin metabolites did not promote opening of channel S–an exit for hydroxylated products. It appears that the undesired changes in the behavior of access channels prevail over the hydroxylating potential of CYP1A2 towards TCDD and the favorable distances, ultimately trapping the metabolites at the enzyme’s active site.

Current state and future perspectives of cytochrome P450 enzymes for C–H and C=C oxygenation

Cytochrome P450 enzymes (CYPs) catalyze a series of C–H and C=C oxygenation reactions, including hydroxylation, epoxidation, and ketonization. They are attractive biocatalysts because of their ability to selectively introduce oxygen into inert molecules under mild conditions. This review provides a comprehensive overview of the C–H and C=C oxygenation reactions catalyzed by CYPs and the various strategies for achieving higher selectivity and enzymatic activity. Furthermore, we discuss the application of C–H and C=C oxygenation catalyzed by CYPs to obtain the desired chemicals or pharmaceutical intermediates in practical production. The rapid development of protein engineering for CYPs provides excellent biocatalysts for selective C–H and C=C oxygenation reactions, thereby promoting the development of environmentally friendly and sustainable production processes.

Computational Study of the Peroxygenase Mechanism Catalyzed by Hemoglobin Dehaloperoxidase Involved in the Degradation of Chlorophenols

The biochemical evidence showed that hemoglobin dehaloperoxidase (DHP B) from Amphitrite Ornata is a multifunctional hemoprotein that catalyzes both dehalogenation and hydroxylation of halophenols via the peroxidase and peroxygenase mechanism, respectively, which sets the basis for the degradation of halophenols. In the peroxygenase mechanism, the reaction was previously suggested to be triggered either by the hydrogen atom abstraction by the Fe═O center or by the proton abstraction by His55. To illuminate the peroxygenase mechanism of DHP B at the atomistic level, on the basis of the high-resolution crystal structure, computational models were constructed, and a series of quantum mechanical/molecular mechanical calculations have been performed. According to the calculation results, the pathway (Path a) initiated by the H-abstraction by the Fe═O center is feasible. In another pathway (Path b), His55 can abstract the proton from the hydroxyl group of the substrate (4-Cl-o-cresol) to initiate the reaction; however, its feasibility depends on the prior electron transfer from the substrate to the porphyrin group. The rate-limiting step of Path a is the OH-rebound, which corresponds to an energy barrier of 14.7 kcal/mol at the quartet state. His55 acts as an acid-base catalyst and directly involves in the catalysis. Our mutant study indicates that His55 can be replaced by other titratable residues. These findings may provide useful information for further understanding of the catalytic reaction of DHP B and for the design of enzymes in the degradation of pollutants, in particular, halophenols.

MD simulations and QM/MM calculations reveal the key mechanistic elements which are responsible for the efficient C-H amination reaction performed by a bioengineered P450 enzyme

An enzyme which is capable of catalyzing C–H amination reactions is considered to be a dream tool for chemists due to its pharmaceutical potential and greener approach. Recently, the Arnold group achieved this feat using an engineered CYP411 enzyme, which further undergoes a random directed evolution which increases its efficiency and selectivity. The present study provides mechanistic insight and the root cause of the success of these mutations to enhance the reactivity and selectivity of the mutant enzyme. This is achieved by means of comprehensive MD simulations and hybrid QM/MM calculations. The study shows that the efficient C–H amination by the engineered CYP411 is a combined outcome of electronic and steric effects. The mutation of the axial cysteine ligand to serine relays electron density to the Fe ion in the heme, and thereby enhances the bonding capability of the heme-iron to the nitrogen atom of the tosyl azide. In comparison, the native cysteine-ligated P450 cannot bind the tosyl azide. Additionally, the A78V and A82L mutations in P411 provide ‘bulk’ to the active site which increases the enantioselectivity via a steric effect. At the same time, the QM/MM calculations elucidate the C–H amination by the iron nitrenoid, revealing a mechanism analogous to Compound I in the native C–H hydroxylation by P450.

Computer simulation method reveals the mechanism of C–H amination reaction due to a single site mutation.

Theoretical Investigation on the Elusive Reaction Mechanism of Spirooxindole Formation Mediated by Cytochrome P450s: A Nascent Feasible Charge-Shift C-O Bond Makes a Difference

Spirooxindoles are pivotal biofunctional groups widely distributed in natural products and clinic drugs. However, construction of such subtle chiral skeletons is a long-standing challenge to both organic and bioengineering scientists. The knowledge of enzymatic spirooxindole formation in nature may inspire rational design of new catalysts. To this end, we presented a theoretical investigation on the elusive mechanism of the spiro-ring formation at the 3-position of oxindole mediated by cytochrome P450 enzymes (P450). Our calculated results demonstrated that the electrophilic attack of CpdI, the active species of P450, to the substrate, shows regioselectivity, i.e., the attack at the C9 position forms a tetrahedral intermediate involving an unusual feasible charge-shift C9^δ+-O^δ- bond, while the attack at the C1 position forms an epoxide intermediate. The predominant route is the first route with the charge-shift bonding intermediate due to holding a relatively lower barrier by >5 kcal mol^-1 than the epoxide route, which fits the experimental observations. Such a delocalized charge-shift bond facilitates the formation of a spiro-ring mainly through elongation of the C1-C9 bond to eliminate the aromatization of the tricyclic beta-carboline. Our theoretical results shed profound mechanistic insights for the first time into the elusive spirooxindole formation mediated by P450s.

Computer-aided understanding and engineering of enzymatic selectivity

Enzymes offering chemo-, regio-, and stereoselectivity enable the asymmetric synthesis of high-value chiral molecules. Unfortunately, the drawback that naturally occurring enzymes are often inefficient or have undesired selectivity toward non-native substrates hinders the broadening of biocatalytic applications. To match the demands of specific selectivity in asymmetric synthesis, biochemists have implemented various computer-aided strategies in understanding and engineering enzymatic selectivity, diversifying the available repository of artificial enzymes. Here, given that the entire asymmetric catalytic cycle, involving precise interactions within the active pocket and substrate transport in the enzyme channel, could affect the enzymatic efficiency and selectivity, we presented a comprehensive overview of the computer-aided workflow for enzymatic selectivity. This review includes a mechanistic understanding of enzymatic selectivity based on quantum mechanical calculations, rational design of enzymatic selectivity guided by enzyme-substrate interactions, and enzymatic selectivity regulation via enzyme channel engineering. Finally, we discussed the computational paradigm for designing enzyme selectivity in silico to facilitate the advancement of asymmetric biosynthesis.

Multiscale Simulations on the Catalytic Plasticity of CYP76AH1

The catalytic promiscuity and fidelity of cytochrome P450 enzymes are widespread in the skeletal modification of terpenoid natural products and have attracted much attention. CYP76AH1 is involved in key modification reactions in the biosynthetic pathway of tanshinone, a well-known medicinal norditerpenoid. In this work, classical molecular dynamic simulations, metadynamics, and DFT calculations were performed to investigate the protein conformational dynamics, ligand binding poses, and catalytic reaction mechanism in wide-type and mutant CYP76AH1. Our results not only reveal a plausible enzymatic mechanism for mutant CYP76AH1 leading to various products but also provide valuable guidance for rational protein engineering of the CYP76 family.

PCPD: Plant cytochrome P450 database and web-based tools for structural construction and ligand docking

Plant cytochrome P450s play key roles in the diversification and functional modification of plant natural products. Although over 200,000 plant P450 gene sequences have been recorded, only seven crystalized P450 genes severely hampered the functional characterization, gene mining and engineering of important P450s. Here, we combined Rosetta homologous modeling and MD-based refinement to construct a high-resolution P450 structure prediction process (PCPCM), which was applied to 181 plant P450s with identified functions. Furthermore, we constructed a ligand docking process (PCPLD) that can be applied for plant P450s virtual screening. 10 examples of virtual screening indicated the process can reduce about 80% screening space for next experimental verification. Finally, we constructed a plant P450 database (PCPD: http://p450.biodesign.ac.cn/), which includes the sequences, structures and functions of the 181 plant P450s, and a web service based on PCPCM and PCPLD. Our study not only developed methods for the P450-specific structure analysis, but also introduced a universal approach that can assist the mining and functional analysis of P450 enzymes.

Tryptophan-75 Is a Low-Energy Channel-Gating Residue that Facilitates Substrate Egress/Access in Cytochrome P450 2D6

CYP2D6 is a major drug metabolizing enzyme with a buried active site. Channels leading to the active site from various enzyme surfaces are believed to facilitate ligand egress and access to the active site. The present study used molecular dynamics (MD) and in vitro studies with CYP2D6*1 and a Trp75-to-Ala mutant to examine channel gating in CYP2D6 by Trp75. MD simulations measured energy landscapes of Trp75 conformations and simulated substrate passage within channel 2b using bufuralol as a model substrate. Trp75 alternated between multiple stable states that supported substrate transport along channel 2b with low-energy barriers between states (∼ -1 kcal/mol). Trp75 conformations were stabilized primarily by hydrogen bonding between Trp75 and Glu222, Asn226, Ala225, or Gln72. Energy barriers were low between Trp75 conformations, allowing Trp75 to easily move between various conformations over time and to function in both binding to and moving substrates in the 2b channel of CYP2D6. Michaelis-Menten kinetic studies completed with purified enzyme in a reconstituted system showed overall reduced enzyme efficiency for metabolism of bufuralol and dextromethorphan by the Trp75Ala mutant compared with CYP2D6*1. In stopped-flow measurements, k off for dextromethorphan was decreased in the absence of Trp75. Our results support a role for Trp75 in substrate shuttling to the active site of CYP2D6. SIGNIFICANCE STATEMENT: Using combined molecular dynamics and in vitro assays, this study shows for the first time a role for Trp75 as a channel entrance gating residue in the mechanism of substrate binding/unbinding in CYP2D6. Energy landscapes derived from molecular dynamics were used to quantitate the strength of gating, and kinetics assays showed the impact on enzyme efficiency and k off of a Trp75Ala mutation.

The molecular mechanism of P450-catalyzed amination of the pyrrolidine derivative of lidocaine: insights from multiscale simulations

Nitrogen heterocycles are key and prevalent motifs in drugs. Evolved variants of cytochrome P450_BM3 (CYP102A1) from Bacillus megaterium employ high-valent oxo-iron(iv) species to catalyze the synthesis of imidazolidine-4-ones via an intramolecular C–H amination. Herein, we use multi-scale simulations, including classical molecular dynamics (MD) simulations, quantum mechanical/molecular mechanical (QM/MM) calculations and QM calculations, to reveal the molecular mechanism of the intramolecular C–H amination of the pyrrolidine derivative of lidocaine bearing cyclic amino moieties catalyzed by the variant RP/FV/EV of P450_BM3, which bears five mutations compared to wild type. Our calculations show that overall catalysis includes both the enzymatic transformation in P450 and non-enzymatic transformation in water solution. The enzymatic transformation involves the exclusive hydroxylation of the C–H bond of the pyrrolidine derivative of lidocaine, leading to the hydroxylated intermediate, during which the substrate radical would be bypassed. The following dehydration and C–N coupling reactions are found to be much favored in aqueous situation compared to that in the non-polar protein environment. The present findings expand our understanding of the P450-catalyzed C(sp₃)–H amination reaction.

Nitrogen heterocycles are key and prevalent motifs in drugs.

Molecular Basis for Metabolic Regioselectivity and Mechanism of Cytochrome P450s toward Carcinogenic 4-(Methylnitrosamino)-(3-pyridyl)-1-butanone

As an abundantly present tobacco component, carcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK) has also been detected in atmospheric particulate matter, suggesting the ineluctable exposure risk of this contaminant. NNK metabolic activation by cytochrome P450 enzymes (CYPs) is a prerequisite to exerting its genotoxicity, but the metabolic regioselectivity and mechanism are still unknown. Here the binding feature and regioselectivity of CYPs 1A1, 1A2, 2A6, 2A13, 2B6, and 3A4 toward NNK are unraveled through molecular docking and molecular dynamics (MD) simulations. Binding mode analyses reveal that 1A2 and 2B6 have definite preferences for NNK α-methyl hydroxylation, while the other four CYPs preferentially catalyze α-methylene hydroxylation. The binding affinities between NNK and CYPs evaluated by the binding free energies follow the order 2A13 > 2B6 > 1A2 > 2A6 > 1A1 > 3A4. Density functional theory (DFT) calculations are further performed to characterize the mechanism of NNK biotransformation. Results show that the α-hydroxyNNK generated from α-hydroxylation may undergo nonenzymatic decomposition to form genotoxic diazohydroxide and aldehyde, and further oxidation by P450 to yield nitrosamide, which mainly contributes to NNK toxification capacity. Meanwhile the pyridine N-oxidation and denitrosation of C_α-radical intermediate play an important role in detoxifying NNK. Overall, the present study provides the molecular basis for CYP-catalyzed regioselectivity and mechanism of NNK biotransformation, which can enable the identification of metabolites for assessing the health risk of individual NNK exposure.

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)–H₂O₂ (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.

Regio- and stereoselectivity in the CYP450BM3-catalyzed hydroxylation of complex terpenoids: a QM/MM study

The selective oxidation of C–H in artemisinin by P450_BM3 variants was disclosed by combining QM/MM and MD simulations.

Insights into the biotransformation of 2,4,6-trinitrotoluene by the old yellow enzyme family of flavoproteins. A computational study

Biodegradation is a cost-effective and environmentally friendly alternative to removing 2,4,6-trinitrotoluene (TNT) pollution. However, mechanisms of TNT biodegradation have been elusive. To enhance the understanding of TNT biotransformation by the Old Yellow Enzyme (OYE) family, we investigated the crucial first-step hydrogen-transfer reaction by molecular dynamics simulations, docking technologies and empirical valence bond calculations. We revealed the significance of the π-π stacking conformation between the substrate TNT and the reduced flavin mononucleotide (FMNH2) cofactor, which is a prerequisite for the aromatic ring reduction of TNT. Under the π-π stacking conformation, the barrier of the hydrogen-transfer reaction in the aromatic ring reduction is about 16 kcal mol-1 lower than that of nitro group reduction. Then, we confirmed the mechanism of controlling the π-π stacking, that is, the π-π interaction competition mechanism. It indicates that the π-π stacking of TNT and FMNH2 occurs only when the π-π interaction between FMNH2 and TNT is stronger than that between TNT and several key residues with aromatic rings. Finally, based on the competition mechanism, the formation of π-π stacking of TNT and FMNH2 can be successfully enabled by removing the aromatic ring of those key residues in enzymes that originally only transform TNT through the nitro group reduction. This testified the validity of the π-π interaction competition mechanism. This work theoretically clarifies the molecular mechanism of the first-step hydrogen-transfer reaction for the biotransformation of TNT by the OYE family. It is helpful to obtain the enzymes that can biodegrade TNT through the aromatic ring reduction.

Use of Computational Biochemistry for Elucidating Molecular Mechanisms of Nitric Oxide Synthase

Nitric oxide (NO) is an essential signaling molecule in the regulation of multiple cellular processes. It is endogenously synthesized by NO synthase (NOS) as the product of L-arginine oxidation to L-citrulline, requiring NADPH, molecular oxygen, and a pterin cofactor. Two NOS isoforms are constitutively present in cells, nNOS and eNOS, and a third is inducible (iNOS). Despite their biological relevance, the details of their complex structural features and reactivity mechanisms are still unclear. In this review, we summarized the contribution of computational biochemistry to research on NOS molecular mechanisms. We described in detail its use in studying aspects of structure, dynamics and reactivity. We also focus on the numerous outstanding questions in the field that could benefit from more extensive computational investigations.

Cytochrome P450-The Wonderful Nanomachine Revealed through Dynamic Simulations of the Catalytic Cycle

This Account addresses the catalytic cycle of the enzyme cytochrome P450 (CYP450) as a prototypical biological machine with automatic features. CYP450 is a nanomachine that uses dioxygen and two reducing and two proton equivalents to oxidize a plethora of molecules (so-called substrates) as a means of supplying bio-organisms with essential molecules (e.g., brain neurotransmitters, sex hormones, etc.) and protecting biosystems against poisoning. An enticing property of CYP450s is that entrance of an oxidizable substrate into the active site initiates a series of events that constitute the catalytic cycle, which functions "automatically" in a regulated sequence of events culminating in the production of the oxidized substrates (e.g., hydroxylated, epoxidized, etc.), oftentimes with remarkable stereo- and regioselectivities. It is timely to demonstrate how theory uses molecular dynamics (MD) simulations and quantum-mechanical/molecular-mechanical (QM/MM) calculations to complement experiments and elucidate the choreography by which the protein regulates the catalytic cycle. CYP450 is a heme enzyme that contains a ferric ion (Fe^III) coordinated by a porphyrin ligand, a water molecule, and a cysteinate ligand that is provided by a strategic residue of the encapsulating protein. While many of the individual steps are sufficiently well-understood, we shall provide here an overview of the factors that cause all of the steps to be sequentially coordinated. To this end, we use examples from three different CYP450 enzymes: the bacterial ones CYP450_BM3 and CYP450_CAM and the mammalian enzyme CYP450_3A4. The treatment is limited to the catalytic cycle, as aspects of two-state reactivity were reviewed previously (e.g., Shaik , S. ; et al. Chem. Rev. 2005 , 105 , 2279 ). What are the principles that govern the seeming automatic feature? For example, how do substrate entrance and binding gate the enzyme? How does the reductase attachment to the enzyme affect the next steps? What triggers the attachment of the reductase? How does the electron transfer (ET) that converts Fe^III to Fe^II occur? Is the ET coordinated with the entrance of O₂ into the active site? What is the mechanism of the latter step? Since the entrance of the substrate expels the water molecules from the active site, how do water molecules re-enter to form a proton channel, which is necessary for creating the ultimate oxidant Compound I? How do mutations that disrupt the water channel nevertheless create a competent oxidant? By what means does the enzyme produce regio- and stereoselective oxidation products? What triggers the departure of the oxidized product, and how does the exit occur in a manner that generates the resting state ready for the next cycle? This Account shows that the entrance of the substrate triggers all of the ensuing events.

The Catalytic Mechanism of Steroidogenic Cytochromes P450 from All-Atom Simulations: Entwinement with Membrane Environment, Redox Partners, and Post-Transcriptional Regulation

Cytochromes P450 (CYP450s) promote the biosynthesis of steroid hormones with major impact on the onset of diseases such as breast and prostate cancers. By merging distinct functions into the same catalytic scaffold, steroidogenic CYP450s enhance complex chemical transformations with extreme efficiency and selectivity. Mammalian CYP450s and their redox partners are membrane-anchored proteins, dynamically associating to form functional machineries. Mounting evidence signifies that environmental factors are strictly intertwined with CYP450s catalysis. Atomic-level simulations have the potential to provide insights into the catalytic mechanism of steroidogenic CYP450s and on its regulation by environmental factors, furnishing information often inaccessible to experimental means. In this review, after an introduction of computational methods commonly employed to tackle these systems, we report the current knowledge on three steroidogenic CYP450s—CYP11A1, CYP17A1, and CYP19A1—endowed with multiple catalytic functions and critically involved in cancer onset. In particular, besides discussing their catalytic mechanisms, we highlight how the membrane environment contributes to (i) regulate ligand channeling through these enzymes, (ii) modulate their interactions with specific protein partners, (iii) mediate post-transcriptional regulation induced by phosphorylation. The results presented set the basis for developing novel therapeutic strategies aimed at fighting diseases originating from steroid metabolism dysfunction.

Open-Source, Python-Based Redevelopment of the ChemShell Multiscale QM/MM Environment

ChemShell is a scriptable computational chemistry environment with an emphasis on multiscale simulation of complex systems using combined quantum mechanical and molecular mechanical (QM/MM) methods. Motivated by a scientific need to efficiently and accurately model chemical reactions on surfaces and within microporous solids on massively parallel computing systems, we present a major redevelopment of the ChemShell code, which provides a modern platform for advanced QM/MM embedding models. The new version of ChemShell has been re-engineered from the ground up with a new QM/MM driver module, an improved parallelization framework, new interfaces to high performance QM and MM programs, and a user interface written in the Python programming language. The redeveloped package is capable of performing QM/MM calculations on systems of significantly increased size, which we illustrate with benchmarks on zirconium dioxide nanoparticles of over 160000 atoms.

Computational Understanding of the Selectivities in Metalloenzymes

Metalloenzymes catalyze many different types of biological reactions with high efficiency and remarkable selectivity. The quantum chemical cluster approach and the combined quantum mechanics/molecular mechanics methods have proven very successful in the elucidation of the reaction mechanism and rationalization of selectivities in enzymes. In this review, recent progress in the computational understanding of various selectivities including chemoselectivity, regioselectivity, and stereoselectivity, in metalloenzymes, is discussed.

Quantum chemical simulations revealed the toxicokinetic mechanisms of organic phosphorus flame retardants catalyzed by P450 enzymes

The metabolic fate and toxicokinetics of organic phosphorus flame retardants catalyzed by cytochrome P450 enzymes (CYPs) are here investigated by in silico simulations, leveraging an active center model to mimic the CYPs, triphenyl phosphate (TPHP), tris(2-butoxyethyl) phosphate and tris(1,3-dichloro-2-propyl) phosphate as substrates. Our calculations elucidated key main pathways and predicted products, which were corroborated by current in vitro data. Results showed that alkyl OPFRs are eliminated faster than aryl and halogenated alkyl-substituted OPFRs. In addition, we discovered a proton shuttle pathway for aryl hydroxylation of TPHP and P = O bond-assisted H-transfer mechanisms (rather than nonenzymatic hydrolysis) that lead to O-dealkylation/dearylation of phosphotriesters.

Is CYP1B1 involved in the metabolism of dioxins in the pig?

2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is the most difficult to biodegradate and the most toxic dioxin congener. Previously, we demonstrated in silico the ability of pig CYP1A1 to hydroxylate 2,7-dichlorodibenzo-p-dioxin (DiCDD), but not TCDD. To increase our knowledge concerning the low effectiveness of TCDD biodegradability, we analyzed in silico the binding selectivity and affinity between pig CYP1B1 and the two dioxins by means of molecular modeling. We also compared the effects of TCDD and DiCDD on CYP1B1 gene expression (qRT-PCR) and catalytic (EROD) activity in porcine granulosa cells. It was found that DiCDD and TCDD were stabilized within the pig CYP1B1 active site by hydrophobic interactions. The analysis of substrate channel availability revealed that both dioxins opened the exit channel S, allowing metabolites to leave the enzyme active site. Moreover, DiCDD and TCDD increased the CYP1B1 gene expression and catalytic activity in porcine granulosa cells. On the other hand, TCDD demonstrated higher than DiCDD calculated affinity to pig CYP1B1, hindering TCDD exit from the active site. The great distance between CYP1B1's heme and TCDD also might contribute to the lower hydroxylation effectiveness of TCDD compared to that of DiCDD. Moreover, the narrow active site of pig CYP1B1 may immobilize TCDD molecule, inhibiting its hydroxylation. The results of the access channel analysis and the distance from pig CYP1B1's heme to TCDD suggest that the metabolizing potential of pig CYP1B1 is higher than that of pig CYP1A1. However, this potential is probably not sufficiently high to considerably improve the slow TCDD biodegradation.

Computational Insight Into Vitamin K1 ω-Hydroxylation by Cytochrome P450 4F2

Vitamin K₁ (VK1) plays an important role in the modulation of bleeding disorders. It has been reported that ω-hydroxylation on the VK1 aliphatic chain is catalyzed by cytochrome P450 4F2 (CYP4F2), an enzyme responsible for the metabolism of eicosanoids. However, the mechanism of VK1 ω-hydroxylation by CYP4F2 has not been disclosed. In this study, we employed a combination of quantum mechanism (QM) calculations, homology modeling, molecular docking, molecular dynamics (MD) simulations, and combined quantum mechanism/molecular mechanism (QM/MM) calculations to investigate the metabolism profile of VK1 ω-hydroxylation. QM calculations based on the truncated VK1 model show that the energy barrier for ω-hydroxylation is about 6-25 kJ/mol higher than those at other potential sites of metabolism. However, results from the MD simulations indicate that hydroxylation at the ω-site is more favorable than at the other potential sites, which is in accordance with the experimental observation. The evaluation of MD simulations was further endorsed by the QM/MM calculation results. Our studies thus suggest that the active site residues of CYP4F2 play a determinant role in the ω-hydroxylation. Our results provide structural insights into the mechanism of VK1 ω-hydroxylation by CYP4F2 at the atomistic level and are helpful not only for characterizing the CYP4F2 functions but also for looking into the ω-hydroxylation mediated by other CYP4 enzymes.

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.