Machine learning-aided engineering of a cytochrome P450 for optimal bioconversion of lignin fragments.

Using machine learning, molecular dynamics simulations, and density functional theory calculations we gain insight into the selectivity patterns of substrate activation by the cytochromes P450. In nature, the reactions catalyzed by the P450s lead to the biodegradation of xenobiotics, but recent work has shown that fungi utilize P450s for the activation of lignin fragments, such as monomer and dimer units. These fragments often are the building blocks of valuable materials, including drug molecules and fragrances, hence a highly selective biocatalyst that can produce these compounds in good yield with high selectivity would be an important step in biotechnology. In this work a detailed computational study is reported on two reaction channels of two P450 isozymes, namely the O-deethylation of guaethol by CYP255A and the O-demethylation versus aromatic hydroxylation of p-anisic acid by CYP199A4. The studies show that the second-coordination sphere plays a major role in substrate binding and positioning, heme access, and in the selectivity patterns. Moreover, the local environment affects the kinetics of the reaction through lowering or raising barrier heights. Furthermore, we predict a site-selective mutation for highly specific reaction channels for CYP199A4.

Engineering Electron Transfer Pathway of Cytochrome P450s.

Cytochrome P450s (P450s), a superfamily of heme-containing enzymes, existed in animals, plants, and microorganisms. P450s can catalyze various regional and stereoselective oxidation reactions, which are widely used in natural product biosynthesis, drug metabolism, and biotechnology. In a typical catalytic cycle, P450s use redox proteins or domains to mediate electron transfer from NAD(P)H to heme iron. Therefore, the main factors determining the catalytic efficiency of P450s include not only the P450s themselves but also their redox-partners and electron transfer pathways. In this review, the electron transfer pathway engineering strategies of the P450s catalytic system are reviewed from four aspects: cofactor regeneration, selection of redox-partners, P450s and redox-partner engineering, and electrochemically or photochemically driven electron transfer.

Two Cytochrome P450s, CYP709B1 and CYP704C1, Play Essential Roles in Metabolism-Based Multiple Herbicide Resistance in American Sloughgrass (Beckmannia syzigachne (Steud.) Fernald).

This study confirmed a field population of American sloughgrass (Beckmannia syzigachne (Steud.) Fernald) that developed simultaneously high levels of resistance (resistance index >10) to three divergent modes of action herbicides: fenoxaprop-P-ethyl, mesosulfuron-methyl, and isoproturon. The resistance phenotype observed in this population was not attributed to target-site alterations; rather, the resistant plants exhibited a significant increase in the activity of cytochrome P450s (P450s) and enhanced metabolism rates for all three herbicides. RNA sequencing revealed significant upregulation of two P450s, CYP709B1 and CYP704C1, in the resistant plants both before and after herbicide treatments. Molecular docking predicted that the homology models of these P450s should exhibit a binding affinity for a range of herbicides. The heterologous expression of the identified P450s in yeast cells indicated improved growth in the presence of all three of the aforementioned herbicides. Collectively, the increased expression of CYP709B1 and CYP704C1 likely contributed to the P450s-mediated enhanced metabolism, thereby conferring multiple herbicide resistance in B. syzigachne.

Structure-Function Analysis of the Essential Mycobacterium tuberculosis P450 Drug Target, CYP121A1.

Cytochrome P450 CYP121A1 is a well-known drug target against Mycobacterium tuberculosis, the human pathogen that causes the deadly disease tuberculosis (TB). CYP121A1 is a unique P450 enzyme because it uses classical and non-classical P450 catalytic processes and has distinct structural features among P450s. However, a detailed investigation of CYP121A1 protein structures in terms of active site cavity dynamics and key amino acids interacting with bound ligands has yet to be undertaken. To address this research knowledge gap, 53 CYP121A1 crystal structures were investigated in this study. Critical amino acids required for CYP121A1's overall activity were identified and highlighted this enzyme's rigid architecture and substrate selectivity. The CYP121A1-fluconazole crystal structure revealed a novel azole drug-P450 binding mode in which azole heme coordination was facilitated by a water molecule. Fragment-based inhibitor approaches revealed that CYP121A1 can be inhibited by molecules that block the substrate channel or by directly interacting with the P450 heme. This study serves as a reference for the precise understanding of CYP121A1 interactions with different ligands and the structure-function analysis of P450 enzymes in general. Our findings provide critical information for the synthesis of more specific CYP121A1 inhibitors and their development as novel anti-TB drugs.

Mechanistic and predictive studies on the oxidation of furans by cytochrome P450: A DFT study.

Furan-containing compounds distribute widely in food, herbal medicines, industrial synthetic products, and environmental media. These compounds can undergo oxidative metabolism catalyzed by cytochrome P450 enzymes (CYP450) within organisms, which may produce reactive products, possibly reacting with biomolecules to induce toxic effects. In this work, we performed DFT calculations to investigate the CYP450-mediated metabolic mechanism of furan-ring oxidation using 2-methylfuran as a model substrate, meanwhile, we studied the regioselective competition of another hydroxylation reaction involving methyl group of 2-methylfuran. As a result, we found the toxicological-relevant cis-enedione product can be produced from O-addition directly via a concerted manner without formation of an epoxide intermediate as traditionally believed. Moreover, our calculations demonstrate the kinetic and thermodynamic feasibility of both furan-ring oxidation and methyl hydroxylation pathways, although the former pathway is a bit more favorable. We then constructed a linear model to predict the rate-limiting activation energies (ΔE*) of O-addition with 11 diverse furan substates based on their adiabatic ionization potentials (AIPs) and condensation Fukui functions (CFFs). The results show a good predictive ability (R2=0.94, Q2CV=0.87). Therefore, AIP and CFF with clear physichem meanings relevant to the mechanism, emerge as pivotal molecular descriptors to enable the fast prediction of furan-ring oxidation reactivities for quick insight into the toxicological risk of furans, using just ground-state calculations.

Rational Design of Triazinone Derivatives with Low Bee Toxicity Based on the Binding Mechanism of Neonicotinoids to Apis mellifera.

Bees, one of the most vital pollinators in the ecosystem and agriculture, are currently threatened by neonicotinoids. To explore the molecular mechanisms of neonicotinoid toxicity to bees, the different binding modes of imidacloprid, thiacloprid, and flupyradifurone with nicotinic acetylcholine receptor (nAChR) α1ß1 and cytochrome P450 9Q3 (CYP9Q3) were studied using homology modeling and molecular dynamics simulations. These mechanisms provided a basis for the design of compounds with a potential low bee toxicity. Consequently, we designed and synthesized a series of triazinone derivatives and assessed their bioassays. Among them, compound 5a not only displayed substantially insecticidal activities against Aphis glycines (LC50 = 4.40 mg/L) and Myzus persicae (LC50 = 6.44 mg/L) but also had low toxicity to Apis mellifera. Two-electrode voltage clamp recordings further confirmed that compound 5a interacted with the M. persicae nAChR α1 subunit but not with the A. mellifera nAChR α1 subunit. This work provides a paradigm for applying molecular toxic mechanisms to the design of compounds with low bee toxicity, thereby aiding the future rational design of eco-friendly nicotinic insecticides.

Bioinspired axial S-coordinated single-atom cobalt catalyst to efficient activate peroxymonosulfate for selective high-valent Co-Oxo species generation.

The efficient activation and selective high-valent metal-oxo (HVMO) species generation remain challenging for peroxymonosulfate (PMS)-based advanced oxidation processes (PMS-AOPs) in water purification. The underlying mechanism of the activation pathway is ambiguous, leading to a massive dilemma in the control and regulation of HVMO species generation. Herein, bioinspired by the bio-oxidase structure of cytochrome P450, the axial coordination strategy was adopted to tailor a single-atom cobalt catalyst (CoN4S-CB) with an axial S coordination. CoN4S-CB high-selectively generated high-valent Co-Oxo species (Co(IV)=O) via PMS activation. Co(IV)=O demonstrated an ingenious oxygen atom transfer (OAT) reaction to achieve the efficient degradation of sulfamethoxazole (SMX), and this allowed robust operation in various complex environments. The axial S coordination modulated the 3d orbital electron distribution of the Co atom. Density functional theory (DFT) calculation revealed that the axial S coordination decreased the energy barrier for PMS desorption and lowered the free energy change (ΔG) for Co(IV)=O generation. CoN4S-PMS* had a narrow d-band close to the Fermi level, which enhanced charge transfer to accelerate the cleavage of O-O and O-H bonds in PMS. This work provides a broader perspective on the activator design with natural enzyme structure-like active sites to efficient activate PMS for selective HVMO species generation.

Molecular Basis for Chemoselectivity Control in Oxidations of Internal Aryl-Alkenes Catalyzed by Laboratory Evolved P450s.

P450 enzymes naturally perform selective hydroxylations and epoxidations of unfunctionalized hydrocarbon substrates, among other reactions. The adaptation of P450 enzymes to a particular oxidative reaction involving alkenes is of great interest for the design of new synthetically useful biocatalysts. However, the mechanism that these enzymes utilize to precisely modulate the chemoselectivity and distinguishing between competing alkene double bond epoxidations and allylic C-H hydroxylations is sometimes not clear, which hampers the rational design of specific biocatalysts. In a previous work, a P450 from Labrenzia aggregata (P450LA1) was engineered in the laboratory using directed evolution to catalyze the direct oxidation of trans-ß-methylstyrene to phenylacetone. The final variant, KS, was able to overcome the intrinsic preference for alkene epoxidation to directly generate a ketone product via the formation of a highly reactive carbocation intermediate. Here, additional library screening along this evolutionary lineage permitted to serendipitously detect a mutation that overcomes epoxidation and carbonyl formation by exhibiting a large selectivity of 94 % towards allylic C-H hydroxylation. A multiscalar computational methodology was applied to reveal the molecular basis towards this hydroxylation preference. Enzyme modelling suggests that introduction of a bulky substitution dramatically changes the accessible conformations of the substrate in the active site, thus modifying the enzymatic selectivity towards terminal hydroxylation and avoiding the competing epoxidation pathway, which is sterically hindered.

Copper-catalysed dehydrogenation or lactonization of C(sp3)-H bonds.

Cytochrome P450 enzymes are known to catalyse bimodal oxidation of aliphatic acids via radical intermediates, which partition between pathways of hydroxylation and desaturation1,2. Developing analogous catalytic systems for remote C-H functionalization remains a significant challenge3-5. Here, we report the development of Cu(I)-catalysed bimodal dehydrogenation/lactonization reactions of synthetically common N-methoxyamides through radical abstractions of the γ-aliphatic C-H bonds. The feasibility of switching from dehydrogenation to lactonization is also demonstrated by altering reaction conditions. The use of a readily available amide as both radical precursor and internal oxidant allows for the development of redox-neutral C-H functionalization reactions with methanol as the sole side product. These C-H functionalization reactions using a Cu(I) catalyst with loading as low as 0.5 mol.% is applied to the diversification of a wide range of aliphatic acids including drug molecules and natural products. The exceptional compatibility of this catalytic system with a wide range of oxidatively sensitive functionality demonstrates the unique advantage of using a simple amide substrate as a mild internal oxidant.

A fungal P450 enzyme from Fusarium equiseti HG18 with 7ß-hydroxylase activity in biosynthesis of ursodeoxycholic acid.

Cytochrome P450 enzyme with 7ß-hydroxylation capacity has attracted widespread attentions due to the vital roles in the biosynthesis of ursodeoxycholic acid (UDCA), a naturally active molecule for the treatment of liver and gallbladder diseases. In this study, a novel P450 hydroxylase (P450FE) was screen out from Fusarium equiseti HG18 and identified by a combination of genome and transcriptome sequencing, as well as heterologous expression in Pichia pastoris. The biotransformation of lithocholic acid (LCA) by whole cells of recombinant Pichia pastoris further confirmed the C7ß-hydroxylation with 5.2% UDCA yield. It was firstly identified a fungal P450 enzyme from Fusarium equiseti HG18 with the capacity to catalyze the LCA oxidation producing UDCA. The integration of homology modeling and molecular docking discovered the substrate binding to active pockets, and the key amino acids in active center were validated by site-directed mutagenesis, and revealed that Q112, V362 and L363 were the pivotal residues of P450FE in regulating the activity and selectivity of 7ß-hydroxylation. Specifically, V362I mutation exhibited 2.6-fold higher levels of UDCA and higher stereospecificity than wild-type P450FE. This advance provided guidance for improving the catalytic efficiency and selectivity of P450FE in LCA hydroxylation, indicative of the great potential in green synthesis of UDCA from biologically toxic LCA.

A Cytochrome P450 Catalyzes Oxidative Coupling Formation of Insecticidal Dimeric Indole Piperazine Alkaloids.

Cytochrome P450 (CYP450)-catalyzed oxidative coupling is an efficient strategy for using simple building blocks to construct complex structural scaffolds of natural products. Among them, heterodimeric coupling between two different monomers is relatively scarce, and the corresponding CYP450s are largely undiscovered. In this study, we discovered a fungal CYP450 (CpsD) and its associated cps cluster from 37208 CYP450s of Pfam PF00067 family member database and subsequently identified a group of new skeleton indole piperazine alkaloids (campesines A-G) by combination of genome mining and heterologous synthesis. Importantly, CYP450 CpsD mainly catalyzes intermolecular oxidative heterocoupling of two different indole piperazine monomers to generate an unexpected 6/5/6/6/6/6/5/6 eight-ring scaffold through the formation of one C-C bond and two C-N bonds, illuminating its first dimerase role in this family of natural products. The proposed catalytic mechanism of CpsD was deeply investigated by diversified substrate derivatization. Moreover, dimeric campesine G shows good insecticidal activity against the global honeybee pest Galleria mellonella. Our study shows a representative example of discovering new skeleton monomeric and dimeric indole piperazine alkaloids from microbial resources, expands our knowledge of bond formation by CYP450s and supports further development of the newly discovered and engineered campesine family compounds as potential biopesticides.

An Enzymatic Carbon-Carbon Bond Cleavage and Aldol Reaction Cascade Converts an Angular Scaffold into the Linear Tetracyclic Core of Ochraceopones.

Meroterpenoids of the ochraceopones family featuring a linear tetracyclic scaffold exhibit exceptional antiviral and anti-inflammatory activities. The biosynthetic pathway and chemical logic to generate this linear tetracycle, however, remain unknown. In this study, we identified and characterized all biosynthetic enzymes to afford ochraceopones and elucidated the complete biosynthetic pathway. We demonstrated that the linear tetracyclic scaffold of ochraceopones was derived from an angular tetracyclic precursor. A multifunctional cytochrome P450 OchH was validated to catalyze the free-radical-initiated carbon-carbon bond cleavage of the angular tetracycle. Then, a new carbon-carbon bond was verified to be constructed using a new aldolase OchL, which catalyzes an intramolecular aldol reaction to form the linear tetracycle. This carbon-carbon bond fragmentation and aldol reaction cascade features an unprecedented strategy for converting a common angular tetracycle to a distinctive linear tetracyclic scaffold in meroterpenoid biosynthesis.

Biosynthesis of a new skyllamycin in Streptomyces nodosus: a cytochrome P450 forms an epoxide in the cinnamoyl chain.

Activation of a silent gene cluster in Streptomyces nodosus leads to synthesis of a cinnamoyl-containing non-ribosomal peptide (CCNP) that is related to skyllamycins. This novel CCNP was isolated and its structure was interrogated using mass spectrometry and nuclear magnetic resonance spectroscopy. The isolated compound is an oxidised skyllamycin A in which an additional oxygen atom is incorporated in the cinnamoyl side-chain in the form of an epoxide. The gene for the epoxide-forming cytochrome P450 was identified by targeted disruption. The enzyme was overproduced in Escherichia coli and a 1.43 Å high-resolution crystal structure was determined. This is the first crystal structure for a P450 that forms an epoxide in a substituted cinnamoyl chain of a lipopeptide. These results confirm the proposed functions of P450s encoded by biosynthetic gene clusters for other epoxidized CCNPs and will assist investigation of how epoxide stereochemistry is determined in these natural products.

Multiscale Study on the Intramolecular C-S Bond Formation Catalyzed by P450 Monooxygenase CxnD Involved in the Biosynthesis of Chuangxinmycin: The Critical Roles of Noncrystal Water Molecule and Conformational Change.

Cytochrome P450 monooxygenase CxnD catalyzes intramolecular C-S bond formation in the biosynthesis of chuangxinmycin, which is representative of the synthesis of sulfur-containing natural heterocyclic compounds. The intramolecular cyclization usually requires the activation of two reaction sites and a large conformational change; thus, illuminating its detailed reaction mechanism remains challengeable. Here, the reaction pathway of CxnD-catalyzed C-S bond formation was clarified by a series of calculations, including Gaussian accelerated molecular dynamics simulations and quantum mechanical-molecular mechanical calculations. Our results revealed that the C-S formation follows a diradical coupling mechanism. CxnD first employs Cpd I to abstract the hydrogen atom from the imino group of the indole ring, and then, the resulted Cpd II further extracts another hydrogen atom from the thiol group of the side chain to afford a diradical intermediate, in which a noncrystal water molecule entering into the active site after the formation of Cpd I was proved to play an indispensable role. Moreover, the diradical intermediate cannot directly perform the coupling reaction. It should first undergo a series of conformational changes leading to the proximity of two reaction sites. It is the flexibility of the active site of the enzyme and the side chain of the substrate that makes the diradical coupling to be successful.

Enzymatic Assembly of Diverse Lactone Structures: An Intramolecular C-H Functionalization Strategy.

Lactones are cyclic esters with extensive applications in materials science, medicinal chemistry, and the food and perfume industries. Nature's strategy for the synthesis of many lactones found in natural products always relies on a single type of retrosynthetic strategy, a C-O bond disconnection. Here, we describe a set of laboratory-engineered enzymes that use a new-to-nature C-C bond-forming strategy to assemble diverse lactone structures. These engineered "carbene transferases" catalyze intramolecular carbene insertions into benzylic or allylic C-H bonds, which allow for the synthesis of lactones with different ring sizes and ring scaffolds from simple starting materials. Starting from a serine-ligated cytochrome P450 variant previously engineered for other carbene-transfer activities, directed evolution generated a variant P411-LAS-5247, which exhibits a high activity for constructing a five-membered ε-lactone, lactam, and cyclic ketone products (up to 5600 total turnovers (TTN) and >99% enantiomeric excess (ee)). Further engineering led to variants P411-LAS-5249 and P411-LAS-5264, which deliver six-membered δ-lactones and seven-membered ε-lactones, respectively, overcoming the thermodynamically unfavorable ring strain associated with these products compared to the γ-lactones. This new carbene-transfer activity was further extended to the synthesis of complex lactone scaffolds based on fused, bridged, and spiro rings. The enzymatic platform developed here complements natural biosynthetic strategies for lactone assembly and expands the structural diversity of lactones accessible through C-H functionalization.

Characterization and heterologous reconstitution of Taxus biosynthetic enzymes leading to baccatin III.

Paclitaxel is a well known anticancer compound. Its biosynthesis involves the formation of a highly functionalized diterpenoid core skeleton (baccatin III) and the subsequent assembly of a phenylisoserinoyl side chain. Despite intensive investigation for half a century, the complete biosynthetic pathway of baccatin III remains unknown. In this work, we identified a bifunctional cytochrome P450 enzyme [taxane oxetanase 1 (TOT1)] in Taxus mairei that catalyzes an oxidative rearrangement in paclitaxel oxetane formation, which represents a previously unknown enzyme mechanism for oxetane ring formation. We created a screening strategy based on the taxusin biosynthesis pathway and uncovered the enzyme responsible for the taxane oxidation of the C9 position (T9αH1). Finally, we artificially reconstituted a biosynthetic pathway for the production of baccatin III in tobacco.

Dimer Stabilization by SpyTag/SpyCatcher Coupling of the Reductase Domains of a Chimeric P450 BM3 Monooxygenase from Bacillus spp. Improves its Stability, Activity, and Purification.

The vast majority of known enzymes exist as oligomers, which often gives them high catalytic performance but at the same time imposes constraints on structural conformations and environmental conditions. An example of an enzyme with a complex architecture is the P450 BM3 monooxygenase CYP102A1 from Bacillus megaterium. Only active as a dimer, it is highly sensitive to dilution or common immobilization techniques. In this study, we engineered a thermostable P450BM3 chimera consisting of the heme domain of a CYP102A1 variant and the reductase domain of the homologous CYP102A3. The dimerization of the hybrid was even weaker compared to the corresponding CYP102A1 variant. To create a stable dimer, we covalently coupled the C-termini of two monomers of the chimera via SpyTag003/SpyCatcher003 interaction. As a result, purification, thermostability, pH stability, and catalytic activity were improved. Via a bioorthogonal two-step affinity purification, we obtained high purity (94 %) of the dimer-stabilized variant being robust against heme depletion. Long-term stability was increased with a half-life of over 2 months at 20 °C and 80-90 % residual activity after 2 months at 5 °C. Most catalytic features were retained with even an enhancement of the overall activity by ~2-fold compared to the P450BM3 chimera without SpyTag003/SpyCatcher003.

Computational identification and analysis of deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) in the human POR gene: a structural and functional impact.

Cytochrome P450 oxidoreductase (POR) protein is essential for steroidogenesis, and POR gene mutations are frequently associated with P450 Oxidoreductase Deficiency (PORD), a disorder of hormone production. To our knowledge, no previous attempt has been made to identify and analyze the deleterious/pathogenic non-synonymous single nucleotide polymorphisms (nsSNPs) in the human POR gene through an extensive computational approach. Computational algorithms and tools were employed to identify, characterize, and validate the pathogenic SNPs associated with certain diseases. To begin with, all the high-confidence SNPs were collected, and their structural and functional impacts on the protein structures were explored. The results of various in silico analyses affirm that the A287P and R457H variants of POR could destabilize the interactions between the amino acids and the hydrogen bond networks, resulting in functional deviations of POR. The literature study further confirms that the pathogenic mutations (A287P and R457H) are associated with the onset of PORD. Molecular dynamics simulations (MDS) and essential dynamics (ED) studies characterized the structural consequences of prioritized deleterious mutations, representing the structural destabilization that might disrupt POR biological function. The identified deleterious mutations at the cofactor's binding domains might interfere with the essential interactions between the protein and cofactors, thus inhibiting POR catalytic activity. The consolidated insights from the computational analyses can be used to predict potential deleterious mutants and understand the disease's pathological basis and the molecular mechanism of drug metabolism for the application of personalized medication. HIGHLIGHTSNADPH cytochrome P450 oxidoreductase (POR) mutations are associated with a broad spectrum of human diseasesIdentified and analyzed the most deleterious nsSNPs of POR through the sequence and structure-based prediction toolsInvestigated the structural and functional impacts of the most significant mutations (A287P and R457H) associated with PORDMolecular dynamics and PCA-based FEL analysis were utilized to probe the mutation-induced structural alterations in PORCommunicated by Ramaswamy H. Sarma.

Microbial P450 repertoire (P450ome) and its application feasibility in pharmaceutical industry, chemical industry, and environmental protection.

Cytochrome P450s (CYPs) are heme-thiolated enzymes that catalyze the oxidation of C-H bonds in a regio- and stereo-selective manner. CYPs are widely present in the biological world. With the completion of more biological genome sequencing, the number and types of P450 enzymes have increased rapidly. P450 in microorganisms is easy to clone and express, rich in catalytic types, and strong in substrate adaptability, which has good application potential. Although the number of P450 enzymes found in microorganisms is huge, the function of most of the microorganism P450s has not been studied, and it contains a large number of excellent biocatalysts to be developed. This review is based on the P450 groups in microorganisms. First, it reviews the distribution of P450 groups in different microbial species, and then studies the application of microbial P450 enzymes in the pharmaceutical industry, chemical industry and environmental pollutant treatment in recent years. And focused on the application fields of P450 enzymes of different families to guide the selection of suitable P450s from the huge P450 library. In view of the current shortcomings of microbial P450 in the application process, the final solution is the most likely to assist the application of P450 enzymes in large-scale, that is, whole cell transformation combined with engineering, fusion P450 combined with immobilization technology.

Enhancing ferryl accumulation in H2O2-dependent cytochrome P450s.

A facile strategy is presented to enhance the accumulation of ferryl (iron(IV)-oxo) species in H2O2 dependent cytochrome P450s (CYPs) of the CYP152 family. We report the characterization of a highly chemoselective CYP decarboxylase from Staphylococcus aureus (OleTSA) that is soluble at high concentrations. Examination of OleTSA Compound I (CpdI) accumulation with a variety of fatty acid substrates reveals a dependence on resting spin-state equilibrium. Alteration of this equilibrium through targeted mutagenesis of the proximal pocket favors the high-spin form, and as a result, enhances Cpd-I accumulation to nearly stoichiometric yields.