A Paradigm Shift in the Catalytic Cycle of P450: The Preparatory Choreography during O2 Binding and Origins of the Necessity for Two Protonation Pathways

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

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

Two Cytochrome P450 Enzymes Form the Tricyclic Nested Skeleton of Meroterpenoids by Sequential Oxidative Reactions

Meroterpenoid clavilactones feature a unique benzo-fused ten-membered carbocyclic ring unit with an α,β-epoxy-γ-lactone moiety, forming an intriguing 10/5/3 tricyclic nested skeleton. These compounds are good inhibitors of the tyrosine kinase, attracting a lot of chemical synthesis studies. However, the natural enzymes involved in the formation of the 10/5/3 tricyclic nested skeleton remain unexplored. Here, we identified a gene cluster responsible for the biosynthesis of clavilactone A in the basidiomycetous fungus Clitocybe clavipes. We showed that a key cytochrome P450 monooxygenase ClaR catalyzes the diradical coupling reaction between the intramolecular hydroquinone and allyl moieties to form the benzo-fused ten-membered carbocyclic ring unit, followed by the P450 ClaT that exquisitely and stereoselectively assembles the α,β-epoxy-γ-lactone moiety in clavilactone biosynthesis. ClaR unprecedentedly acts as a macrocyclase to catalyze the oxidative cyclization of the isopentenyl to the nonterpenoid moieties to form the benzo-fused macrocycle, and a multifunctional P450 ClaT catalyzes a ten-electron oxidation to accomplish the biosynthesis of the 10/5/3 tricyclic nested skeleton in clavilactones. Our findings establish the foundation for the efficient production of clavilactones using synthetic biology approaches and provide the mechanistic insights into the macrocycle formation in the biosynthesis of fungal meroterpenoids.

On the engineering of reductase-based-monooxygenase activity in CYP450 peroxygenases

Recent bioengineering of CYP450OleT shows that peroxide-based CYP450OleT can be converted to a reductase-based self-sufficient enzyme, which is capable of showing efficient hydroxylation and decarboxylation activity for a wide range of substrates. The so-generated enzyme creates several mechanistic puzzles: (A) as CYP450 peroxygenases lack the conventional acid-alcohol pair, what is the source of two protons that are required to create the ultimate oxidant Cpd I? (B) Why is it only CYP450OleT that shows the reductase-based activity but no other CYP members? The present study provides a mechanistic solution to these puzzles using comprehensive MD simulations and hybrid QM/MM calculations. We show that the fusion of the reductase domain to the heme-binding domain triggers significant conformational rearrangement, which is gated by the propionate side chain, which constitutes a new water aqueduct via the carboxylate end of the substrate that ultimately participates in Cpd I formation. Importantly, such well-synchronized choreographies are controlled by remotely located Tyr359, which senses the fusion of reductase and communicates to the heme domain via non-covalent interactions. These findings provide crucial insights and a broader perspective which enables us to make a verifiable prediction: thus, the catalytic activity is not only limited to the first or second catalytic shell of an enzyme. Furthermore, it is predicted that reinstatement of tyrosine at a similar position in other members of CYP450 peroxygenases can convert these enzymes to reductase-based monooxygenases.

Multiscale QM/MM modelling of catalytic systems with ChemShell

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.

Mechanistic Insights into Photochemical CO2 Reduction to CH4 by a Molecular Iron-Porphyrin Catalyst

Iron tetraphenylporphyrin complex modified with four trimethylammonium groups (Fe-p-TMA) is found to be capable of catalyzing the eight-electron eight-proton reduction of CO₂ to CH₄ photochemically in acetonitrile. In the present work, density functional theory (DFT) calculations have been performed to investigate the reaction mechanism and to rationalize the product selectivity. Our results revealed that the initial catalyst Fe-p-TMA ([Cl-Fe(III)-LR₄]⁴⁺, where L = tetraphenylporphyrin ligand with a total charge of -2, and R₄ = four trimethylammonium groups with a total charge of +4) undergoes three reduction steps, accompanied by the dissociation of the chloride ion to form [Fe(II)-L^••2-R₄]²⁺. [Fe(II)-L^••2-R₄]²⁺, bearing a Fe(II) center ferromagnetically coupled with a tetraphenylporphyrin diradical, performs a nucleophilic attack on CO₂ to produce the ¹η-CO₂ adduct [CO₂^•--Fe(II)-L^•-R₄]²⁺. Two intermolecular proton transfer steps then take place at the CO₂ moiety of [CO₂^•--Fe(II)-L^•-R₄]²⁺, resulting in the cleavage of the C-O bond and the formation of the critical intermediate [Fe(II)-CO]⁴⁺ after releasing a water molecule. Subsequently, [Fe(II)-CO]⁴⁺ accepts three electrons and one proton to generate [CHO-Fe(II)-L^•-R₄]²⁺, which finally undergoes a successive four-electron-five-proton reduction to produce methane without forming formaldehyde, methanol, or formate. Notably, the redox non-innocent tetraphenylporphyrin ligand was found to play an important role in CO₂ reduction since it could accept and transfer electron(s) during catalysis, thus keeping the ferrous ion at a relatively high oxidation state. Hydrogen evolution reaction via the formation of Fe-hydride ([Fe(II)-H]³⁺) turns out to endure a higher total barrier than the CO₂ reduction reaction, therefore providing a reasonable explanation for the origin of the product selectivity.

How does multiple substrate binding lead to substrate inhibition of CYP2D6 metabolizing dextromethorphan? A theoretical study

CYP2D6 is one of the most important metalloenzymes involved in the biodegradation of many drug molecules in the human body. It has been found that multiple substrate binding can lead to substrate inhibition of CYP2D6 metabolizing dextromethorphan (DM), but the corresponding theoretical mechanism is rarely reported. Therefore, we chose DM as the probe and performed molecular dynamics simulations and quantum mechanical calculations on CYP2D6-DM systems to investigate the mechanism of how the multiple substrate binding leads to the substrate inhibition of CYP2D6 metabolizing substrates. According to our results, three gate residues (Arg221, Val374, and Phe483) for the catalytic pocket are determined. We also found that the multiple substrate binding can lead to substrate inhibition by reducing the stability of CYP2D6 binding DM and increasing the reactive activation energy of the rate-determining step. Our findings would help to understand the substrate inhibition of CYP2D6 metabolizing the DM and enrich the knowledge of the drug-drug interactions for the cytochrome P450 superfamily.

Mechanistic Insight into Peptidyl-Cysteine Oxidation by the Copper-Dependent Formylglycine-Generating Enzyme

The copper-dependent formylglycine-generating enzyme (FGE) catalyzes the oxygen-dependent oxidation of specific peptidyl-cysteine residues to formylglycine. Our QM/MM calculations provide a very likely mechanism for this transformation. The reaction starts with dioxygen binding to the tris-thiolate Cu^I center to form a triplet Cu^II -superoxide complex. The rate-determining hydrogen atom abstraction involves a triplet-singlet crossing to form a Cu^II -OOH species that couples with the substrate radical, leading to a Cu^I -alkylperoxo intermediate. This is accompanied by proton transfer from the hydroperoxide to the S atom of the substrate via a nearby water molecule. The subsequent O-O bond cleavage is coupled with the C-S bond breaking that generates the formylglycine and a Cu^II -oxyl complex. Moreover, our results suggest that the aldehyde oxygen of the final product originates from O₂ , which will be useful for future experimental work.

Efficient Production of Epoxy-Norbornane from Norbornene by an Engineered P450 Peroxygenase

Epoxy-norbornane (EPO-NBE) is a crucial building block for the synthesis of various biologically active heterocyclic systems. To develop an efficient protocol for producing EPO-NBE using norbornene (NBE) as a substrate, cytochrome P450 enzyme from Pseudomonas putida (CYP238A1) was examined and its crystal structure (PDB code: 7X53) was resolved. Molecular mechanism analysis showed a high energy barrier related to iron-alkoxy radical complex formation. Therefore, a protein engineering strategy was developed and an optimal CYP238A1^NPV variant containing a local hydrophobic "fence" at the active site was obtained, which increased the H₂ O₂ -dependent epoxidation activity by 7.5-fold compared with that of CYP238A1^WT . Among the "fence", Glu255 participates in an efficient proton transfer system. Whole-cell transformation using CYP238A1^NPV achieved an EPO-NBE yield of 77.6 g ⋅ L^-1 in a 30-L reactor with 66.3 % conversion. These results demonstrate the potential of this system for industrial production of EPO-NBE and provides a new biocatalytic platform for epoxidation chemistry.

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.

Computations reveal a crucial role of an aromatic dyad in the catalytic function of plant cytochrome P450 mint superfamily

Enzymes are highly specific for their native functions, however with advances in bioengineering tools such as directed evolution, several enzymes are being repurposed for the secondary function of contemporary significance(Khersonsky and Tawfik, 2010 [1]⁾. Due to the functional versatility, the Cytochrome P450 (CYP450) superfamily has become the ideal scaffold for such bioengineering. In the current study, using MD (molecular dynamics) simulations and hybrid QM/MM (Quantum mechanics/molecular mechanics) calculations, we have studied the mechanism of spontaneous emergence of a secondary function due to a single site mutation in two plant CYP450 enzymes from the mint family. The MD simulations of WT (wild type) CYP71D18 and CYP71D13 enzymes and their variants show a crucial gating mechanism by aromatic dyad formed by Phe121 and Phe363 which regulates the substrate recognition. The QM/MM calculations reveal that the hydroxylation reactions at C3 and C6 positions in WT CYP71D18 and CYP71D13 enzymes as well as their variants follow a hydrogen atom transfer (HAT) followed by a single electron transfer (SET) mechanism, which is different from the typical rebound mechanism shown by most of the CYP450 enzymes.

Single-Site Mutation Induces Water-Mediated Promiscuity in Lignin Breaking Cytochrome P450GcoA

Cytochrome P450_GcoA is an enzyme that catalyzes the guaiacol unit of lignin during the lignin breakdown via an aryl-O-demethylation reaction. This reaction is intriguing and is of commercial importance for its potential applications in the production of biofuel and plastic from biomass feedstock. Recently, the F169A mutation in P450_GcoA elicits a promiscuous activity for syringol while maintaining the native activity for guaiacol. Using comprehensive MD simulations and hybrid QM/MM calculations, we address, herein, the origin of promiscuity in P450_GcoA and its relevance to the specific activity toward lignin-derived substrates. Our study shows a crucial role of an aromatic dyad of F169 and F395 by regulating the water access to the catalytic center. The F169A mutation opens a water aqueduct and hence increases the native activity for G-lignin. We show that syringol binds very tightly to the WT enzyme, which blocks the conformational rearrangement needed for the second step of O-demethylation. The F169A creates an extra room favoring the conformational rearrangement in the 3-methoxycatechol (3MC) and second dose of the dioxygen insertion. Therefore, using MD simulations and complemented by thorough QM/MM calculations, our study shows how a single-site mutation rearchitects active site engineering for promiscuous syringol activity.

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.

Local Electric Fields Dictate Function: The Different Product Selectivities Observed for Fatty Acid Oxidation by Two Deceptively Very Similar P450-Peroxygenases OleT and BSβ

Cytochrome P450 peroxygenases use hydrogen peroxide to hydroxylate long-chain fatty acids by bypassing the use of O₂ and a redox partner. Among the peroxygenases, P450_OleT uniquely performs decarboxylation of fatty acids and production of terminal olefins. This route taken by P450_OleT is intriguing, and its importance is augmented by the practical importance of olefin production. As such, this mechanistic choice merits elucidation. To address this puzzle, we use hybrid QM/MM calculations and MD simulations for the OleT enzyme as well as for the structurally analogous enzyme, P450_BSβ. The study of P450_OleT reveals that the protonated His85 in the wild-type P450_OleT plays a crucial role in steering decarboxylation activity by stabilizing the corresponding hydroxoiron(IV) intermediate (Cpd II). In contrast, for P450_BSβ in which Q85 replaces H85, the respective Cpd II species is unstable and it reacts readily with the substrate radical by rebound, producing hydroxylation products. As shown, this single-site difference creates in P450_OleT a local electric field (LEF), which is significantly higher than that in P450_BSβ. In turn, these LEF differences are responsible for the different stabilities of the respective Cpd II/radical intermediates and hence for different functions of the two enzymes. P450_BSβ uses the common rebound mechanism and leads to hydroxylation, whereas P450_OleT proceeds via decarboxylation and generates terminal olefins. Olefin production projects the power of a single residue to alter the LEF and the enzyme's function.

Mechanism of Sugar Ring Contraction and Closure Catalyzed by UDP-d-apiose/UDP-d-xylose Synthase (UAXS)

Uridine diphosphate (UDP)-apiose/UDP-xylose synthase (UAXS) is a member of the short-chain dehydrogenase/reductase superfamily (SDR), which catalyzes the ring contraction and closure of UDP-d-glucuronic acid (UDP-GlcA), affording UDP-apiose and UDP-xylose. UAXS is a special enzyme that integrates ring-opening, decarboxylation, rearrangement, and ring closure/contraction in a single active site. Recently, the ternary complex structure of UAXS was crystallized from Arabidopsis thaliana. In this work, to gain insights into the detailed formation mechanism of UDP-apiose and UDP-xylose, an enzyme-substrate reactant model has been constructed and quantum mechanical/molecular mechanical (QM/MM) calculations have been performed. Our calculation results reveal that the reaction starts from the C4-OH oxidation, which is accompanied by the conformational transformation of the sugar ring from chair type to boat type. The sugar ring-opening is prior to decarboxylation, and the deprotonation of the C2-OH group is the prerequisite for sugar ring-opening. Moreover, the keto-enol tautomerization of the decarboxylated intermediate is a necessary step for ring closure/contraction. Based on our calculation results, more UDP-apiose product was expected, which is in line with the experimental observation. Three titratable residues, Tyr185, Cys100, and Cys140, steer the reaction by proton transfer from or to UDP-GlcA, and Arg182, Glu141, and D337 constitute a proton conduit for sugar C2-OH deprotonation. Although Thr139 and Tyr105 are not directly involved in the enzymatic reaction, they are responsible for promoting the catalysis by forming hydrogen-bonding interactions with GlcA. Our calculations may provide useful information for understanding the catalysis of the SDR family.

Pinene-Based Oxidative Synthetic Toolbox for Scalable Polyester Synthesis

Generation of renewable polymers is a long-standing goal toward reaching a more sustainable society, but building blocks in biomass can be incompatible with desired polymerization type, hampering the full implementation potential of biomaterials. Herein, we show how conceptually simple oxidative transformations can be used to unlock the inherent reactivity of terpene synthons in generating polyesters by two different mechanisms starting from the same α-pinene substrate. In the first pathway, α-pinene was oxidized into the bicyclic verbanone-based lactone and subsequently polymerized into star-shaped polymers via ring-opening polymerization, resulting in a biobased semicrystalline polyester with tunable glass transition and melting temperatures. In a second pathway, polyesters were synthesized via polycondensation, utilizing the diol 1-(1'-hydroxyethyl)-3-(2'-hydroxy-ethyl)-2,2-dimethylcyclobutane (HHDC) synthesized by oxidative cleavage of the double bond of α-pinene, together with unsaturated biobased diesters such as dimethyl maleate (DMM) and dimethyl itaconate (DMI). The resulting families of terpene-based polyesters were thereafter successfully cross-linked by either transetherification, utilizing the terminal hydroxyl groups of the synthesized verbanone-based materials, or by UV irradiation, utilizing the unsaturation provided by the DMM or DMI moieties within the HHDC-based copolymers. This work highlights the potential to apply an oxidative toolbox to valorize inert terpene metabolites enabling generation of biosourced polyesters and coatings thereof by complementary mechanisms.

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.

Dual-function enzyme catalysis for enantioselective carbon-nitrogen bond formation

Chiral amines can be made by insertion of a carbene into an N-H bond using two-catalyst systems that combine a transition metal-based carbene-transfer catalyst and a chiral proton-transfer catalyst to enforce stereocontrol. Haem proteins can effect carbene N-H insertion, but asymmetric protonation in an active site replete with proton sources is challenging. Here we describe engineered cytochrome P450 enzymes that catalyse carbene N-H insertion to prepare biologically relevant α-amino lactones with high activity and enantioselectivity (up to 32,100 total turnovers, >99% yield and 98% e.e.). These enzymes serve as dual-function catalysts, inducing carbene transfer and promoting the subsequent proton transfer with excellent stereoselectivity in a single active site. Computational studies uncover the detailed mechanism of this new-to-nature enzymatic reaction and explain how active-site residues accelerate this transformation and provide stereocontrol.

Developing C2-Aroyl Indoles as Novel Inhibitors of IDO1 and Understanding Their Mechanism of Inhibition via Mass Spectroscopy, QM/MM Calculations and Molecular Dynamics Simulation

Indoleamine-2,3-dioxygenase (IDO1) and tryptophan dioxygenases are two heme based metalloenzymes that catalyze the tryptophan oxidation reaction by inserting molecular dioxygen to cleave the pyrrole ring. The mechanism of such ring cleavage reaction is of carcinogenic importance as the malignant tumors recruit this mechanism for immune invasion. In the presence study, we have synthesized a Novel C2 aroyl indoles inhibitor, 8d, which shows significant inhibition of 180 nM at IC₅₀ scale. The binding and conformational changes that transpire after inhibitor binding were thoroughly studied by molecular docking and MD simulations. The subsequent QM/MM (Quantum Mechanical/Molecular Mechanical) calculations were used to proposed the mechanism of inhibition. The QM/MM calculations show that the reaction proceeds via multistep processes where the dioxygen insertion to the substrate 8a is the rate determining process. Theoretical mechanism is further supported by mass spectroscopy, and drug metabolism/pharmacokinetics study (DMPK) and metabolic stability of compound 8d was investigated in rat and human liver microsomes.